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ABSTRACT 

This report presents the' results and conclusions of a series of experi- 
ments conducted for a Self-Contained Navigation System (SCNS) carried aboard 
the Applications Technology Satellite, ATS-III. The system, consisting of a 
flight model star sensor and auxiliary ground based data processing equipment, 
utilized star field mapping techniques to determine satellite attitude. In 
addition sightings made on both an ejected sphere and the moon were used to 
determine satellite orbital parameters. 

The ATS-III satellite was launched into synchronous orbit on 6 November 
1967. The SCNS equipment was initially operated on 2 March 1968. During the 
period between 13 March and 4 October 1968, 10 separate experimental opera- 
tions were conducted. 

Experiments were operated successfully for the determination of attitude. 
With the sun 45 degrees off the optical axis, useful data was obtained in full 
daylight when the scanning motion of the satellite directed the optical axis 
away from the sun. Scattered light from the sun masked star transit observa- 
tions during one-half of the scan although the photomultiplier tube recovery 
rate permitted detection of stars as dim as visual magnitude +3.30 over the 
other half of the 0.6 second scan period. 

Automatic reduction of the telemetered data, including star identifi- 
cation and attitude computation was accomplished. The standard deviation of 
the residuals for individual star transits was observed to be typically 60 arc 
seconds. The standard deviation of the derived attitude quantities varied 
within the range of 14 to 20 arc seconds for cases of good geometry. Position 
of star images within the slit were shown to be interpolative to a small 
fraction of the 7.6 arc minute slit width. 

Difficulties were experienced in obtaining the precision required to 
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calculate orbital parameter from sightings of an ejected sphere. However, 
moon sightings were used as a substitute in conducting the navigation experi- 
ment. Satellite orbital position computed from moon transits was determined 
within an error band from -50 to + 40 kilometers over a 28 day interval when 
compared to a moon ephemeris . 
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FOREWORD 


This report was prepared by the Research Division of Control Data 
Corporation for the National Aeronautics and Space Administration, Goddard 
Space Flight Center, under Contract NAS5-9683. The Control Data Project 
Engineer was William J. Lewis; however, during the data acquisition phase 
of the project Bruce D. Vannelli served as the Project Engineer. Control 
Data Corporation desires to particularly acknowledge the technical guidance 
of Forrest H. Wainscott of the ATS Project Office of the National Aeronautics 
and Space Administration, Technical Project Officer. 

This report has been prepared in four chapters. Chapter 1 is concerned 
with the data acquisition system; Chapter 2, flight equipment performance; 
Chapter 3, determination of satellite attitude; and Chapter 4, determination 
of satellite orbit. 
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INTRODUCTION 


SELF-CONTAINED NAVIGATION EXPERIMENT 
INTRODUCTION 


The Self-Contained Navigation System experiment had the following 
ob j ectives : 

(1) To determine the three axis attitude of the ATS-III space- 
craft from star transit times as obtained from an on-board 
electro-optical star sensor. Data to be telemetered to a 
ground station for processing by Control Data. 

(2) To determine the position of the spacecraft as a function 
of time by intermittently observing an ejected sphere 
against the background of stars over a period of two or 
more orbits. Also, to use the moon as a substitute for 
or in addition to the ejected spheres to determine the- 
orbital parameters. 


The SONS operation equipment consisted of a flight model star sensor 
mounted within the ATS-III satellite and a ground-based data processing unit 
called a Star Data Formatter. Reduction of the star sensor raw data to pro- 
vide satellite attitude navigational information was accomplished using 
Control Data 3600 and 6600 computers located at Control Data Corporation, 
Minneapolis , Minnesota . 

fl 21 

Earlier documentation 1 * ’ J describes the flight and ground equipment 
and the applicable operational procedures. However, the following brief 
description of the SCNS operation principle is included here as an aid to the 
understanding of this report. 


The flight equipment contains a star sensor and particle ejector. 
Basically, the sensor consists of a lens assembly, slit reticle and photo- 
multiplier detector. Analog electronics process and deliver the photomulti- 
plier output to the satellite telemetry system. 


Radiation from stellar targets contained within the star sensor field 
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of view passes through the sensor aperture and is focussed by an optical 
system on an opaque reticle containing three transparent hair width slits. 

As the satellite rotates, the starlight is gated through the slits onto the 
cathode of the photomultiplier detector which, with its associated electronics, 
produces an analog pulse that is subsequently digitized and transmitted by wide 
band telemetry Digital data, corresponding to ingress and egress times 
for each pulse and necessary identifying remarks, are recorded on magnetic 
tape. The data is subsequently reformatted to be compatible with the FORTRAN 
language of the attitude computer program. This data is used as input to the 
computer program to provide attitude solutions. 

Whereas, star sightings alone are necessary for attitude determination, 
ejected particle sightings provide the additional data required for orbit 
determination. The flight equipment contains three identical ball ejection 
mechanisms for this purpose. 

The Star Data Formatter processes real-time data from the flight unit. 
This equipment routes analog data to an Ampex recorder, converts analog data 
to digital data, monitors data, and transfers digital data through the CTEG 
(Command Test Evaluation Console) to the SDS-910 Computer for preliminary 
processing and storage. Figure 1-1 contains a flow diagram for the data ac- 
quisition of the ATS-III/SCNS experiment. Figure 1-2 is a pictorial repre- 
sentation of the satellite and SCNS equipment. 

The ATS-III satellite was launched into synchronous orbit on 
6 November 1967. The SCNS equipment was initially operated on 2 March 1968. 
During the seven month period between March and October ten separate experi- 
ments were conducted. A summary of operational conditions is included as 
Appendix A. 


1. During the experiment the Data Acquisition Facility at Rosman, North 
Carolina was utilized. 
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Figure I-l: Flow Diagram for Data Acquisition of ATS-III/.SCNS 
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CHAPTER 1 

EXPERIMENT DESCRIPTION 
I. SCNS OPERATIONAL EQUIPMENT 

The SCNS operational flight equipment is shown in Figure 1-1 and con- 
sists of four basic assemblies: the optical system, the detector, particle 

ejection assembly, and electronics assembly. 

A. Optical System 

The optical system is designed to focus star images on its image 
surface (slit plane) which is opague except for the three narrow slits 
shown in Figure 1-2. The light from images focused on any portion of 
these slits is transmitted through the opening to the cathode of the photo- 
multiplier tube. The output voltage of the detector indicates the presence 
of a star image in the slit. Since the satellite rotates, the star images 
move across the slit plane normal to the center slit and at an inclined 
angle to the outer slits. The output signal for each star which crosses 
the slit pattern consists of three pulses with their separation dependent 
upon the zenith position of the star. This is illustrated in Figure 1-3. 


The optical system includes the following items, each of which is 
shown in Figure 1-2. 

Cone: A fiberglass baffle included to prevent stray light from 

illuminating the lens assembly. The baffle configuration was 
constrained to fit within the skin of the spacecraft. As a con- 
sequence, the general shape became that of a frustrum of a cone. 

Objective Lens: An f/0.87 Super-Farron objective lens having a 

30 degree field, a 76-mm focal length and a 40-mm diameter image 
plane . 

The lens assembly was an off-the-shelf item manufactured by the 
Farrand Optical Company. A computer automated ray trace was per- 
formed based upon patent disclosure information. The blur spot 
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Figure 1-3: Pulse Separation for Two Candidate Stars Transiting 

the Slit Pattern 
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(image) diameter for an on-axis source was 5 arc minutes; for 
7.5 degrees off-axis the blur spot diameter was 20 arc minutes 
and for 10 degrees it was 80 arc minutes . 

, • Slit Reticle: A vacuum-deposited opaque layer on a glass sub- 

strate having three translucent slits. The slits were designed 
to cover the entire 30 degree (40-irm) field of the objective and 
form an angle of 10.5 degrees having a vertex located 1,0675 
inches above the optical axis. 

Actual measurements of the slit reticle indicated a value of 
10.454 degrees (+ one arc minute) for the vertex angle and 
1.075 + .003 inches for the vertex-optical axis displacement. 

Based upon actual measurements the slit-to-slit separation angle 
is 11?102 + .02°. The rotational slit width at a point midway 
from the vertex is 7.6 arc minutes. A detailed analysis describ- 
ing slit characterization is included in Appendix D. 

# Radiation shield: An optically flat window of fused silica used 

to protect the objective lens from electron radiation damage. 

# Condenser lens: Consists of two plano-ampheric lenses having a 

40-mm diameter and a 34-mm focal length. 

m Reflector: A 69-rnm diameter, 46-mm focal length aspheric con- 

densor having a reflective aluminum coating on its spherical con- 
vex back surface. 

Sunshade: The sunshade, designed as a cover for the open end of 

the optical system, was removed from the flight equipment prior to 
actual launch. This was done to preclude the possibility of a 
"fail-closed" condition and thus wipe out the experiment. 

The sunshade assembly consisted of a thin mylar shade (configured 
to roll up when released), a release mechanism, and an indicator 
switch. Its intended use was to prohibit direct sunlight from be- 
coming incident upon the cathode during the de-spun launch phase 
of the satellite. Although referred to in References [1] and [2], 
it is not shown in Figure 1-2. 

Data taken during the experiment revealed a loss of sensor sensi- 
tivity to off-axis sources. Relatively bright stars having visual 
magnitudes around 1.5 lying five to ten degrees off the optical 
axis were not regularly detected. This was due principally to 
the severe vignetting of the Super-Farron lens. The effect is 
observed in Figure 1-4 which shows the light transmission char- 
acteristic of the SONS optical system with off-axis position. 
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Figure 1-4: Transmission Characteristic of the SONS Optical System 
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B . Detector 

The detector selected was the model 541-E-01-14 photomultiplier tube 

manufactured by Electro-Mechanical Research, Inc. The PMT was used- with an 

integral high voltage power supply. The detector contains a tri-alkali type 

cathode surface whose quantum efficiency at 4200 Angstroms is typically 0.25. 

Electrons emitted from the cathode are accelerated through a series of dynodes 

which both accelerate and emit electrons to produce a current gain from cath- 

6 7 

ode to collecting anode of 1.5 x 10 to 1.5 x 10 , depending upon voltage 
used. Both gains were employed at different times during the experiment. 

C . ' System Electronics 

The SONS was designed to operate in 'one of two system operating states. 
The first, referred to as System 1, was designed to provide digital informa- 
tion only. The system did this by providing’ a digital -pulse whenever the 
analog star pulse exceeded a fixed detection (or threshold) level. • A step 
attenuator was used to deliver the analog star pulse to the threshold detector, 
which upon ground command, adjusted system gain in discrete steps. It thereby 
served the equivalent function of a variable threshold detector. Additional 
system sensitivity was obtained by adjusting the photomultiplier supply 
voltage . 

System 2 output was analog only, .and essentially reproduced the photo- 
multiplier analog output pulse. Both systems amplified and filtered the 
detector output followed by a filter matched to the average slit dwell time 
of a star. The analog output of System 2 was telemetered to the ground data 
acquisition station where the star pulse information was digitized. The 
ground based data processing equipment circuitry was designed to produce 
digital pulses corresponding to a pre-set threshold crossing of the pulse 
rise and fall. 

During certain phases of the experiment, particularly when the satellite 
passed from the earth's shadow into full sunlight, a significant increase in 
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spatial background noise caused the output to assume a variable baseline 

i , 

characteristic. This resulted in an inability to detect stars using a fixed 
baseline design. Modifications were subsequently .made to the ground based 
data processing equipment to provide a variable threshold capability which 
utilized a baseline following circuit. 

Since it was not possible to modify System 1 (and thereby provide a 

capability to detect stars under daylight conditions), this system.fell irto general 

( 2 ) 

disuse. All of the data obtained subsequent to the initial' operation of 
the SONS equipment was obtained using System 2, modified to include the base- 
line following circuit. 

D . Probe Ejection Mechanism 

The navigation method of. the Self-Contained Navigation System requires 
that a probe be ejected from the satellite at known velocity and direction, 
and positional measurements be taken during the orbit. Optical viewing con- 
straints require that the probe be ejected at a velocity of 3 ft/sec. with an 
elevation angle of -45° south, and an azimuth normal to the satellite velocity 
vector. If the satellite spin velocity is. 100 rpm and the probe release point 
is approximately 31 inches from the spin axis, the resultant tangential velocity 
is 27 ft/sec. The ejector is therefore required to compensate for this velocity 
component. Figure 1-5 shows the geometric interface between satellite and SCNS 
ejection mechanism. 

Several ejector -designs were evaluated, which resulted in the selection 
of the design illustrated in Figure 1-6. The probe, a three inch diameter 
spherical ball, is cupped in the throwing arm prior to and during ejection. 


2. Initial operation occurred during early April 1968 when the satellite 
was in the earth's shadow. The resultant dark sky conditions produced 
a constant (non-varying) baseline. 
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Figure 1-5: Satellite Interface with SCNS Ejection Mechanisir 
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Figure 1-6 
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When the release switch is activated the springs and satellite- centrif- 
ugal force accelerate the ejector to the released position. As illustrated 
in Figure 1-7, the ejector is designed to stop with an angle x = 4.5 degrees 
and with the pivot axis tilted about x axis by t = 4.5 degrees. If the y 
component of the ejector velocity entirely compensates the satellite tan- 
gential velocity, the x and z (parallel to the satellite spin axis) velocities 
are 


v = v. cos t sin y = 2.1 ft/sec. 

X 0 


Cl) 


V 

z 


= v. sin t - 2.1 ft/sec. 
0 


( 2 ) 


and the resultant- speed is 3 ft/sec. with an azimuth angle of zero degrees 
and an elevation angle of 45 degrees. 



Figure 1-7: Ball Ejector Force Diagram 
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II. PRE -FLIGHT CHECK-OUT . 

Prior to installation within the ATS-III satellite, tests were -conducted 
to obtain a measurement .of response of the SCNS flight equipment to a cali- 
brated artificial star and to a star field simulation. 

The artificial star provided a close approximation to an AO-5 star in 
the range of wavelengths where SCNS is most sensitive. The desired intensity 
for a given magnitude was obtained by setting the distance from the artificial 
source aperture to the SCNS aperture. Star intensities were variable to 4th 
magnitude. 'The optical path for the source energy (adjustable to obtain 
desired intensity) reflects' from a rotating mirror used to simulate the 
satellite scan motion. SCNS response to a simulated star at various off-axis 
positions within the field of view was evaluated by varying the vertical posi- 
tion of the source. The on-axis position corresponded to a declination angle of -25 
degrees on the celestial sphere. A more detailed description of the tests 
conducted and their results is contained in reference [2], Figure 1-8 shows 
a schematic diagram for the test configuration used. 

Tests were .also conducted to obtain a measurement of the SCNS response 
to the simulated stars contained in the Star Field Simulator (SFS). The SFS 
was designed, built and calibrated to provide four artificial stars, each 
having a different magnitude and off-axis position. The SFS was fixed to the 
SCNS optical assembly and’ simulated the satellite scanning motion through 
rotation of a mechanical chopping plate which gated the artificial starlight 
onto the SCNS optical system. The SFS also had the capability to adjust star 
intensities within the range of -1.5 to + 3.7 visual magnitude. 

In addition to performing functional check-out and determining character- 
istic behavior for the SCNS flight equipment, the SFS pre-flight calibration 
and testing provided the theoretical (laboratory determined) limits of detect- 
ability for the different operational states possible. Information relative 
to SCNS response versus star magnitude was compiled as a nomograph which was 
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contained as Figure 1-19 of Reference [1]. During flight operation it was 
noted that actual limit of detectability was slightly degraded over the lab- 
oratory determined value. For dark sky conditions the faintest star detected 
was generally 0.1 to 0.4 magnitudes brigher than the theoretical value 
derived from the pre-flight calibration tests. A- summary of sensor detect- 
ability performance for various operational conditions is included as Appen- 
dix B. 
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III. GROUND BASED DATA PROCESSING 

Figure 1-1 illustrates the procedural steps involved in the acquisition 
of data from the ATS-III/SCNS equipment. It indicates the dual path available 
to process the ra-w data in both analog and digital form. During the experi- 
ment operations attitude solutions were obtained utilizing an algorithm named 
SCAT which provides a solution in four unknowns . A second algorithm evolved 
during the course of the project. It contained the capability of solving a 
set of equations containing up to nine unknowns . The program named SCOOT is 
discussed along with program SCAT in Chapter 3, Section I. 

The original Star Data Formatter .(SDF) design was governed by the design 
parameters of the Self-Contained Navigation System, the CTEC/SDS 910 interface 
and computer, and interface to Ampex FR607 tape recorder. The intention was 
to process real time data as it occurred and simultaneously record the analog 
information (System 1 or System -2) on tape. The tape system was considered to 
be a backup system running on a coarser resolution time base compared to the 
Star Data Formatter. The Star Data Formatter was designed to process real 
time data off line (off line being the processing of recorded data from the 
Ampex FR607) . 

The SDF analog electronics was designed to accept video from either 
the SCNS sensor via transmission link or the FR607 data recorder. The SDF 
contained a level detector which provided a positive going output pulse 
during the detected star transit. The leading edge of the pulse corresponded 
to the star-slit ingress time and the trailing edge to the egress time. The 
SDF digital electronics converted the analog pulse to time encoded leading, 
and trailing edges and also tagged the data. The digital electronics also 
recorded the time occurrence of the sun pulse generated by the ATS-III sun 
sensors. Also, the SDF interfaced with the necessary timing signals (1 me, 

1 pps) . Operation of the equipment when the satellite was in the earth's 
shadow permitted an evaluation of sensor performance under dark sky condi- 
tions. A photograph of typical analog output taken while the satellite was 
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in the earth's eclipse is shown in Figure 1-9. The characteristics -of this 
output are a straight baseline, with photomultiplier- and background noise 
superimposed, plus occasional star and dark current pulses.' 

When bright objects, such as the sun, moon, or sunlit earth is in or 
near the field of view, fhq character of the analog output .changes markedly. 
Figure 1-10 shows a typical output for a complete satellite scan with the sun 
illuminating the shield and/or lens for approximately 180 degrees of the scan. 
Large overload and recovery transits are present, and photomultiplier noise is 
increased appreciably in those areas of the scan containing star pulses. The 
straight baseline has disappeared and it is impossible to detect all of the 
desired star pulses with a fixed trigger level. The same situation- exists 
when the moon is in the field of view, as shown in Figure 1-11. 

A level-follower circuit was designed and added to the SDF that measured 
the voltage of the baseline and from it provided a component of the bias input 
to the level-detector. The level-follower had the ability to follow non- 
linear and variable frequency 'excursions of the baseline. An indication of 
the performance of the circuit is shown in Figure 1-12. 

It should be noted when comparing Figure 9 with 10, 11, and 12 that 
Figure 9 has an expanded time scale. The three star pulses shown emerging 
from the straight baseline noise correspond to those appearing near the 
extreme right portion of Figure 10. 



EXPERIMENT DESCRIPTION 



Figure 1-9: Typical Dark Sky Conditions Occurring Within 

Earth's Eclipse Showing Straight Baseline 
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Figure 1-10: Bright Object Interference Showing Varying Baseline 
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Figure 1-11: Lunar Transits Superimposed on Varying Baseline 



Figure 1-12: Output Characteristic of Baseline Following Circuit 
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CHAPTER 2 

FLIGHT EQUIPMENT PERFORMANCE 


INTRODUCTION 

For approximately sixty minutes each day during both spring and fall 
eclipses, the ATS-III satellite was in the earth's shadow. Under those con- 
ditions, with the exception of an illuminated crescent at the earth’s south 
polar region^, the star sensor viewed the nighttime sky. Typical sensor 
output is reproduced in Figure 2-1. 

The data of Figure 2-1 are reproduced portions of a single scan analog 

output at a specific value of sensor gain. ■ A higher gain would produce 

stronger star pulses at the expense of an increased noise level. As expected, 

the brighter stars produce stronger pulses. Sirius, the brightest star in the 

( 2 ) 

sky, has a visual magnitude of -1.58. Antares , the sixteenth brightest^ ' and 
Fomalhaut, the seventeenth brightest, have visual magnitudes of +1.22 and +1.29 
respectively. Both possess comparable pulse heights and both possess weaker 
pulses relative to Sirius. However, it is also observed that a given star 

does not appear to regularly produce an equally strong output pulse as it 

crosses each of the three reticle slits. This can be explained with the aid 
of Figure 2-2. Due to the off-axis attenuation or vignetting of the optical 

system when a star crosses the outer two slits, its radiation is reduced by a 

factor greater than the center slit. The center pulse should, therefore, 
possess a greater height and this is illustrated in Figure 2-2B. However, 


1. The upper edge of the SCNS field of view scanned within 1.5 degrees of 
the earth's limb. A marked rise in background noise level was noted as 
the field scanned past the earth. The halo effect is produced by at- 
mospheric scattering of the sun's rays. 

2. The number designation given a star is predicated upon its visual magni- 
tude, and is ranked higher according to decreasing brightness for all 
stars over the entire celestial sphere. Star number 16 (Antares) is 
the sixteenth brightest star (visual) in the sky. 
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Figure 2-1: Reproduced Portions of Typical' Analog Output for Dark Sky Conditions . 
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SLIT 2 


OPTICAL FIELD 
VIEW 



Stars crossing the outer slits (1, 3) do so at a greater distance from the 
optical axis (a > b). Pulse heights are less than center slit (2) pulse due 
to off-axis attenuation of starlight within the optical system (vignetting). 
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Figure 2-2: Statistical Nature of a Star Pulse 
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due to the statistical behavior of photon, arrival rate the outer slit pulse 
heights could be unequal, as shown in Figure 2-2C. 

The analog data illustrated in Figure 2-1 shows clearly only the 
brighter stars contained within the region of sky swept out by the optical 
field of view . Although not easily discernable, pulses from fainter stars are 
contained within the noise spectrum. Digital filtering techniques implemented 
as routines of the ground based computer automated data processing program 
successfully extracted .and identified stellar targets as faint as visual 
magnitude 3.8^ . 

Figure 2-3 graphically locates those stellar targets which were both 
(4) 

regularly detected ' and occasionally detected. The' grid represents that 
portion of the southern sky swept out as an annulus by the field of view and 
corresponds to a dark sky condition occurring when the satellite was in the 
earth's shadow during the fall eclipse. 

The stellar limiting magnitude (or faintest star detectable) based upon 
detectability alone was, therefore, 3.8. The theoretical limiting magnitude 
is 4.5. However, determination of vehicle attitude was obtained using only 
the eight stellar targets which were regularly detected. Figure 2-3- identi- 
fies these stars and also indicates the excellent geometric spacing of the 
stars over the entire scan - a condition which generally leads to ‘a more, ac- 
curate attitude solution. The limiting magnitude, then, based upon useful 
star information corresponds to the faintest regularly detected .star, Dschubba, 
having a visual magnitude of 2.54. 


3. Criterion was that the star be detected in at least one of the three 
slits 10% of the time. 

4. Criterion was that the star to be detected in all three slits (as trip- 
lets) at least 25% of the time. 
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NO. 

STAR NAME 

MAG 

NO. 

STAR NAME 

MAG. 

1 

SIRIUS 

-1 58 

38 

WEZEN 

1.98 

16 

ANTARES 

1 22 

39 

MIRZAM 

1.99 

17 

FOMALHAUT 

1 29 

48 

NUNK 1 

2.14 

23 

ADHARA 

1 63 

81 

DSCHU8BA 

2.54 


RIGHT ASCENSION (degrees) 
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When other than dark sky conditions prevailed, stellar target detection 
was degraded. These conditions occurred whenever the satellite came out of 

r 

the earth's shadow into' full sunlight and also, but to a lesser extent, when- 
ever the full or nearly full moon entered the field of view. When the satel- 
lite was in full sunlight, 180 to 200 degrees of the scanned star field was 
blanked out by the sun's interference, even though the sun was never closer 
than 28 degrees to the edge of the field of view. When the full moon entered 
the field of view an additional 50 to 70 degrees of scanned star field was 
wiped out. When both conditions occurred, a reduced number of stars were 
regularly detectable. A less accurate attitude solution resulted and in 
certain instances, due to the geometric spacing of detected targets, no solu- 
tion was possible. Figure 2-4 illustrates the effect produced by the combined 
sun and moon interference. In this instance, due to the fact that Fomalhaut 
was regularly detected, attitude solution was possible. The limiting magnitude, 
based upon daylight detectability and using the criterion that the target is 
detected occasionally (less than ten percent of time in any one slit), is 2.91 
visual. If only "usable targets" is assumed as the criterion, the limiting 
magnitude again becomes 2.54. 

During daylight operation of the star sensor a significant rise in the ' 
noise level was observed. Reproduced portions of the analog output for typi- 
cal sunlit conditions are included in Figure 2-5. In addition to the obvious 
increase in noise due to sun and moon interference, the baseline assumes a 
varying characteristic, and contains large overload and recovery transients. 

The straight baseline, typical of dark sky viewing, shown in Figure 2-1 has 
disappeared. It was for this reason that the ground based data processing 
equipment was modified to provide a baseline following circuit described 
earlier which permitted detection of stellar targets not possible with a fixed 
threshold level. 

An additional observation from Figure 2-5 is the characteristic behavior 
of the electronics in recovering from saturation following an incidence of high 
energy radiation. 
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Figure 2-5: Reproduced Portions of Analog Output for Typical Sunlit 

Operation Showing Varying Baseline Characteristic 
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I. STAR SENSOR PERFORMANCE 

A . Detectability 

The primary method of evaluating the detectability of the SONS sensor 
is to plot the detection frequency (the ratio of the number of times a 
particular star/slit was detected to the number of maximum possible detections) 
against the instrument magnitude ^ for a given sensor configuration. Data 
required for the preparation of detection frequency versus instrument magni- 
tude curves are obtained from an analysis of the raw digital data tapes in 
conjunction with the computed satellite attitude. Each star/slit in the 
listing of available stars ordered by visual magnitude to 4.0 is considered 
as a potential candidate able to provide a detectable pulse. Each computed 
star/slit transit time is compared with all of the measured transit times for 
each scan, Whenever a comparison produces agreement to within + 200 micro- 
seconds^, the measured time is associated with the known star/slit. This 
is done for all stars in the field of view. The instrument magnitude which 
depends upon the star color-temperature and its off-axis angle is identified 
with each star candidate. The number of times a star/slit combination is 
theoretically scanned is determined by counting the number of times each 
— degree of right ascension is scanned, deducting from the total those scans 

missed during the computer input/output (tape write) processing time for each 
data set. In general, the brighter on-axis stars are detected with higher 
frequencies than dimmer off-axis stars. Some noteable exceptions are apparent 
but in most cases any apparent detection anomally can be explained. 


5. The instrument magnitude relates stellar magnitudes to the spectral 
response of the SCNS sensor. A discussion of the significance of in- 
strument magnitude and a listing of star corrections is included as 
Appendix C . 

6. At a nominal satellite spin rate of 600 degrees per second, the sensor 
rotational slit width of 7.6 arc minutes is equivalent to a slit transit 
time of 211 microseconds. 
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Figure 2-6 is typical of the. detection frequency versus instrument mag- 
nitude characteristic. The curve was drawn from data taken when the satellite 
was within the earth's shadow. Figure 2-7 was also compiled from data taken 
during dark sky conditions but separated from the data of Figure 2-6 by ap- 
proximately six months. Three observations are made: 


1. There is a great deal of scatter in the data as seen in 
the two figures, a condition which one would expect even 
under the best of detection conditions. However, some 
bright star/slit combinations, noteably star number 1 
(Sirius), in slit 3 appears to be less detectable than 
its predicted frequency. This can be- explained , with 
the aid of Figure 2-8, by the fact that the target is 

so bright that a systematic shift, in the actual detection 
transit time occurs. The cause is fundamental to the. 
electronics. This shift makes the measured and computed 
transit times disagree enough so that occasionally this 
star/slit crossing will fall outside the acceptable 
tolerance about the computed transit time, thereby af- 
fecting a reduced detectability. 

2. In many cases, a dimmer star has a higher detection 

• frequency than its magnitude might predict . This is 
generally caused by a conflict with one or more ad- 
ditional stars. A conflict exists whenever more than 
one star appears in a given slit within a transit time 
interval. For the sake of clarity, most of the points 
below a detection frequency of about .2 have been de- 
leted. 

3. A comparison of the results obtained during spring and 
autumn eclipses reveals an apparent loss of sensitivity 
over the six month time interval involved. This becomes 
apparent when noting the shift to a higher gain setting 
(from 5 to 6) required to detect a given magnitude star. 

Some portion of the apparent loss of sensitivity can be 
accounted for by ground based equipment modifications and 

by some minor variations in the operational methods employed 
at each of the data acquisition intervals. However, it is’ 
apparent that a deterioration of approximately one magnitude 
equivalent sensitivity did occur over six months of oper- 
ational usage. 
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Figure 2-6: Typical Star Detection Frequency Under 

Dark Sky Conditions (Spring Eclipse) 
Gain Setting of 5. 
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Figure 2-7: Typical Star Detection Frequency Under Dark Sky 

Conditions (Autumn Eclipse) Gain Setting of 6 
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- Figure 2-8: Missed Detections of Sirius by Trailing Slit 
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B. Bright Source Recovery 

The recovery characteristic of the star sensor after an encounter with 
very bright sources was a particularly important and unanswered question 
before the ATS-III satellite was launched. Prior to satellite launch, tests 
were conducted relative to photomultiplier response to bright object inter- 
ference. The tests were conducted by exposing the photocathode surface to 
bright objects placed in close proximity. Unfortunately, it was' not possible 
to predict performance with any degree of certainty and, in* particular, be 
able to describe the recovery characteristics of the photomultiplier or its 
behavior when scanning near the sun or moon. 

Discussion of the recovery characteristics of the detector is based 
upon operational data obtained when the satellite was in full sunlight and 
when a nearly full moon was in the sensor field of view. In addition, the 
brightly illuminated halo about the south polar region of the earth was only 
confined to the time period centered about local satellite midnight so that 
the spacecraft, earth, and sun were approximately in line. The relative 
positions of the sun and moon for this operating condition are shown in 
Figure 2-9 along with a representation of the background noise conditions 
expressed as a function of right ascension. The noise pulse frequency results 
were obtained from the original digital data by deleting all of the identified 
stellar detections to visual magnitude 4.0, and subsequently concluding that 
the remaining "detections" are included within the background noise distri- 
bution. The unit of the ordinate scale is expressed as a frequency 
but can be interpreted to mean the number of noise pulses contained within 
each degree of the total scan per scan. For example, a value of 0.1 (on the 
frequency scale) means one noise pulse per degree for every 10 scans . 

The moon's position has been corrected by approximately two degrees for 
parallax to represent apparent right ascension as viewed from the satellite. 

A sun gate data lock-out was operative and was positioned from approximately 

7. Lock-out is defined as the disabling of the ground based data processing 

equipment when energized by pulses from the on-board sun -sensor. 
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Figure 2-9: Typical Noise Pulse Frequency for Daylight Conditions 


FLIGHT EQUIPMENT PERFORMANCE 






FLIGHT EQUIPMENT PERFORMANCE 


( 8 ) 

202 degrees to 23 degrees with the sun at 110 degrees. Under the operat- 
ing condition, the gate had been fully extended to prevent digitizing of the 
■noise pulses- with subsequent overloading of the computer memory and digital 
tape. It is apparent from the data that bright source pulses were processed 
at the extremities of the lock-out gate.. At a later date the sun gate was 
modified to increase -the extent of lock-out to include even greater portions 
of any noise dominated scan.' 

The region. from 236 to 202 degrees is sun associated -noise and could be 
locked out. The large spike at 236 degrees is the sun reflecting off the 
fully extended ejector arm. The burst of noise in the region about 345 to 
324 degrees is noise generated by the moon's -illumination reflecting off the 
side of the sunshield just prior to its entrance into the center of the field 
of view. The large spike at 312 degrees is a very repeatable pulse which 
shows up regularly between the three moon transits. It has not been ade- 
quately explained, although it is conjectured to be a phenomenon related to 
the sensor electronics caused by the high intensity surge of energy from the 
moon^^ , It is also noteworthy to point out that the scan area immediately 
following the moon in right ascension (310 to 280 degrees) is noise free, 
indicating that the photomultiplier tube and electronics are in a recovery 
phase. In fact, this recovery interval seems to extend approximately 42 
degrees* beyond the moon's center as evidenced by the slow rise in background 
noise up to approximately 266 degrees-. 

The positions of five regularly detected stars are also shown in 
Figure 2-9 to provide indication of the difficulty in making positive 
identifications of targets buried within background noise caused by the sun 


8. Since the scan direction was counterclockwise around the North Pole, it 
is shown going right to left on the right ascension scale. 

9. This unexplained spike can be seen lying between the moon/slit 1 and 
moon/slit 2 crossings of Figure 2^5. 
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and moon. No detections of Sirius were made due to the sun gate lock-out but 
other star/slit combinations were detected with the detection frequencies 
indicated in Table. 2-1. The detection frequency is the ratio of the number of 
times a star was detected in the slit to the maximum number possible. 

Star 48 (Nunki) - appears in Table 2-1 in slit 3 only and at an extremely 
low detection frequency. This is due to the 'proximity of Star 48 to the moon. 
Slit 3 is slightly farther (-28?9) from the moon than slits 2 and 1 (-24?5 and 
-20? 1 respectively) which, of course, would be obliterated by the moon's il- 
luminance . . 


The only star/slit combination that appears out of place in Table 2-1 is 
16/2; however, there is some rationale for its low position. For the parti- 
cular set of operating conditions at that time , the spacecraft sun pulse which takes 
priority over star transits , occurred somewhere between 16/1 and 16/3 
transits. This occasionally caused the loss of slit 2 transits. 

A further illustration of the effect of the moon on star detection 
frequency is shown in Figure 2-10 in which a detection frequency paramater *S 
is identified as a function of moon separation angle. In an attempt to remove 
some of the minor detection frequency variations caused by changes in the 
threshold level from one data set to the next and to account for the influence 
of instrument magnitude on detection frequency, a normalized parameter *S was 


chosen where 


oO) 

*s<W = !»__ 

n Cm) 

■ S 23 


o( m ) 

S 23 - 
g (263) 


and 


s (m) 

n 


is the detection frequency of the nth star on day m. 


10 ‘ 


A condition established within the ground based data processing equip- 
ment . 
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’ . TABLE 2-1 

TOPICAL STAR/SLIT DETECTION FREQUENCIES 
FOR DAYLIGHT CONDITIONS 


Scan 

From 



Star/Slit 

-56.4 

16/1 

36.2 

17/2 

-65-4 

16/3 

30.9 

17/3 

41.5- 

17/1 

-68.2 

81/2 

-71.7 

81/3 

-60.9 

16/2 

-64.7. . 

81/1 

-28.9 

‘ 48/3 


Detection Frequency 

.884 

.873 

.784 

.747 

.646 

.571 

.398 

.313 

.064 

.013 


Instrument Magnitude 

1.72 
1.39 
1.72 
1.60 
• 1.60 
2.30 
2.41 
1.56 
2.41 
2.09 


Normalized data is presented in Figure 2-10 using a particular day 


(Day 263) and the star/slit combination 23/2 as the basis for normalization. 


Day 263 was used as a basis for normalization to account for variations 
in detectability due to the moon's presence because on that day the moon was 
well outside the field of view. The curves give indication that for stars 
well away from the moon, the detection frequency of most stars is increased by 
the moon's presence in the field of view. This is not intuitive and may be 
related to the characteristic nature of the baseline following circuits. 

The detection frequency of star 23 was used for normalization to minimize 
threshold variations because it- has- an .instrument magnitude comparable to that 
for stars 16 and 17 (and thus the sensor should respond to these stars in much 
the same manner). More importantly, it was well removed from the moon's 


11. Scan separation in degrees is positive when the star appears in front of 
the moon (in the direction of the scan), negative when the star occurs 
after the moon. 
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influence for each of the data acquisition periods employed. Three of the 
stars, 16, 48 and 81, were observed on several nights when the moon was both 
behind (negative scan separations) and ahead (positive scan separations) of 
the stars. Using these and their ranges of separations from the moon just 
prior to and shortly after encountering 5 the moon, it is shown that the detec- 
tion cutoff point is in the range of from 0 to 16 degrees in front of the 
center of ’the moon's position. The detection turn-on point is located from 
20 to 45 degrees behind the moon's center. As would be expected, the recovery 
period extends over a broader range than the more abrupt turn-off region. By 
deleting the moon-no moon’ normalization factor S^ 63 ^/S a slower recovery 

(after lunar encounter) to normal detectability for the two fainter stars 48 
and 81 is evident. This is shown for the data normalized using 23/2 only. 

C . Spatial Noise 

A measure of SCNS instrument performance as a star detection device is 
obtainable by comparing the star detection frequencies described in Section A 
with the noise background from which such detections are -determined . Figure 
2-11 shows the noise pulse frequency under typical dark sky conditions just 
as Figure 2-9 did for daylight conditions. The noise pulse frequency histo- 
grams indicates the position of the sun as observed from the satellite. 

Because the results shown were obtained from the eclipse data, the sensor 
field of view does not scan past the sun but sees only the halo about the 
earth. The earth directly blocks any incident sunlight from illuminating the 
satellite. The geometry of the sun-earth-satellite relationship is given in 
Figure 2-12. It is apparent, since no direct sunlight illuminates the satel- 
lite, that the high noise pulse frequencies observed in Figure 2-11 attributed 
to solar Illumination of the earth's atmosphere “which appears as a halo around 
the earth in the region of the poles.' 

A second observation which may be made from a study of the background 
noise is the increase in noise pulses per scan in the region of the Milky Way, 
Notice that except for star 17 (Fomalhaut) other regularly detectable stars, 

1, 16, 48 and 81, all fall within the Milky Way. 
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Figure 2-11: Typical Noise Pulse Frequency for Dark Sky Conditions 
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Figure 2-12: Geometric Relationship Occurring During Eclipse 

Operations 


U. Measurement Accurracy 

Transit time accuracy of star detections for the SONS sensor have been 
analyzed, utilizing three stars having widely separated magnitudes. The 
selected stars and their characteristics are listed in Table 2-II. Each was 
regularly detected in slit 2. 


TABLE 2 -II 

SELECTED STARS USED IN DETECTION ACCURACY ANALYSIS 


Star Number 
1 

16 

48 


Star Name 

Sirius 
Ant ares 
Nunki 


Visual 

Magnitude 

-1.58 

1.22 

2.14 


Instrument 

Magnitude 

- 0.8 

1.6 

2.0 
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Figure 2-13 indicates the transit time position -within a particular scan 
for each of the stars listed in Table 2-II for dark sky conditions occurring 
in both spring and autumn eclipses and for similar sensor operating conditions 
Broken lines are used to connect points where detections - from intermediate 
scans are missing. The solid vertical line for each star/slit combination 
represents the averaged transit time (mean). Also indicated for each graph 
is the standard deviation (cr) of the transits around the mean. Transit times 
were obtained by averaging the leading and trailing edges of each pulse (in 
effect adding one-half the pulse width to the leading edge threshold crossing 
time) . % 

One technique developed in conjunction with the project activity was a 
method of filtering star transits on the basis of both position within the 
scan and its pulse width. Transits were discarded which deviated by more than 
one standard deviation from a mean transit time and which had pulse widths 
exceeding a predetermined narrow band-width. Filtering was done for all the 
regularly detectable stars in order to produce a set of highly "selective 
filtered transit times. When these times are used as input data to the atti- 
tude algorithm, a. reliable and accurate solution with small functional residu- 
als results. 

E . Lunar Detectability 

When the full moon which subtends slightly over 0.5° is centered on one 

of the slits (each having a 7.6 arc minute projected width) the light input to 

4- 

the photomultiplier is approximately 10 times that from Sirius, the brightest 
star in the. sky. When the moon is detected the photomultiplier output elec- 
tronics are severely overdriven. Recovery of the electronics from this con- 
dition exhibits a complex pattern. This is observable- in Figure 2-5. The 
net effect of this behavior, as it relates to the interpretation of the lunar 
transit data, is that the leading edges of the apparent moon pulses are actual 
ly recovery transits from the preceding pulses. The pulses result from moon 
crossings of a preceding slit and/or moon reflections off the cone. 
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Figure 2-13: Star Measurement Accuracies for Spring and Autumn Eclipse Conditions 
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Moon simulated laboratory tests were performed using the spare flight 
SONS model. The results led to the conclusion that the trailing edges of the 
lunar pulses correspond to the optical event associated with the instant of 
moon/slit egress. This was corroborated by analyzing over thirty separate 
data sets containing moon transits. The expected average lunar trailing edge 
detection was calculated for each data set and for each slit, taking into 
account the parallax in the moon's apparent direction and extent, as seen 
from the spacecraft. Because the moon is an extended source and the slit 
crossings occur at various angles, each slit exits from different points on 
the moon's illuminated surface. This will give rise to a small differential 
in the slit-to-slit separation times. For the days under consideration, the 
slits exited from the moon's limb, not the terminator. Thus, lunar phase 
does not affect the lunar detection problem. 

f41 f41 

Accurate ephemerids for the position of the sun L J , the moon L J , and 
the spacecraft^* and geocentric semi -diameter of the moon^^ were fitted to 
polynominal interpolation expressions covering the time of the SONS data ac- 
quisition. Then, for the time of each data set, the apparent celestial 
direction of the. center of the moon as seen from the spacecraft was determin- 
ed. The precise spacecraft attitude solution (including values for the 
second and. third time derivative of the azimuth angle s and s’ for the group 
of data sets under consideration was then used to calculate the time when the 
trailing edge of each slit would cross the moon's limb. 

If n denotes the -time varying normal to the ith slit plane; m the ap- 
parent celestial direction of the moon as seen from the spacecraft; cp the 
apparent lunar radius as seen from the spacecraft; and £ the half -width of 
the slit, the condition that the moon begins its egress from the trailing 
edge of slit i is describable by the equation: 

m • n^ = cos ^90° + (cp + £)J. (I) 

Since this equation is a single condition in t^, the time of the trailing 
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edge crossing, it may be 'solved to determine t^. 

For each data set, each of the 15 seconds duration, the average lunar 
trailing edge transit times -were computed for the measured data. These 
measured detection times were then compared with the predicted detection 
times. If; for each slit, these differences were constant from set-to-set 
and from day-to-day, they could be used' as correction factors to be added to 
the measured detection times to obtain the actual trailing edge detection 
times . 

The results of this analysis are shown in Figure 2-14 where, for each 
day, the average of the difference between the predicted and measured detec- 
tion times is plotted against the phase of the moon as seen from the space- 
craft. The standard deviations of the points range from one to seven micro- 
seconds with an average of four microseconds . The electronics of the space- 
craft are overloaded in such a way that the trailing edge detection by the 
first slit is not as consistent as the trailing edge detection by slits 2 and 
3. For five data sets the correction factors for slit 1 turned out to be 
very large (i.e., .over 500 msec). These results were not included in the 
data presented. 

There appears to be a strong correlation between the correction factor 

and the phase of the moon. The correction is larger at full moon (0 degrees 

o 

phase) than at new moon (180 phase). This effect may be qualitatively ex- 
plained as follows-: The output of the photomultiplier tube is accumulated 

in energy storing .elements within the electronics. When the input voltage 
falls below a given level, the energy balance in the first stages- can no 
longer be maintained, and the relative input polarity to the last stage 
reverses. When this occurs, the output state changes almost d is continuously 
from the saturation state of maximum illumination to the opposite saturation 
state which results in zero output. The trailing edge of the pulse is then 
recorded. When the moon is crescent, the integrated input from the 
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photomultiplier saturates the electronics for a shorter time than when the 
moon is full and the input continues to be integrated by the energy storage 
elements while the pulse goes through its peak and starts down its trailing 
edge. Only when the input voltage falls below the level needed to maintain 
the energy balance do the polarity reversal and output state changes occur. 
However, when the moon is full, the energy storage elements are in equili- 
brium, so that at the instant the early stages come out of saturation, the 
polarity reversal occurs. Hence, the trailing edge detections at full moon 
should occur earlier than the trailing edge detections when the moon is 
crescent. Thus, the correction factor, which is the difference between the 
predicted and measured detection times, should be larger at full moon than 
when the moon is crescent. Such behavior is, indeed, observed. 

The simplest relationship between the phase and the correction 
factor is a linear relation between the two. The effect on the solution of 
the navigation problem of applying such a correction factor to each measure- 
ment of the trailing edge of the moon is described in Chapter 4. 

F. Ejected Probe Detectability 

Two probe ejections were conducted with neither providing adequate 
data to solve the navigation problem. It was concluded that the first 
launch exhibited an error in timing but near nominal mechanical performance; 
in the second launch, timing was correct but the ejection mechanism malfunc- 
tioned. 

Even though the malfunction prevented probe sightings from being suc- 

( 12 ) 

cessfully utilized to determine satellite orbital position, a measure of 
the detectability of the sensor-probe combination was obtained. The second 
probe was tracked by the sensor to a limit of 8.6' miles, equivalent to a star 
of 1.92, visual, and 2.88, instrument magnitude. 


12. Launch error for the second probe prevented the probe from re-entering 
the sensor field of view on its return orbit 24 hours after launch. 
Sightings were only possible up to 2.5 hours after .launch. 
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The parallax in the SCNS sightings of the probe on the first few scans 
after ejection was used to determine the initial conditions of launch. The 
specific conditions are the time of release, t Q , the azimuth of release \ 
relative to a radial line through the spin axis and the release point, the 
elevation of release E relative to the azimuth plane, and the speed of ejec- 
tion Av relative to the spacecraft. 

For the second probe, the parallax analysis has been carried through 
and the results are summarized in Table 2-III. 

TABLE 2 -III 

INITIAL CONDITIONS FOR SECOND PROBE EJECTION - 
RESULTS FROM PARALLAX ANALYSIS 


Time of release (20 July 1968): 
Relative azimuth of release: 
Elevation of release: 

Release speed: 


t = Oh 0m 0.7626 seconds 
o 

\ = -12?89 

E = 23?02 

Av =4.54 feet/second 


It should be noted here that the navigation equations described in Chapter 4 
could not be solved for a precise value of the release speed, particularly 
for the short time base - two and one-half hours - over which probe data was 
acquired. Thus, the parallax analysis was the only method which yielded a 
direct determination of Av. 

Once the release speed was determined, it was used in connection with 
the navigation equations to describe the apparent direction and separation 
distance from the satellite to the probe as a function of time.' The results 
of this characterization are summarized in Table 2-IV for three times at 
which probe detections were acquired on 20 July (Day 202). 
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TABLE 2 -IV 

SECOND EJECTED PROBE - .* 
DIRECTION AND SEPARATION '.FROM THE SATELLITE 


Time - 

Right Ascension 

Declination 

Separation 

1:48 

261?85 

-20?51 

5.9 miles 

2:15 

262?65 

-19?44 

7.6 miles 

2:30 

263?29 

. -18?79 • 

8.6 miles 


The data in the last column of Table 2-IV represent one of the factors 

in a separation - apparent magnitude relation for the probe; such a relation 

provides the basic measure of probe detectability. 

To complete the detectability characterization, the apparent magnitudes 
corresponding to the separations were found. Actually, since the -probe 
brightness follows an inverse square relationship with distance, it is suffi- 
cient to determine the apparent magnitude corresponding to one known separa- 
tion. At 1:48, the detection frequency of the probe in each slit was slightl 
below the detection frequencies of the two brightest accessible stars in the 
field (Fomalhaut and Antares) but still well above the noise level. At 1:48 
the instrument magnitude of the probe, corrected to the- optical axis, was 
found to be 1-.73 with a spread of only 0.09 magnitude among the three slits. 

The instrument magnitude corrected to the optical axis takes the par- 
ticular spectral response characteristics of the SONS instrument into account 
If it is assumed that the probe reflectivity' is essentially constant with 
wavelength, then the color index of the probe was essentially that of the 
sun. Since the sun is a rather red star, as seen by SONS, the instrument 
magnitude for a star of this color (B-V = 0.62) is 0.38 units fainter than it 
visual magnitude. Thus, at 1:48 the visual magnitude of the probe is 1.35. 

As was noted above, the probe brightness can be predicted for other 
separation values by employing the inverse square law. Table 2-V presents 
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this information in equivalent form in terms of time. In the righthand 
column, the magnitude has been corrected for both spectral response and 
location of the probe in the field of view; this is the instrument magnitude 
of the probe as it appears in the center one of the three slits . 


TABLE 2-V 

MAGNITUDES OF SECOND EJECTED PROBE 


Time 

(20 July 1968) 

Visual 

Magnitude 

Instrument Magnitude 
(on the optical axis) 

Apparent Instrument 
Magnitude 

(at the actual off- 
axis crossing for 
slit 2) 

1:48 

1.35 

1.73 

1.89 

2:15 

1.90 

2.28 

2.54 

2:30 

2.17 

2.55 

2.88 


In Table 2 -VI the observed magnitude, as indicated in the last column 
of Table 2-V is compared with the predicted magnitude. The predicted magni- 
tude was obtained before conducting the experiment. 


TABLE 2-VI 

COMPARISON OF OBSERVED AND PREDICTED MAGNITUDES 


Time 

Observed Magnitude 

Predicted Magnitude 

Difference 

1:48 

1.89 

.99 -j 

I 

2:15 

2.54 

1.64 

} 0.90 

2:30 

2.88 

! . 98 

1 


Thus , the predicted magnitude was determined to be less than one magnitude 
difference from the observed magnitude . 
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CHAPTER 3 

DETERMINATION OF SATELLITE ATTITUDE 


INTRODUCTION 

The processing of star data makes use of computer programs designed to 
solve a set of equations ■which describe satellite motion. The development of 
the attitude method presented here, includes a description of the coordinate 
system used and defines the basic constraint equation in the unknowns of the 
problem. Several computer programs to obtain attitude solutions were written 
during the performance of the experiment. An early program designated BOOTS 
was designed to accept hand identified star transits. A fully automatic 
program designated SCAT, needing no hand identification, was subsequently 
written which solves for four parametric unknowns. A later program was de- 
signed to accommodate an increased number of unknowns . Each of these pro- 
grams is discussed in some detail within the chapter. The SCAT program which 
served as the basic program in determining satellite attitudes has each of 
its routines fully described in Section I B. This description has been 
.presented in sufficient detail to serve as a users manual. Section II docu- 
ments the numerical results and for the most part tabulates the values of 
each of the solved for unknowns according to the environmental conditions of 
each experiment. A measure of instrument accuracy is expressed in the tabu- 
lation of the residuals for each of the unknowns , and are also included in 
the section. 
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CONCLUSIONS 

Chapter 3 relates the results of the' data -processing methods and the 
accuracy of the attitude solutions. The most significant observations from 
the results are included here as conclusions. 


1. Attitude solutions obtained from data taken when the 
satellite was in the earth's eclipse, contained residuals 
which were typically 60 arc seconds. The standard devia- 
tions of the computed spin axis directions' were generally 
less than 10 arc seconds, and for the azimuthal rotation 
about the spin axis less than 15 arc seconds. 

2. Attitude solutions were obtained through the use of a 
completely automatic computer program, requiring no hand 
identification of stellar targets. Two programs were 
written. The basic program. (SCAT) solved for four un- . 
knowns (two pointing angles, azimuthal position and spin 
velocity). A more complete program, based- upon SCAT, was 
used to solve the eight unknown problem. 

3. Attitude solutions obtained -from data taken when the 
satellite was in full sunlight contained residuals which 
were typically 90 arc seconds. The standard deviation 
of the computed spin axis directions were generally 100 
arc seconds and for the azimuthal rotation about the spin 
axis 40 arc seconds. Degraded accuracy was due to the 
reduction in the number of regularly detected targets and ■ 
the resultant poorer geometry. 

4. The parameters which define the location of the slit 
intersection with respect to the satellite, figure axis 
were observed to change with time and particularly with 
each maneuver. Slight changes in spin axis position, most 
probably due to loss of mass, was concluded as the cause. 

5. The degree of satellite nutation (and precession) was 
observed to be negligible in comparison to the ability 
of the star sensor to detect the motion. The coning 
angle was found' to be less than one minute of arc which 
permitted the utilization of a mathematical model which 
assumed rotational motion about a fixed axis. 

6. It was observed that the satellite in passing from full 
sunlight into the earth's shadow increased in spin 
velocity, necessitating the inclusion -of the spin accel- 
eration term as an additional unknown. 
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7. Over the 28 days when the satellite was in. the earth’s 
fall eclipse the direction of the spin axis changed by 
a -negligible but constant' rate . The total angular . ' . 

changes for each of the two parameters over the interval 
were less than 12 arc minutes. The most probable cause 
was due to magnetic, gravitational, and gas leakage torques. 
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I. STAR DATA PROCESSING 

A. Development of the Attitude Determination Method 

The general problem is to find -the orientation as a function of time, 
tj of a coordinate system fixed in a. spinning body as measured in a station- 
ary coordinate system. The basic input measurements which yield this orien- 
tation, or attitude, are the times of transit of known stars across a slit. 

A slit is etched on the focal plane of an optical system as shown in 
Figure 3-1. If the slit is a straight line segment and the optical system 



PHOTOETCHED SLIT 


Figure 3-1: The Slit Plane and Its Normal n 
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is distortionless, then a portion of a plane will be defined which contains 
the slit, and the nodal point of lens system. Given a distant bright point 
source, this source will be sensed by a detector behind the slit if, and only 
if, it lies on the plane defined by the slit and optical system. If a point 
source crosses the plane, the source transits the slit. The instant the 
point source lies in the plane is called the transit time. 

For any transit time of a star, the following equation may be written: 

n . S = 0 (1) 


where 


n = unit vector normal to the plane defined by the optical system 

A 

S = unit vector in the star's direction. 

Since the optical system is fixed within the satellite, n is also fixed 
within the satellite. Its components are therefore most easily written in a 
coordinate system fixed in the satellite. On the other hand, S is most easily 
written in a celestial coordinate system (a system in which the star directions 
are cataloged) . In order that these vectors may be written in the same 
coordinate system, parameters which specify the orientation of the coordinate 
system fixed in the satellite with respect to the celestial system must be 
introduced. However, these are the very parameters which are to be determin- 
ed . 

In order to define parameters which allows n and S to be written in the 
same coordinate system, the physics which govern the motion of the satellite 


(1) At this point it is assumed that the star has been identified although 
such identification is not a trivial problem. Target identification is 
treated in Section I B. 
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must be considered. Two models were applied to the ATS-III satellite atti- 
tude dynamics as follows* 

(1) Spinning, rigid, nearly symmetric, untorqued body 

(2) Body spinning about a fixed axis. 

The first assumption permits a precessional motion of the spacecraft. 
It was discovered during data processing that the angle between the satel- 

c 

lite's angular momentum vector and the principal axis of its largest moment 
of inertia was extremely small (i.e. less than a minute of arc). This ob- 
servation was also made independently through the use of an on-board accel- 
erometer. It was apparent that the motion of the spacecraft principal axis 
was too small to be reliably detected by the star sensor. Consequently, the 
precessional model was abandoned and the simpler fixed axis model was used 

to obtain the majority of the results which follow. A discussion of the 

f 2 1 

analysis of the precessional motion model may be found in reference 1, . 


An analysis of the model assuming a body spinning about a fixed 
axis is included below. A six coordinate frame system is used to . 
define satellite attitude. The development of the coordinate system is 
presented. 

1. Coordinate Systems and Transformations 

Figure 3-2 defines the coordinate .systems S.. , S„ , and S„. The system 

A J. £ J 

S^, is the celestial system, i^ is in the direction of the First Point of 
Aries, k^ is in the direction of the North Pole. The system Sg, is such 
that kg is in the direction of the spin axis of the spacecraft. Hence, for 

A 

ATS-III satellite kg has. approximately the direction of the South Pole. The 
figure, however, is drawn with kg closer to the North Pole for ease of il- 
lustration. 

For the general case, two angles are required to define the direction 
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of the spin axis, k^, with respect to the coordinate system S^. These angles 
are chosen to be called and £ 2 , and are defined as follows: 

* » , c- <■ 

k^-> rotation about 1^ = 1^ 

& 3 rotation g 2 about j 2 j 3 • - 

It is more conventional to locate the direction of an axis with respect 
to S 1 by its right ascension and declination rather than by ^ and g 2 « However, 
if the 'axis is close to the North or South Pole, its' right -ascension is poorly 
defined and undefined if the axis has the direction of the North or South Pole. 
Since the satellite spin axis always lies close to the South Pole, its dir- 
ection is chosen to be defined by an< ^ ?2‘ i s done to eliminate a near 

singularity which would be present in the conventional formulation. However, 

A 

the choice of 5^ and ^ introduces a signularity if the spin axis, k 3> were 
at + i^. 

The system S^ has one axis fixed in the satellite, k^. However, S^ is 
an inertial system if the satellite is spinning about a fixed axis. Another 
system, S^, is fixed in the satellite such that 

A A A 

i. -* i, a rotation s about k = k. . 

3 4 3 4 

The choice of the .position of is not arbitrary but will be defined by the 
position of the slits. . 

Three slits were designed and fabricated to intersect at a point so the 
center slit bisects the angle formed by the outer two slits as shown in Fig- 
ure 3-3. The point of intersection, of course, rotates with the satellite. The 
i direction is chosen such that the (2^, k^) plane contains this common point. 

Let the angle from k^ to the intersection of the slits be a, and let & 5 be 
the unit vector to this point. The nominal position of the center slit is in th 
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FIXED IN SATELLITE 

Figure 3-3: The Orientation of the Slits as- Defined by e, T, cr 
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^ A 

(i , k^) plane. Its actual position may then be defined by a rotation of e 
about k,_. Finally, the two outer slits each make an angle r with the center 
slit. The angles a, e, iTj and ar ® usually treated as constants in- 
dependent of time. The angle s however is time dependent and is assumed to 
be of the form 

s = s + St = % st^ + s t . .. (2) 

0 6 

( 2 ) 

This form is essential if the spacecraft is shaded from the sun , otherwise 
we may usually set s = *s =0. 


Now, the system S^_ is defined so that 

A A A A 

k, -» k c a rotation a about j, = . 

4 5 J 4 5 


System S , j = 1, 


2, 3 are defined by each of the three slits so that 


A 




a rotation 


A 

about, kj. 



here = e + T. 
p 2 = e 

p 3 = e - r 


A A 

Thus, the plane defined by the jth slit is the (i^ , kg .) plane. 


The transformations between the various systems may now be given as; 




(celestial to spin) 


(2) For reasons to be explained in Section I-B. 
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(Spin to satellite fixed) 



1 4 

B (<?, p.)l (satellite fixed to slit) 


where 


H = 



sin § 1 sin § 2 
cos. 

-sin £. _ 

=1 cos g 2 



R = 



sin s 
cos s 
0 


B. 

J 


/cos 

n . COS 

o 

sin 


-cos 

p j 

sin 

a 

f . 
-sin 

i 

0 . cos 

a 

cos 

P i 

sin 

p j 

sin 

G 


sin a 


cos O' 


2. The Constraint Equation 

It is now possible to write Equation (1) in a usable form, and when so 
written, becomes the constraint equation because it offers, a constraint on the 
attitude at the instant of transit of any known star. 


A 

The components of S are simply: 
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S = cos 6 cos a + cos 6 sin a j + sin 6 k^ 


where a 
6 


star's right ascension (known) 
star's declination (known) 


( 3 ) 


while those of 



the unit normal pertinent to the jth slit is simply 



where B„. is the second row of B.. Thus, the constraint equation becomes: 
zj j 

0 = B 2j R H s (4) 

A 

where S is a column vector whose components are given by (3). 

If the following are assumed: 

(1) The satellite is- spinning about a fixed axis so that and £ 2 

are constant. 

♦ 2 3 

(2) s=s 0 +st+%st + 1/6 *s t , where s 0 , s, s, and *s are con- 
stant . 

(3) The optical system is fixed in the satellite and- is distortion 
free. 

(4) The three slits intersect, in a common point, and the angle between 
the center slit bisects the angle formed' by the outer slits; 

then (4) represents the implications of this model and is one equation in the 
9 unknowns ^ 2> s 0 , s, s’, s, e, r, <7. 
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The components of S and the transit time, t, are measured. If a large number 
of measurements are taken, an overdetermined system is defined and a "least 
squares" solution may be sought. 

Obtaining a solution is a non-trivial although standard problem and 
need not be discussed here. Finally, it is noted that s 0 , s, s, and 

s' define the attitude of the satellite, and e, r> and a are parasitic para- 
meters which define the orientation of the slits within the satellite. 

B. Computer Program 

With the absence of both spatial noise and a high incidence of spurious 
transits, the stellar target identification process would be relatively 
simple and straight forward. Almost every star detected would appear as a 
unique triplet of star signals. Time separations of triplet pulses for a 
given star could be compared unambiguously to known separations of candidate 
stars, and a matching of triplets to corresponding stars accomplished . 

Section C of Chapter 2 describes the nature of spatial noise observed 
in the SCNS output data. The high incidence of noise necessitated the use of 
digital filtering techniques which were accomplished with the computer as a 
part of the ground based data processing. 

A program was written to perform the necessary digital filtering and 
to calculate the attitude parameters. The program utilizes signals from 
stars which have been identified by the digital filtering process through 
star transits across known slits. 

A program called SCAT 1 determines the three axis attitude of a 
spinning spacecraft from star transit times obtained from an on-board electro- 
optical star sensor. 

SCAT reads a central crossing time and a measure of the transit width 
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of transit events obtained from the star sensor and written on tape in 
FORTRAN -compatible Binary Coded Decimal (BCD) format. It then identifies as 
many of these transits as possible which originate from- a known star and a 
known instrument slit. The program then uses these identified transits to 
calculate the four most important attitude parameters, other attitude, 
parameters being set equal to their nominal values. 

An earlier program designated BOOTS was designed to accept hand 
identified star transits in a single data set and obtain an attitude solution 
on the basis of the identifications. The program then would "bootstrap" 
through all subsequent sets, identifying stars and computing attitudes. Use 
of this approach was required whenever there were fewer than three solidly 
detected star triplets - the minimum number required for automatic identi- 
fication. 

SCAT requires no such hand identification. SCAT solves for the two 
angles - and f^^which define satellite pointing direction, the initial 
azimuth position s 0 and the spin rate s .. These parameters constitute the 
four unknowns . 

A second fully automatic program, called SCOOT, was designed to accom- 
modate additional unknowns, including the slit parameters e, F, and a and the 
terms s and s', and is essentially an expansion of the program SCAT. SCAT, 


(3) The acronym SCAT was given to the automatic computer program written 
to process star data' which leads to the solution of the four attitude 
parameters § , £ , s 0 and s. The letters used stem from scanning 
attitude and the fact that the solution is obtained quickly, and is 
necessary for purposes of program identification. 
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therefore, is the basic program utilized in the processing of star data. 

The sequence from raw data to final attitude consists of the following 
broad substeps ; 

1. From the raw data, determine^ the scan period' to within 
approximately one part in 10. 

2. Modulo all transit data within a data set consisting of 
approximately 25 scans with respect to this scan period. 

In effect, this operation performs a multi-scan corre- 
lation on the data. 

3. Detect clusters of repeated pulses in the moduloed data. 

This picks out events which are highly correlated with 
the scan azimuth and rejects the transits from spurious 
noise events or from very weakly detectable stars. 

4. Search through the detected pulse clusters for triplets 
of clusters which represent possible transits of a star 
across the three slits. 

5. Compute the separation angles between triplets as these 
angles are projected on the celestial sphere; match these 
separations with the separations on known pairs of stars. 

6. Link the separation star pairs into a closed spherical 
polygon, a star pattern as seen on the celestial sphere. 

If the process is successful to this point, each pulse 
cluster used in a linked triplet may be identified with 
the corresponding star and instrument slit. 

7. Return to the original data and match these identified 
stars to the individual transit times in each scan. 

8. Perform a least squares attitude solution using these 
identified transits. 

9. Using the period calculated in the' attitude solution, 
remodulo the transit data. 

10. Given the above 'attitude solution, generate transit times 
for all stars in the field of view for all three slits. 

11. Retain moduloed peaks that .correspond to generated transit 
times. This picks up peaks of only one or two slits. 

Each peak may be identified with the corresponding star 
and instrument slit. 
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12. Return to original transit data and keep only measured 
transits that can be matched within a specified tolerance 
to. one of the identified, remoduloed peaks. 

13. Recalculate the four attitude parameters using these 
identified, measured transits. 

These substeps are executed by the main program and various- subroutines 
as shown in the flow chart of Figure 3-4. A description of each of the 
routines is included below: 

SCAT 

The program SCAT begins by reading in a star catalog con- 
taining a maximum of 460 stars ordered by visual magnitude. 

For each star the serial number, the general catalog number, 
right ascension, declination, and magnitude are read in'. The 
direction cosines of each star are calculated’ and put in the 
array SDC. The limits on -the field of view are calculated from 
the estimated pointing direction (RA, DE) and only stars brighter 
than a specified magnitude (XLM) and within the field of view are 
considered throughout the remainder of the program. The index 
numbers of ‘ these NLS stars in the star catalog are put in the 
list IAS. The angular separations between these stars are then 
calculated and put in the array SEP. 

A header card for a data set is now read in and the input 
tape is searched for this data set. The centers and widths of 
the NTT transits are put in the arrays TIQ and IWD, respective- 
ly. A maximum -of 4,000 transits is read in. 

Pulses from the. sun sensor have been written on the tape as' 
having a width of -100. These are now detected and eliminated. 

The subrouting PERIOD is now called' to determine an esti- 
mate of the scan period given estimates of the upper and lower 
limits of the period and the measured transit times. 
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SCAT 



Print list of detected 
trans its 


Print list of triplets 



Print list of all 
identified triplets. 



MODULO 

Modulo all transit data with 
respect to this scan, period. 
Histogram the data and choose 
the bins that have more pulses 
than a specified threshold as 
good transits. 


SLMAT 

Search through detected peaks for 
triplets of peaks which represent 
possible transits of a star across 
three slits. 


IDONT 


Cycle through triplets, 3 at a time 
in order of decreasing total count, 
calculating separations. 

NO 


Successful? 


J?ES 

Try to match separations of all 
remaining triplets to brightest 
identified triplet with matrix of 
known separations of stars to corre- 
sponding brightest identified star. 
Each triplet thus matched may be 
identified with the corresponding 
star. 


LINK 

Attempt to link 
the three separa- 
tions into a closed 
spherical polygon 
which is identical 
to a closed 
spherical polygon 
composed of the 
separations of 
three known stars. 
Each of the three 
triplets may then 
be identified with 
the corresponding 
star. 


Figure 3-4: Computer Program, SCAT -Flow Chart 
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All transit data are moduloed with respect to this scan 
period by subrouting MODULO. The centers of the ITT detected 
peaks are returned in the array AVTT and the number of transits 
in each peak in the array NAV, 

The subroutine SIMAT is now called to search through the 
detected peaks and attempt to identify as many as possible as 
members of triplets of peaks having equal spacing. The centers 
of the three peaks belonging to each of the NOT identified 
triplets are returned in the array GT. 

If less than three triplets have been detected, the 
program quits this data set and reads in a new data set. 

If at least three triplets have been detected, the sub- 
routine IDONT identifies as many as possible of these triplets 
with known stars from the catalog . The triplets that have been 
identified are tagged by setting ID(n) equal to the index of 
the corresponding catalog star. 

The identified triplets are printed with the corresponding 
star and the unidentified triplets are printed with the comment; 
"Unidentified triplet". The identified peaks are put in the 
array TTI and the corresponding star and instrument slit in 
the arrays MST and MSL, respectively. If no triplets have been 
identified, the program quits this data set and reads in a new 
data set. 

An estimate of the origin, SO, of the azimuth angle is made 

P 

using the first identified star crossing the counter slit. The 
subroutine COMPAR returns to the original transit data, matches 
the measured transit times to the identified peaks, and retains 
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only those transit times that are within a specified tolerance 
of an identified peak. These matched transit times are put in' 
the artay TTI and are used by the subroutine ATTUD4 to solve the 
four unknown problem and calculate the four attitude -parameters . 

The subroutine MODULO is again called to remodulo the 
transit data using the period determined by the attitude 
solution. The ITT detected peaks are put in the array AVTT. 

Using the calculated attitude parameters, the subroutine 
GENTIME generates a transit time for each star in the field of 
view going past all three slits. These generated transit times 
are put -in the array TTI and are compared with the remoduloed 
peaks by the subroutine COMPAR. Only those peaks that are 
within a specified tolerance of the generated transit times 
are retained. They are stored in- the array TTI. This pro- 
cedure picks up peaks that are not in triplets but originate 
as the transit of only one or two slits by a star. The corre- 
sponding known star and slit are stored in the arrays MST and 
MSL, respectively. 

The subroutine SIMAT again detects as many triplets as 
possible and these triplets are printed out. 

The origin of the azimuth angle is again estimated from 
the first identified star. 

The subroutine COMPAR returns to the original transit times 

and matches the measured transit times to the identified, remodu- 
' « 

loed peaks. Only those transits that are within a specified 
tolerance of the identified peaks are kept; These are returned 
in the array TTI and the corresponding stars and slits in the 
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arrays MSTT and MSLT. These matched transit times are then used 
by the subroutine ATTUD4 to recalculate the four attitude parame- 
ters . 

PERIOD 

The routine PERIOD uses the r aw transit data and estimates 
of the lower and upper bounds (TMIN, TMAX) of the spin period to 
determine a reasonable approximation of the period. 

The routine cycles through all of the transit times, forming 
the difference between the transit time and all succeeding transit 
times. The differences which lie between TMIN and TMAX are stored 
in the array D-T . These differences are sorted into order by the 
subroutine SORT. The program then cycles through the differences , 
counting how many fall into each bin of a histogram whose bin 
width is specified (TP), and lower and upper bounds are TMIN - TP 
and TMAX + TP. 

The mean of the differences in the bin containing the 
largest number of differences is chosen as the period and the 
variance is calculated. 

MODULO 

This routine moduloes all transit times to within one scan 
period. 

The transit data (TT) is cycled through, each transit being 
moduloed with respect to the period (PER). The moduloed transits 
are then sorted into order and put into a histogram whose bin 
width is twice a specified quantity, SW. The lower edges of the 
overlapping bins .begin at intervals of SW and the range of the 
histogram is from 0 to PER. 
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The bins that contain more transits than a given threshold 
are considered to be solidly detectable peaks. The threshold is 
equal to a specified percentage (THRESH) times the number of 
effective scans. If more than one acceptable bin are adjacent, 
only the bin with the largest number of transits is chosen. The 
mean positions of the transit times in the .detected bins are 
.stored in the array AVTT and the number of transits in the corre- 
sponding bin is stored in NAV. 

SIMAT 

The routine SIMAT searches through the detected peaks for 
triplets of peaks separated by equal spacing. The lower and 
upper limits of the allowed spacing (HDTMIN and MDTMAX) are 
equal to the minimum and maximum possible differences of transit 
times of one star across two slits. Each detected triplet 
represents a possible transit of a star across the three slits. 

The routine overlaps the first few peaks onto the end to 
make a continuous band in order to detect triplets which overlap 
the zero mark. The routine then cycles through these peaks (WT) , 
attempting with each peak to find two other peaks which are 
equally spaced (i.e., equal within the tolerance DT) . These 
spacings must be between HDTMIN and HDTMAX, The positions of the 
three peaks of the detected triplets are stored in the array 
GT(N, 3). The number of transits in the detected peaks are 
stored in NE(N, I, 1), the corresponding identified star in 
NE(N, I, 2) and the corresponding slit in NE(N, I, 3). 

IDONT 

The routine IDONT identifies as many triplets as possible 
with known stars . 

The routine assumes that each triplet is the result of 
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transits of a target star across all three slits. Using the. slit 
parameters (SIG, GAM, EPS), the spin rate (W) , and the transit 
times across the first and third slits (TT(N, 1), TT(N, 3) for 
each triplet, the direction cosines of the line from the center 
of the celestial sphere to the target star with respect to the 
spin axis is calculated for each triplet. These direction 
cosines are stored in the array U(N, 3). 

All of the triplets are now cycled through, three at a 
time, each time taking the three triplets having the largest 
total count of transits and not yet used together. For each 
set the separations (TRIP) of the projections on the celestial 
sphere of each of the three triplets from each other are cal- 
culated from the previously calculated direction cosines of each 
triplet. The subroutine LINK is now called to attempt to link 
the separation star pairs into a closed spherical polygon which 
is identical to a closed spherical polygon composed of the 
separations of three known stars. If this linkage is successful, 
the triplets may be identified with the corresponding stars . 

If no such linkage is possible for any set of three triplets, 
the routine returns to the main program. 

If -a linkage has been established, the remaining triplets 
are cycled through. For each triplet the separations (Z) between 
the projection of this triplet and the projections of the first 
three identified triplets on the celestial sphere are calculated. 
These are compared to the matrix (SEP) of separations of the known 
stars with the first three identified stars. If, for a given 

i 

known star, these separations match within a specified tolerance 
(ES), the triplet is identified with the corresponding known star. 

The indices of the IKK identified triplets are stored in the 
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array ND and the indices of the corresponding known stars are 
stored in the array ID. 

LINK 

The routine LINK attempts to link the three separations of 
the projections on the celestial sphere of three triplets into 
a closed spherical polygon x^hich is identical to a closed spher- 
ical polygon composed of the separations of three known stars. 

If this procedure is successful, each of the three triplets may 
be identified with the corresponding known star. 

The matrix (SEP) of the separations of the known stars is 
searched for all star pairs whose separation is similar to any 
of the separations (TRIP) of the three triplet pairs. The 
indices of all star pairs whose separations, are within a speci- 
fied tolerance (El) of one of the separations of the three trip- 
let pairs are placed in a table ITT, each of whose three columns 
is headed by one of the triplet pairs. The difference betxtfeen 
the separations of the known star and triplet pairs is placed in 
the table ETT. NIS(J) is the number of entries placed in each 
column of the table. 

The entries of the table are then cycled through and a 
continuous chain of indices from column one to column two, from 
column two to column three, and from column three to column one 
is searched for. If only one such linkage is found,® the knox^n 
star indices can be identified x<rith the corresponding triplets 
heading the table. If more than one linkage is found, the link- 
age having the smallest root mean square error is accepted. The 

* 

indices of the identified stars are stored in the array ID, 

COMPAR 

The routine COMPAR compares an ordered list (TTC) of transits 

* 
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or peaks having NC elements with an ordered list (TIM) of 

transits or peaks having NM elements. Only the elements of 

% 

TTC which are within a specified tolerance (TOL) of any 
element of TTM are retained. The NCP retained elements are 
stored in the array TTM. The subscript's of the star index 
number (MST) and slit number (MSL) corresponding to a retained 
element of TTC are adjusted to correspond to the new position 
of the element in the list TTM. 

ATTUD4 

The routine ATTUD4 uses the identified, measured transits 

and initial guesses of the attitude parameters to calculate the 

attitude parameters. The unknown attitude parameters are: XII 

XI2 which are the pointing angles, S° which is the origin of 

the 1 , axis at time t = 0, and SD which is the angular spin 
4 

rate. The routine solves the basic constraint equation by an 

iterative scheme. In each iteration the best solution is found 

by a least square fit of all of the identified transits. The 

procedure is described in detail in the Feasibility Study Report 

[3] 

for Self-Contained Navigation System, pp. 255-257 

The routine first calls HMAT to calculate the matrix H 
which transforms coordinates from the astronomic coordinate 
system S^ to the inertial system S^ whose k^ unit vector is 
parallel to the angular momentum vector. 


For each transit the following steps are carried out. 

The Euler matrix R is calculated by the subroutine RMAT. 

The matrix R transforms coordinates from the S^ system to 
the coordinate system which is fixed to the spinning body 
and whose k^ unit vector lies parallel to the principal axis 
of the body. The B matrix which transforms coordinates from 
the system to the coordinate system whose unit 
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vector lies perpendicular to the jth slit plane has already been 
calculated for each slit. The products of the second row of B 

t 

with the matrix R and with the derivative of R are' calculated 

and stored in B2R and B2RS. The products of the ‘direction 

\ 

cosines (SDC) of the star corresponding to the transit under • 
consideration with the H matrix and its partial derivatives are 
calculated and stored in HS , HIS, and H2S. These quantities 
are then used to calculate the various coefficients which are 
normalized and converted into time units before being, stored in . 
the matrix V. 


The matrix 


[(V^V) 


-WIT ... 


J is calculated by the subroutine 
GENINV and the corrections (DX) to the attitude parameters .are 
calculated. After the attitude parameters are updated, the 
correction terms are compared with the error .tolerances (EPS). 
If any of the correction terms are larger than its correspondin 
error tolerance, the new attitude parameters are used as input 
and the calculation is* repeated. Up to 50 such iterations are 
allowed in searching for a converging solution. 


O' 

o 


After a converging solution has been achieved, the transits 
with residuals (ATK) greater than twice the RMS residual (SIG) 
are eliminated and the entire calculation is repeated; 


Finally, the standard deviations of the four unknowns are 
calculated. 

GENTIME 

Given the spin axis direction (included in the H matrix), 
the initial azimuth (SO) and the azimuth rate (SD) from the 
attitude solution, GENTIME generates expected transit times for 
all of the stars in the field of view (LS) going past all three 
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slits . 

The direction cosines of each star are transformed from the 
coordinate system to the coordinate system. These are 
stored in the array HS. For each star, the azimuth (S) corre- 
sponding to the transit time for each slit is calculated and 
the transit t'ime corresponding to the azimuth is stored in 
TTC. The corresponding star index and slit number are stored 
in the arrays MTX and MLX, respectively. 

These generated transits are sorted in order of occur- 
ence. All transit pairs -which are closer together than the 
rotational slit width (ROTSW) may lead to ambiguous transits 
since the transits would overlap. The generated transits are 
searched for such pairs. They ar.e then printed out as having 
possible transit time conflicts. 

GENINV 

The routine GENINV takes the matrix A, which has NR rows 
and NC columns and calculates the matrix A', the generalized 
inverse of A. The transpose of A' is stored in A and is 
returned to the calling program. 


• B2VEC 

The B matrix transforms coordinates from the system, in 

A 

which the k^ unit vector is parallel to the spin axis, to the Sg_. 
system in which the j ^ unit vector is normal to the slit plane. 
For each slit there’ will be a separate B matrix. Since the com- 
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always equal to zero, the easiest set of equations to solve will 
involve only the second row of the B matrix. These equations 
will be 


0 = B 2 S 6 ’ 

l 

where is the second row of B. For each slit the elements of 
the second row are 

- sin p cos cr, cos p, sin p sin c, 

where p and cr are the slit parameters. 

The routine B2VEC calculates the second row of the B matrix 
for each slit and stores, the results in the array B2(3, 3). 


HMAT 

Given the angles XII and.XI2, the routine HMAT calculates 
the elements of the matrix H and the partial derivatives .(HP(1, N, N), 
HP(2, N, N)) of these elements. The matrix H transforms coordinates 

A 

from the astronomic coordinate system S. to the inertial system 

A A 

whose unit vector is parallel to the angular momentum 
vector. The elements 'of the matrix are 


H =4 


V 


cos § 2 sin sin § 2 
cos §1 

-sin cos 5 2 


0 

sin § 2 


-cos ^ sin § 2 
sin § 1 

cos § 1 cos 5 2 


RMAT 

Given the angle. SPF, the routine RMAT calculates the elements 
of the Euler matrix R. and the partial derivatives (RPO) df these 
elements. The matrix R transforms coordinates from the inertial 
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A A 

system S~ , whose k„ unit vector is parallel to the angular 
momentum vector, to the coordinate system S^, which is fixed in 
the spinning body and whose k'^ unit vector lies parallel to the 
principal axis of the body. The elements of R are 


R = 


cos g sin p 0 
sin j3 cos p 0 
0 0 1 


where p is the azimuth angle SPF. 


SORT 

The routine SORT (or MORT) sorts the N elements of the list X 
from the smallest to largest (algebraically) and returns the list 
in the array X to the calling routine. The list KEY, which is 
also returned,, contains the position in the original list of each 
element in the new list. 

A dictionary of selected variables used in SCAT is included as Table I. 
A listing of parameters which must be specified at the beginning of the 
program along with typical values for each (shown in parentheses) is included 
as Table II. 


The first input to the program is a star catalog containing a maximum 

of 460 stars — one card for each star. The cards are read in with the for- 

mat (13, 17, 2F20.13, F10.2). Each card contains the following information 
for one star: serial number of the star, general catalog number, right 

ascension in radians, declination in radians and magnitude. 

For each data set on the input tape for which an attitude determination 
is desired, a header card is read in. The header card contains the data set 

number and a comment. The header cards are read in with the format (15, IX, 

13, 31X, 5A8). The last such card should be a blank card. 
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SDC (N , 3) 
SEP (N , N) 
IAS (N) 

NLS 
TTQ (N) 
IWD(N) 
NTT 

AVTT(N) 

NAV(N) 

ITT 

GT(N, 3) 

NGT 

ID(N) 

TTI(N) 

MST(N) 

NSL(N) 

NLP 

PER 

SD 

SO 

XII 

XI2 


TABLE 3-1 

DICTIONARY OF SELECTED VARIABLES IN SCAT 


direction cosines of catalog stars 
separation matrix of catalog stars 

list of indices of stars in catalog that are in field of 
view 

number of stars in catalog that are in field of view 
center of measured transit 
width of measured transit 
number of measured transits 

til 

average of transit times in N good moduloed peak 

til 

number of transits in N good peak 
number of good peaks 

til 

positions of peaks of N triplet 
number of triplets 

til 

index of corresponding catalog star for N identified 
triplet 

til 

position of N identified peak or transit 

star corresponding to N^ identified peak or transit 
> til 

slit corresponding to N identified peak or transit 
number of matched peaks or transits 
period 
spin rate 

origin of axis at time t = 0 
pointing angle of angular momentum vector 
pointing angle^ of angular momentum vector 
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NS 

XIM 

MAXND 

T0L2 

T0L3 

GAP 

TGAP 

PERSLT 

HDTMIN 

HDTMAX 

THRESH 

EID 

MIN 

MAX 

SW 

XLD 

XHD 

XI1D 

XI 2D 

RAD 

DED 

SIGMAD 

CGAMD 

EPSD 

ROTSW 
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TABLE 3 -II 

LISTING OF SCAT PARAMETERS 
TO BE SPECIFIED 

number of last catalog star to be read in (460) 

magnitude of faintest catalog star desired in list IAS (4) 

maximum allowed width for acceptance of transit (1000) 

tolerance for comparing identified with measured transit 
times (0.000110) 

final tolerance in slit matching to obtain triplets (0.000220) 

if difference between two consecutive transit times are 
greater than this quantity, a writing gap is assumed to 
have occurred (1.200) 

actual length of writing gap (1.260) 

bin size for determination of scan period (0.000050) 

minimum peak-to-peak separation for triplet search (0.012000) 

maximum peak-to-peak separation for triplet search (0.012000) 

threshold of relative count frequency for detection of 
repeatable triplets (0.30) 

tolerance in degrees for comparing the separation of 
measured triplets with star separations (0.20) 

minimum limit on the scan period (0.564800) 

maximum limit on the scan period (0.565800) 

bin width for selection of repeatable transits (0.000110) 

limits (in degrees) of field of view (50.0 and 80.0) 

initial angles of rotation involved in transformation from 
system to system (180.0 and 0.0) 

right ascension (in degrees) of initial pointing direction 
of angular momentum vector (135.0) 

declination (in degrees) of initial pointing direction of 
angular momentum vector (-89.5) 

slit parameters in degrees (84.43651, -11.08676, 0.91898) 

rotational slit width in arc minutes (7.6) 
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The data is written on the input tape in data sets consisting of all 
the transits in approximately 25 scans. For each, data set, the first record 
contains the data set number, the number of transits in the set, and the 
starting time of the set in hours, minutes, and seconds in the format. (6X, 
13, 16, 215, F12.6), The next records contain the list of transit times 
and transit widths — eight transits per record. The records are written 
on the tape with the format (8(F10.6, 15). The widths of the transits 
corresponding to sun pulses have been written on the tape as -100. 
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II. ATTITUDE DETERMINATION 

The numerical results of satellite attitudes are documented within this 
section. The results were obtained from data acquired at various times be- 
tween 13 March and 4 October 1968. 

Figure 3-5 illustrates the geometric relationship between the SONS field 
of view and sources of bright object interference. The effectiveness of the 
instrument field of view was influenced by bright object constraints when, 
under certain conditions, the number of solidly detected targets was reduced 
thereby affecting the accuracy of the attitude solution. The field of view 
was further affected by the loss of sensitivity for off-axis sources -- a 
condition brought about by the nature of the optical system. This was dis- 
cussed in Section. I of Chapter 1. 

■ A. Precessional Motion Model 

. On Day 73 (13 March 1968) the satellite was in the earth's eclipse for 
approximately one hour. During that interval star sensor measurements were 
acquired which utilized a mathematical model based upon the following assump- 
tions ; 

(1) The satellite is torque free, rigid, and nearly symmetric. 

(2) The angular momentum of the satellite is nearly parallel to 
the principal axis of largest moment of inertia. 

For reasons to be explained later, this model was applied to Day 73 data only. 
The model, termed the Precessional Motion Model, considers precession of the 
satellite principal axis about the angular momentum vector . 

The results from measurements taken over an interval of approximately 
13 seconds are shown in Table 3-TII. All angles listed in the table are 
measured in degrees and all rates are in degrees/seconds. Figure 3-6 is 
included as an aid in understanding the notations used in Table 3 -III. In 
each of the distinct cases shown in Table 3-III no star dimmer than magnitude 
2.2 was utilized. If all measurements were perfect, these cases would yield 
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Figure 3-5: Geometric Relationship of the SONS Field of View 




TABLE 3 -III 

SOLUTION FOR PRECESSION MODEL (DAY 73 DATA, 10 UNKNOWNS) 


oo 

vo 


Initial 
Spin Angle 


Precession 

Rate 


Half Cone 
(Nutation) 
Angle 


Initial 

Precession 


Spin Rate 


Angles Which 
Define Direction of 
Angular Momentum 


Right Ascension and 
Declination of 
Angular Momentum Vector 


t 

.A- 




Angles Which 
Define Location 
of Slits 


_A_ 


~\r 


Case 

Ho. 

e 




1 

h 

S 2 

a 

r 

e 

O' 

6 

1 

.040 

- 15.656 

- 11.222 

252.672 

587.328 

180.224 

.094 

84.264 

- 11.155 

.185 

67,36 

- 89.757 

2 

.008 

- 68.720 

- 21.935 

- 54.255 

598.042 

180.197 

.093 

84.272 

- 11.057 

.192 

64.81 

- 89.782 

3 ** 

.005 

106.187 

- 22.589 

117.065 

598.595 

180.202 

.105 

84.373 

- 11.122 

.186 

62.42 

- 89.772 

4 * 

.003 

51.419 

- 18.211 

185.589 

594.318 

180.213 

.093 

84.294 

- 11.170 

.186 

66.53 

- 89,768 

5 ** 

.003 

104.132 

- 22.445 

119.117 

594.451 

180.197 

.119 

84.391 

- 11.114 

.184 

60.43 

- 89.773 

6 

.017 

- 186.766 

- 13,423 

48.939 

589.435 

180.200 

.111 

84.401 

- 11.105 

.184 

61.07 

- 89.771 


*, ** Indicate cases which had common initial times 
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AAA A 

i 3 5 ] 3 »k 3 - INERTIA SYSTEM WITH k 3 IN DIRECTION OF ANGULAR 

MOMENTUM OF SATELLITE 

AAA A 

'6 > k 6 “ B0DY FIXED SYSTEM WITH kg IN DIRECTION OF 

PRINCIPAL AXIS OF GREATEST MOMENT. 


Figure 3-6: Relationship of Body-Fixed Reference 

System to Angular Momentum Frame 
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the same values for 8 .of the unknowns. The values of i|f 0 and § 0 , however, 
would vary because not all cases shared a common initial time. The initial 
time was the same for Cases 2, 4 and 3,5, Thus, except for the effect of 
errors. Cases 2 and 4 would process a common \jf 0 and a common 8 This is also 
true for Cases 3 and 5. In Table 3-IV the six cases are briefly described. 


TABLE 3-IV 

DESCRIPTION OF CASES LISTED IN TABLE 3 -III 



Number of Stellar 
Targets Used^^ 

Time Interval (sec) 


Case No. 

Begin 

■ 1 ■ 

Comments 

1 

3 

.888455 

13.989168 

Cases 1 and 2 utilize the 

2 

3 

.263583 

13.364247 

same stars, but no common 
transit times. 

3 

6 

.265618 

13.421220 

Uses all point sources 
reliabily detected. 

4 

4 

.263583 

13.989168 

Same as Case 2 except 
Mercury used. 

5 

5 

.265618 

13.417741 

Same as Case 3 except 
Mercury not used. 

6 

6 

.267434 

14.044438 

Cases 3 and 6 utilize same 
stars, but no common 
transits 


Table 3-III indicates an extremely small cone angle. Also the angular 
momentum vector of the satellite points very close to the celestial South 
Pole. In addition, the variation of the angles 0, Q > F, and e is 

about one minute of arc. 

For small 0, the angles § and if are poorly defined, and only their sum 
is important. The variations of in Cases 2 and 4 implies that the measure- 
ment of 9 is in the noise, (i.e. -the error in the determination of 9 is about 
the size of 9 itself). 


(4) Venus and Mercury are included as potential targets. 
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From i|f, §, and 9 it is possible to determine the ratio of the largest 
to the smallest moment of inertia 


A 

C 


# cos 6 
+ I cos 9 


(5) 


For each of the six cases, 1.023 < A/C < 1.040, 

Since the angle 9 is extremely small, the precession model was abandoned 
in favor of a model which assumes 0 = 0 (i.e. the satellite is assumed to 
rotate about a fixed axis). However, if 0 were about two minutes of arc or 
larger, the instrument accuracy is sufficient to justify the use of the pre- 
cession model. 

This assumption of a fixed axis of spin is valid because the spacecraft 
contains a precession and nutation damper which apparently restricts the 
angle between the satellite principal axis and angular momentum vector to less 
than one minute of arc . 

B. Fixed Spin Axis Model 

The assumption of a fixed spin axis formulates the mathematical model 
in terms of eight unknowns. The eight unknown problem includes a spin accel- 
eration term, s, introduced by the fact that the satellite, in passing from 
full sunlight into earth's shadow cools. This cooling causes physical 
shrinkage of the satellite with the consequence that the spin velocity 
increases to conserve momentum. Therefore the velocity term is not a constant. 

The following terms define the eight unknowns of the model; 

s g . - initial position of £ body axis 

» 

s - spin rate 

«• 

s - spin acceleration 
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} angles which define direction of spin axis 
with respect to celestial sphere 

} angles which define location of slits with 
respect to body axis 

Attitude solutions obtained from data taken on Days 73, 95, 96 and 97 
is presented in Table 3-V which lists the derived values of the eight unknowns 
s xs in degrees/second, s is in degrees/second , and all other angles are 
expressed in degrees. O' and 6 the right ascension and declination, respect- 
ively, of the spin axis are listed as are the number of transits used; E, the 

_ ^ 

RMS value of the residuals (in seconds x 10 ); and the duration of time over 

which the transits were gathered (in seconds). The three solutions obtained 
for Day 97 result from three distinct sets of transit times with a common 
interval of time. A spacecraft maneuver was performed on 22 March, thereby 
explaining the relatively large change in attitude observed between Days 73 
and 95. Also shown in Table 3-V are estimates of the standard deviation of the 
eight unknowns which were obtained from the residuals of the final solution. 

maximum value of from the three values of E; obtained 
from Day 97 
minimum value 

Using similar definitions for other variables, we obtain 

21,6 sec = 6 
14.4 sec 

2 

,03 sec/sec 
95.4 sec 
36 sec 
1,584 sec 
11,736 sec 



Let ?l = 

h - 




a 


r 

e 
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TABLE 3-V 

SUMMARY OF ATTITUDE RESULTS FROM DATA TAKEN ON DAYS 73, 95, 96 AND 97 


Day 

Si 

s 2 

s o 

s 

s 

r 

e 

a 

O' 

6 

No. of 
Transit 

2 

Duration 

73 

180.221 

.121 

256.694 

575.9005 

.00012 

-11.1688 

.193 

84.280 

61.277 

-89.748 

444 

31.8 

104 

95 

180.443 


126.143 

576.1102 

.00018 

-11.0943 

.196 

84,410 

78.508 

-89.548 

1101 

27.9 

73 

96 

180.441 

.089 

47.497 

576.3982 

.00010 

-11.0834 

.220 

84.436 

78.548 

-89.550 

867 

34.9 

84 

97 

180.457 


126,640 

576.2715 

.00013 

-11.1117 

.227 

84.368 

78.899 

-89.534 

565 

29.0 

100 

97 

180.463 

.086 

197.649 

576.2735 

.00014 

-11.1241 

.219 

84.348 

79.424 

-89.529 

611 

29.2 

74 

!Lj 

180.457 

.087 

-233.360 

576.2715 

.00013 

-11.0976 

.229 

84.392 

79.225 

-89.535 

744 

28.8 

99 


Day 

*«!> 


<7(0 

a(s) 

a(s') 

<j(r) 

<7(s) 

a(a) 

73 

13.08 

8.63 

9.44 

.446 

.0078 

32,26 

29.33 

54.17 

95 

5.90 

4.67 

8.76 

.328 

.0086 

22,85 

18.70 

44.63 

96 

10.78 

5.68 

11.92 

.391 

.0089 

35.15 

26.65 

67.26 

97 

8.62 

6.38 

12.00 

.274 

K&9 

34.12 

27.88 

69.45 

97 

8.31 

6.17 

11.93 

.400 

.0105 

32.36 

26.48 

66.09 

97 

7.52 

5.53 

10,53 

.257 

.0051 

28.94 

23.79 

58.92 


2 -- measured In microseconds 

Duration --in seconds 

Standard deviations -- in seconds of arc 
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Note that all parameters are quite accurately determined except for cr 
(angle which determines the position of the slit concurrence) and a (right 
ascension of the spin direction). The right ascension. O', is poorly deter- 
mined because the spin axis lies close to the South Pole, and a is poorly 
determined because the angle between the center slit and the outer slit is 
a relatively small angle (11?1020).- 

In Figure 3-7, the computed direction of the satellite spin axis is 
plotted as a function of time for Day 201 (19 July), (the day on which the 
second probe was ejected). The satellite on that day was in full sunlight 
over its entire orbit. The results shown in Figure 3-7 were derived using 
much poorer stellar geometry than that used to derive the results of Days 
73, 95, 96 and 97. Poor geometry has a dual cause: (1) Direct blocking 

due to the sun's interference and (2) blocking due to sun's reflections off 
the launch mechanism arm used to eject the first and second probes. Only 
three stars were regularly detected. This fact led to solutions having a 
standard deviation of about 1.5 minutes of arc for the direction. There- 
fore, it is probably not permissible to interpret the variations of the 
spin direction given in Figure 3-7 as real. These variations are almost 
entirely due to errors in the transit times. 

In obtaining the results shown it is assumed that the satellite spins 
uniformly about a fixed axis., and that the parameters which characterize the 
slits are fixed so that T = -11.087°, e = .199°, a = 84.437°. 

Figure 3-8 is similar to Figure 3-7 except that the time interval is 
extended over 28 days. Note, however, that a and 9 are much smoother in Fig- 
ure 3-8 due to the improved stellar geometry on Days 249 to 277. On Day 201 
only three stars were regularly detected. On Days 249 - 277 when the satel- 
lite was in the earth's shadow, eight stars were regularly detected. 

The results of Figure 3-8 indicate a systematic increase with time for 
both the right ascension and declination of the spin axis. This increase 
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is small but real. Note that the direction of the spin axis varies by less 
than 14 minutes of arc from its centroid in the 28 day interval. These small 
changes are the effect of torques acting on the satellite, due most probably 
to magnetic and gravitational effects. 

Table 3-VI provides a detailed tabulation of the computed spacecraft 
attitudes from Day 249 to 277. 

C . Slit Parameters 

The slit parameters or, e, and T defined in Figure .3-3 represent 
parasitic unknowns in the attitude problem. A shift in the spin axis relative 
to the star sensor will change either e, cr, or both; but such a shift will 
not affect T. This, is because e and cr define the orientation of the sensor 
with respect to the satellite spin axis, but T does not. The slit separation, 
T, depends only on the mask on which the slits were etched and the location 
of this mask with respect to the instrument's optical axis. Unless the SCNS 
instrument undergoes significant internal displacements. Twill be independent 
of time. However, e and cr may very well be time dependent for the spacecraft 
is not a rigid body. 

A means of measuring T independently (from that of treating it as an 
unknown in the attitude problem) exits. T was therefore obtained by measuring 
the angle between the slits bn the mask and the location of the mask with 
respect to the optical system of the instrument. The r.esultant value for F 
was -11.102 + .02 which is in good agreement with values obtained . from the 

stellar transits. 

It should be recognized that e and o were not measured by an independent 
means, due to the difficulty of accurately measuring the orientation of the 
SCNS instrument' with respect to the satellite spin axis. A further discussion 
relative to slit characterization is included as Appendix D. 

All results after Day 201, except for Day 263, were obtained by assuming 
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TABLE 3 -VI 

ATTITUDE AND SLIT PARAMETERS OBTAINED FROM DAY 249 to 277 


DAY 

h 

m 

s 



S ° 

s 

3 

s x 10 

8 

S X 10 

€ 

<r 

249 

3 

i 

52. 

179.55261 

-.78828 

-258.37690 

636.76887 

.201162 

-7.26640 

.35356 

84.21268 


3 

5 

46. 

179.54909 

-.78427 

-289.09922 

636.81403 

.183327 

-5.42379 

,37506 

84.20182 


3 

9 

39. 

179.54789 

-.78369 

-226.55999 

636.85543 

.167515 

-2.99884 

.39733 

84.20899 

2S0 

2 

56 

52. 

179.53187 

-.78772 

-152.86261 

636.72988 

.210150 

-6.44431 

.40470 

84.19895 

269 

2 

55 

26. 

179.34523 

-.80512 

-96.22028 

636.94307 

.136820 

.223200 

.38259 

84.20411 


2 

59 

59. 

179.34562 

-.80565 

-85.69124 

636.97998 

.128174 

-.855446 

.38844 

84.20344 


3 

13 

57. 

179.34508 

-.80512 

-133.11407 

637.08027 

.109005 

-'.828800 

.36724 

84.20235 


3 

18 

9. 

179.34644 

-.80588 

-145.44161 

637.10730 

.107780 

-3.57838 

.37354 

84.20657 


3 

22 

14. 

179.34659 

-.80503 

-290.99994 

637.13277 

.101684 

-1.22416 

.37648 

84.20403 

270 

2 

46 

17. 

179.33416 

-.80336 

-350.14623 

636.86344 

.153739 

11.7902 

.39458 

84 . 18103 


2 

56 

1 . 

179.33638 

-.80491 

-275.72152 

636.95077 

.142589 

-8.02144 

.38864 

84.18266. 


3 

6 

0 . 

179.33779 

-.80517 

-318.43541 

637.02830 

.118495 

-1.38603 

.39836 

84.17786 


3 

16 

0 . 

179.33597 

-.80331 

-40.52580 

637.09696 

,104282 

.244014 

,39144 

84.17348 


3 

21 

0 . 

179.33620 

-.80239 

-66.69590 

637.12852 

.102382 

-1.17487 

.38537 

84.17294 

271 

2 

39 

2. 

179.32657 

-.80536 

-322.55184 

636.77110 

.193076 

-10.3370 

.41091 

84.18262 


2 

53 

10. 

179.32573 

-.80398 

-278.3781.4 

636.91178 

.143635 

-2.36739 

,41751 

84.18034 


3 

7 

0 . 

179.32328 

-.80266 

-74.80117 

637.02107 

.122294 

-4.03589 

.41989 

84.18037 


3 

18 

0 . 

179.32422 

-.80366 

-95.73663 

637.09574 

.106830 

-1.75402 

.40270 

84.18009 

272 

2 

40 

2. 

179.31560 

-.80302 

-274.08844 

636.79835 

.184237 

-6.80615 

.42820 

84.17452 


3 

10 

38. 

179.31256 

-.79986 

-135.73530 

637.05952 

.113796 

-2.20786 

.42118 

84.17376 


3 

18 

2. 

179.31566 

-.80294 

-230.36006 

637.10804 

.106853 

-2.39279 

.44390 

84.17316 

273 

2 

31 

2. 

179.31157 

-.81410 

-356.79356 

636.67081 

.233408 

-11.0348 

.47958 

84.18827 


2 

45 

59. 

179.30732 

-.80427' 

-140.35696 

636.84500 

.166588 

-8.41468 

.44847 

84.17956 


3 

6 

0 . 

179.30527 

-.80423 

-183.38722 

637.01312 

.121298 

-2.98056 

.44243 

84.17376 


3 

16 

0 . 

179.30711 

-.80546 

-274.54730 

637.08176 

.108437 

-1.87672 

.43239 

84.17511 

274 

3 

14 

47. 

179.29267 

-.79744 

-134.57631 

637.07120 

.104404 

4 . 70085 

.43263 

84.16640 

276 

2 

37 

47, 

179.27728 

-.80025 

-8.33732 

636.72993 

.200037 

-8.15586 

.38674 

84.18593 


2 

45 

0 . 

179.27690 

-.79926 

•-46.58335 

636.80956 

.170455 

-5.70767 

.35590 

84.17783 


2 

56 

0 . 

179.27918 

-.80086 

-197.61615 

636.91136 

.140322 

-2.70504 

.35141 

84.18361 


3 

7 

0 . 

179.27557 

-.79927 

-286.99748 

636.99751 

' .119918 

-1,67184 

.36614 

84.17696 


3 

14 

45. 

179.27959 

-.79721 

-350.48496 

637.05142 

.105532 

3.86296 

.35768 

84.17520 

277 

2 

55 

2. 

179.27051 

-.79829 

-23.76617 

636.88593 

.145909 

-3.44978 

.41227 

84.21091 


3 

0 

2. 

179.26881 

-.79821 

-111.58610 

636.92806 

.136487 

-3.76851 

.38781 

84.21330 


3 

10 

0 . 

179.26610 

-.79719 

-85.40157 

637.00392 

.119686 

-2.46531 

.41452 

84.20792 


3 

15 

31. 

179.26769 

-.79795 

-190.69802 

637.04229 

.111618 

-1.39404 

.41346 

84.21500 
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T = -11.102°. Day 263 data was reduced by assuming T to be unknown, and a 
result of T = -11.092° was obtained. 

D . Attitude as a Function of Time 

Table 3-VII presents a summary of the attitude parameters for each day 
in which measurements were taken. Except for Days 132, 133, 134, and 135 a 
solution was obtained for each day on which measurements were made. No 
solutions were possible on these days because of the unfavorable position 
of the sun. For these four days, a computer program could not obtain a star 
identification. 

The values shown in Table 3-VII represent average values of parameters 
for each day shown. Horizontal lines are used to separate groups of days 
over which no spacecraft maneuvers were performed. Thus, no maneuvers were 
present from Day 95 to Day 97, resulting in a relatively fixed direction for 
the spin axis. 

Note that e shows a strong tendency to increase with each maneuver. 

An explanation for that behavior is as follows: The sensor is fixed with 

respect to the figure axis of the satellite, but the spin axis of the satel- 
lite will very nearly coincide with the principal axis of greater moment of 
inertia. Hence, as the satellite looses mass, the spin axis may change 
slightly' with respect to the satellite's figure axis. This change will 
affect e and/or a. No systematic change in o with time is noted. However, 
the error in cr is probably two or three times greater than the error in e 
(listed in Table VIIl)c This results because a defines the position of the 
common intersection of the three slits, and this intersection is somewhat 
poorly defined because ( T f ’ is small .compared to 90°. 

From Day 249 to Day 277 a non zero s and s* were assumed because the 
satellite was in the earth's eclipse when the measurements were taken. 

Table 3-VH lists the computed s as always being positive and the computed s' 


100 




101 


TABLE 3 -VII 


ATTITUDE PARAMETERS AS A FUNCTION OF TIME 


Day 

5l 

?2 

• 

s 

" i rA 

S X 10 

*'* i n - "® 

s x 10 

r 

e 

<T 

O' 

6 

73e 

180.221 

.121 

575.9005 


0* 

-11.1688 

.193 

84.280 

61.277 

-89.748 

95e 

180.443 

.090 

576.1102 

1.8 

0* 

-11.0943 

.196 

84.410 

78.508 

-89.548 

9 6e 

180.441 

.089 

576.3982 

1.0 

0* 

-11.0834 

.220 

84.436 

78.548 

-89.550 

97e 

180.460 

.088 

576.2735 

1.3 

0 * 

- 11.1111 

.223 

84.369 

79.283 

-89.533 

131 

180.180 


576.8571 

0* 

0* 

-11.087* 

.199* 84.437* 



132 











133 

No solution 

was obtained for these four days because of the unfavorable 


134 

position of 

the sun. 








135' 











151 

180.829 

-.217 

558.765 

0* 

0* 

-11.101 

.230 

84.453 

104.67 

-89.142 

171 

180.708 

-.731 

597.150 

0* 

0* 

-11.087* 

.199* 

84.437* 

135.92 

-88.981 

193 

180.686 

-.617 

597.891 

0* 

o* 

-11.087* 

.199* 

84.437* 

131.97 

-89.076 

201 

180.051 

-.878 

598.923 

0* 

0* 

-11.087* 

.199* 

84.437* 

176.68 

-89.119 

24 9e 

179.550 

-.785 

636.8128 

1.8 

-5.2 

-11.102* 

.375 

84.208 

209.82 

-89.094 

250e 

179.532 

-.788 

636.7299 

2.1 

-6.4 

-11.102* 

.405 

84.199 

210.70 

-89.082 

263e 

179.403 

-.805 

637.0817 

1.1 

0* 

-11.092 

.458 

84.213 

216.56 

-88.997 

269e 

179.346 

-.805 

637.0487 

1.0 

-1.2 

-11,102* 

.378 

84.204 

219.09 

-88.962 

270e 

179.336 

-.804 

637.0136 

1.2 

1.6 

-11.102* 

.392 

84.178 

219.45 

-88.956 

271e 

179.324 

-.804 

637.7599 

1.1 

-3.7 

-11.102* 

.410 

84.181 

220.06 

-88.948 

272e 

179.315 

-.802 

636.9886 

1.3 

-3.8 

-11.102* 

.431 

84.173 

220.50 

-88.944 

273e 

179.306 

-.806 

636.9027 

1.6 

-6.1 

-11.102* 

.451 

84.179 

220.73 

-88.935 

274e 

179.293 

-.797 

637.0712 

1.0 

4.7 

-11.102* 

.433 

84.166 

221.57 

-88.934 

276e 

179.278 

-.799 

636.9226 

1.1 

-4.4 

-11.102* 

.364 

84.178 

222.10 

-88.922 

277e 

179.268 

-.798 

636.9651 

1.2 

- 2.8 

-11.102* 

.407 

84.212 

222.53 

-88.916 

278e 

179.263 

-.770 

636.9112 

0* 

0 * 

-11.094* 

.460* 

84.212* 

223.73 

-88.933 


* Jlxed value, e Indicates the satellite was eclipsed when the measurements were taken 
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as usually negative. This result is theoretically possible as may be seen from 
the following simple argument: Assume the satellite is spinning about a fixed 

axis , then 


Ttfhere Lq 
C 

but , 0 

0 

since 3 
So , 

•# 

s 

••• 

s 


= C s , s > 0 

= the satellite angular momentum 

= the moment of inertia about the spin axis 

* • •* 

= C s + C s 

= C s + 2 C s + C ‘s ' 

= 0 . 

C s 

s ^ • z. *° . 

= -r (2 c - c c) - 

c 


Now, if the satellite is losing heat due to radiation, C < 0, C > 0, Hence 
s > 0. 

But 

s’ > 0, for 2 C 2 > C C > 0 
s'< 0, for C C > 2 C 2 > 0. 


The above argument is further discussed as Appendix B 


Figure 3-9 plots, in pol’ar coordinates, the declination and right ascen- 
sion of the spin axis for each day on which results were obtained. Note that 
this axis was always less than 1.05° from the South Pole, but its right 
ascension varied by almost 180° in the time interval from Day 73 to Day 278. 
No spacecraft maneuvers were performed from Day 95 to 97, nor from Day 171 to 
278. 


E. Accuracy of Attitude Determination 

A primary reference of higher precision than the star sensor used to 
obtain attitude solutions is, of course, not available. The accuracy of the 
computed attitudes are therefore only estimates. Two methods for estimating 


102 




103 



DETERMINATION OF SATELLITE ATTITUDE 


accuracy are utilized. The first requires two or more sets of independent 
transit times to compute the attitude. The accuracy is then estimated from 
the degree of consistency in the results. The second method makes use of 
the residuals in the least squares solution of the attitude parameters to 
estimate the accuracy of those parameters. Both methods yield about the 
same estimates. 

Table 3-VII lists the estimated standard deviation of the attitude 
parameters for both good and bad geometry obtained on Days 201 and 273 
respectively. 

TABLE 3 -VIII 

STANDARD DEVIATION OF ATTITUDE PARAMETERS 
GOOD AND POOR GEOMETRIES 


Parameter 

Good Geometry (Day 269) 

Poor Geometry (Day 201) 

1 

5.5 

36.7 

1 \ 

4.8 

74.4 


16. 

37.8 


.29 

2.1 

1 

.008 

— 

1 

.0001 

— 

1 

35.7 

-- 


12,0 

— 


In Table 3-VIII the standard deviation of the angles are indicated in 
seconds of arc, and time is measured in seconds. 

Generally speaking, for all data the standard deviation of the mean 
direction of the spin axis is between 6 arc seconds for good geometry and 
100 arc seconds for poor geometry. For the azimuthal direction (rotation 
about the spin axis) the standard deviation is between 15 arc seconds and 
40 arc seconds. 
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CHAPTER 4 

DETERMINATION OF SATELLITE ORBIT 
I, NAVIGATION DATA PROCESSING 

A method of obtaining the orbit of a satellite by sightings of. a second 
satellite is described herein. Also included in the following section is a 
description of a computer program which may be used to calculate the orbit. 

A . Analysis of Navigation Problem 

The general problem is to obtain the time-path of .a satellite by sight- 
ings of a second satellite. In practice, two secondary satellites were used 
as targets; a 3 inch diameter probe ejected from the satellite itself, and 
the moon. However, these targets were never used together during a common 
time interval. Sightings from the satellite to the secondary satellite are 
supplied by three slits fixed in the ATS-III satellite. These same" slits 
were used to obtain the attitude of the satellite as discussed in Chapter 3. 

Generally speaking, the time-path of the satellite is obtained from 
observing parallax in the direction to the secondary satellite as reviewed 
from the primary satellite. This direction is a function of time and is 
defined with respect to a fixed star field. However, it was not measured 
directly by the SCNS instrument on-board the ATS-III satellite. In the fol- 
lowing paragraphs consideration is given to defining the specific nature of 
the measurements. 

1. Constraint Implied by Slit Transits 

A discussion of geometric properties of the SCNS optical instrument 
pertinent to the determination of attitude was presented in Chapter 3. For 
the ATS-III/SCNS experiment three slits were etched on the focal plane of the 
optical system. If each slit is a straight line segment and the optical 
system is distortionless, then a portion of a plane will be defined by each 
slit and the optical system which is illustrated in Figure 3-1. Each such 
plane contains the slit and the nodal point of the lens system. Given a 
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distant bright source, this source will be sensed by a detector behind the 
focal plane if, and only if, it lies on one of the planes defined by a slit 
and the optical system. Since the optical system is fixed with respect to 
the satellite, the slit-planes rotate with respect to the stars, moon, and 
ejected probe. These targets thus intermittently cross the slit-planes. With 
each crossing, a time is sensed and recorded. This time is called the transit 
time of the target. 

Let, (t) = position vector of the satellite at time t. 

R^ (t) = position vector of the secondary satellite's 
center at time t. The secondary being the 
moon or an ejected probe. 



Figure 4-1: Position of Satellite, Secondary, and 

Slit Normal at Transit Time 
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For the moment, assume that the secondary satellite is a point target. 
This assumption is valid for the case of an ejected probe if the time of 
transit is greater than about 5 minutes after ejection. However, it is not 
valid if the secondary is the moon. 


Now, when the secondary transits the slit at time t^, the vector from 
the satellite to the secondary must lie in the slit-plane. Thus, 


hj <h> ‘( R 2 - R 1 = 0 


( 1 ) 


where (t.) = unit vector normal to the jth slit-plane at t. (notation 

J 6j i i 


used in Chapter 3). 


To calculate the dot product in (1), the vectors j^and - R^ must be 

expressed in the same coordinate system. This system is chosen to be with 

in the direction of the First Point of Aries, in the direction of the 

North Pole, and i, = £, x . Since the attitude of the satellite at any 
J 1 1 1 

time is computable by use of star transits, j 
. The earlier notation of Chapter 3 yields 


6j 


(t^) has known components in 



where 


B 2j iS 3 


RCt,) is 


H is 


1x3 matrix (row matrix) 
a 3 x 3 matrix 
a 3 x 3 matrix . 


107 




DETERMINATION OF SATELLITE ORBIT 


Let B„ . R (t.) H = N. (t . ) where N. (t . ) is .1 x 3 matrix and is com- 

2j i J i j i 

putable from the transit times of known stars. Thus (1) may be rewritten in 
matrix notation as 

i 

."i < E i> (*2 (V -h (t i>)’ 0 < 2 > 

where R„ (t,) and R. (t.) are 3x1 column vectors. 

2 i 1 l 

Equation (2) is a constraint equation, for it constrains the position 
of the satellite and the secondary which may occur at the measured transit of 
the secondary satellite. This equation is used if the secondary is an ejected 
probe since the probe may be considered to be a point target at times greater 
than 5 minutes after ejection. 

To derive the constraint equation for the case in which the secondary 
is the moon, two new complications must be considered in the analysis; viz., 
the moon is an extended body, and the slit is not an ideal line, but has a 
non-zero width. A typical sensor output, as the slit crosses the moon, is 
shown in Figure 2-5. It was determined that the most clearly identifiable 
feature in the sensor output from moon crossings occurs as the slit's trailing 
edge leaves the moon's limb. 

The constraint equation obtained if the moon is used as the secondary 
satellite thus becomes 

" R m = N j (t i ' At) [ R 2 (t i " A fc > " \ & ± - At)] (3) 

where 2^t = duration of time a star is within the slit 

R^ = radius of the moon (-R^ appears in the equation because of the 
convention used to define the direction of the positive slit- 
plane normal).. 
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The slits -were constructed so that At is the same for each slit. It. is com- 
putable since the angular extent of the slits and their rate of rotation is 
known. Equation (3) holds only when a slit leaves any portion of an illum- 
inated limb (as opposed to a day-night line) . This was always the case when 
the measurements were taken. 

In summary. Equations (2) and (3) are relations which constrain the 
motion of the secondary satellite and the primary satellite. Equation (2) 
holds if the secondary is an ejected probe, while (3) holds if the secondary 
is the moon. The row vector 1L in. these equations is computable from the 
attitude which in turn is found from star transits. These equations, in 
themselves, are not sufficient to determine the orbit for the positions 
change with each measurement. However, the equations of motions of the 
primary and secondary satellite offer added information which is sufficient 
for a solution, 

2. Analytical statement of the Problem 

The analytical problem to be solved so that the time-path of the 
primary may be determined from transits of the secondary may now be stated. 
First assume the secondary is an ejected probe. Then, the equations of motion 
yield 


• « 

\ (t) = G (R k (t), t), k = 1, 2 
for all t. 

h < - t o ) " ®1 + ^ 

where 


(4) 


(5) 
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tg = time of ejection 
Av = ejection velocity 

0 = Nj <t ± ) (l 2 (t.) - \ (t.) ), j = 1, 2, 3 (6) 

where t = a measured transit time of the probe across the jth slit. 


The 1x3 row vector N^ in (6) is computable from the attitude solution. 


-The problem is to find vectors (t) and R 2 (t) which satisfy the 
following conditions : 


(a) the given differential equations (4) for all t, 

( b ) the incomplete initial condition (5), and 

(c) the constraint equations, (6) , at discrete times i = 1, 2, 
. . . , n. 


. If the moon is used as the secondary, then its time-path may be con- 
sidered as known. Hence, for this case the problem takes the form 

(R 1 (t), t) (7) 

for all t. 

0 = N. (t ± - At) (r 2 (t. - At) - R t (t. - At)) (8) 

1 - 1, 2, 3 

i =1, 2, ..., n. 

3. Method of Solution 

First, assume the secondary satellite is' an ejected probe. Also, 
suppose R 1q (t Q ) and (t^) are chosen such that they are "close to" their 

true values, R^ (t^) anc ^ (tg) . Moreover, assume the ejection velocity Av 
is given. 
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This information is sufficient to yield R^ (t) for all t as the solu- 
tion of the differential equations (4) . Now, the position vectors so obtained 
wiil undoubtedly not satisfy the constraint equation (6) at each measured 

transit time t.. Instead, define t . such that 
x ci 

0 = N j (t ci> 0*20 ^cP " R 10 ^ fc ci0 


and 


It - t I < .5 P where P is the rotational period of the satellite 

i c-x i‘ — 

(P = .6 seconds). Also, let t . - t. = fit.. Now, the 6t.’s are a measure of - 
the error in the original guesses R^ (t Q ) -and (t Q ) and R (t Q ). If each 
of these time deviations is small, then a solution of the problem has been ob- 
tained. If not, then these time deviations may be used to correct the guessed 
initial conditions by use of a linearization method. This point will now be 
considered. 


Let, 


\ < c > 


( > ct) 

\\ & 



, k - 1, 2. 


By use of this substitution, (4) may be written as a first order system of 
differential equations of the form 



F (X k , t) 


k - 1, 2 


(9) 


Given (t^), the initial position and velocity of the probe is determined 
by 


20 


- X 10 


(t 0 ) + 



) 
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where AVq = nominal velocity of probe ejection. Now, for all t, (t) can 

be found as a solution to the system of differential equations 

\ 0 = F (X kQ , t) , k'= 1, 2 (10) 

The solution vectors, (t), constitute the "nominal trajectories" of the 

satellite and probe. 

Now let X k .(t) = X kQ (t) + 6 X k (t) . 

Substituting this expression into (9) and taking advantage of the fact 
that X k Q (t) must satisfy (10), one obtains an approximation equation which 
6X k (t) must satisfy. Thus, 

5X k = VF (X k0 , t) 6X k + ... (11) 



evaluated at X k Q,(t). Here F^, F^, ... Fg are the components of F. 

So, to a first order approximation, fiX^^ satisfies a linear homogeneous dif- 
ferential equation. The coefficients of 6X, are non-constants, being 

K. 

functions of the nominal trajectory. The solution of this system may be 
written as 
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6X k (t) = § k (t ’ V 6X k (t o> 


( 12 ) 


where 

•\ (t ’ v = ™ Ko> e) \ (tj t o ) 

? k ^0’ t 0* > = I » ( 6 x 6 identity) . 

Equation (12) relates position and velocity deviations from the nominal 
trajectory at one time to position and velocity deviations at another time. 
The matrix $ in this equation is called the transition matrix in the language 
of orbital theory, it is called a fundamental solution matrix in the language 
of differential equation theory. 

Equation (6) must also 1 be linearized. To do so, define t such that 

0 = N. (t . ) (r (t .) - R (t .)') . 
j cx V 20 cx 10 ciV 

[ t^ - t_J < .5 P, where P = .6 seconds is the rotation period of the 
satellite. 

Since N, R ori , and R are known functions of time, t . is computable for 
5 20’ 10 5 ci 

each i. .Now set 


t . + fit. = t. 
cx xx 


and rewrite (6) in the form 


0 = N. (t . + 6t.) (R on (t . + 5t.) - R iri (t . + St.) + 

j - cx i' \ 20 v cx i' 10 ' cx i' 

6R 2 (t c . + 6t i ) - SR-,1 (t c .. + St.)) . 
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If higher powers of fit. are neglected, one obtains U. (t .) 6 . « 

«=ci> ( 65 2 <t ci ) _ 3 

where Ik = lk (R 20 - R-^q) + N (R 20 - R 10 ) . 


(13) 


Equation (13) relates deviations of the position of the primary and ejected 
probe from their respective nominal positions to a change in the transit 
time of the probe. In this equation all quantities are known except 

6 ^2 < t ci ) and 6 ^1 < t ci ) ’ 

Substituting (12) into (13) one finds 


U. fit . 
J i 


N.K 

J 


(t ci’ 


V 


fiX 


2 


(t 0 ) - 


«, (t 


cx ' 


‘o> 


6X- 


,<v) 


(14) 


where K = (I, 0), I being the 3x3 identity 

0 being the 3x3 zero number. 


Thus far, no use of the incomplete initial conditions (5) has been 
made. These conditions imply 


§X 2 (t 0 ) = 6 X l (t Q ) + 



where fiAv = Av - Av^ = deviation of ejection velocity from its nominal value. 
Hence, (14) becomes 


U. fit. = N. 
1 i 1 



" a V $2 “ ^ 2 ^ 



(15) 


where, (t, t^) = (o'^ 3 3^)* anc ^ 3^ are 3x3 submatrices of § . 


U., N., cv 1 s and R's are evaluated at t . . 

J J M cx 

Equation (15) represents one equation in nine unknowns. If a suffi- 
cient number of independent measurements are made, these nine unknowns may be 
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computed. After finding SX^ (t^) and fi^v, better reference trajectories may 
be obtained, and the entire process may be repeated. 


For reasons which will be considered later, the ejection speed (Avj, 
cannot be treated as an unknown in the total problem. However, the direction 
of ^v must be considered unknown, for the ejection mechanism is not' to be 
relied upon to furnish a known direction. 

In order to rewrite (15) so that jAvlis considered as known, let 

Av = fiv u 

= Av (cos e cos a i + cos e sin a j^+ sin e kj) 

where 

Av = ejection speed = j Av | 
a = right ascension of ejection velocity 
e = declination of ejection velocity. 

Also let, 

u = Uq + 6 u, Uq the nominal value of u. 


The reference orbits are chosen to satisfy 


X 2 (t 0 ) =X 1 ( V- + AV ( “C 


Hence 


5X 2 (t 0 ) = 6X X (t Q ) + Av 6u 
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Now 


6u' = 


- sm a o - sin e 0 cos a 0 
cos a D - sin e 0 sin a 0 


= M 


0 • 


cos 




cos e 0 6a 0 
6e„ 


, M a 3 x 2 matrix. 


Hence, (15) becomes 


u j (t ci ) 6 ti = N j (q, 2 - h ~ V Av ?2 M) 


/ sx, ct 0 ) \ 


cos e 0 6a 0 


\ * / 


( 16 ) 


The problem now contains the eight unknowns (t Q ), cos e 0 6a 0 , and §e o . 
This is the recommended form to be used with transits of an ejected probe. 
The ejection speed' Av must be given. 


Now, consider the secondary to be the moon. For this case t . will be 

cx 

defined such that 


R m - H j <C cl> O 20 < c ci> - R 10 «=ci>) 


It - t.j < .5 P. Define $t, as 6t. = t . - t. + At. 

ci 1 — 1 1 ci 1 

Assume the orbit of the moon is known, then (13) becomes 

u, <t cl ) St. - - H. (t c .) SR, Ct ci ) (17) 

and (16) becomes 

D J (t ci> 6t l - - far »1> (t 0 ). (18) 
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Hence, for the case in which the secondary is an ejected probe (16) is 
one member of a linear system of equations whose solution yields corrections 
to the initial guess of the spacecraft position and velocity and direction 
of ejection. If the moon is the secondary, (18) represents the correction 
equation. For this case the moon's orbit is assumed given. 

B . Data Processing Program 

The steps required to prepare and run the computer program to solve 
the navigation problem using both the moon and satellite particle are out- 
lined below in which' preparation of satellite attitude, selection and identi- 
fication of particle and moon transits and the filtering and adjusting of the 
data are described. A description of the operation of the actual navigation 
program is included later. The block diagram in Figure 4-2 shows the data 
processing tasks required to determine the satellite orbit from observations 
of the particle. Figure 4-3 shows the tasks required when using moon 
transits . 

1. Probe Navigation Data 

( 1 ) 

Ejection of the second of three particles took place on Day 202 at 
0000 GMT (July 20, 1968). Attitude data had been collected for several 
minutes prior to ejection and continued for some 15 minutes afterwards. 

Transit data was digitized live from one minute before to approximately 
three minutes after ejection and at various times thereafter for six hours. 
Star data was also gathered on the second and third nights of the navigation 
experiment but sighting of the probe could not be made from the first night 

The process of isolating probe and transits from the noise and from star 
transits can be described by the following operations outlined in Figure 4-2: 


The first particle was ejected prematurely, preventing useful naviga- 
tion sightings. The third particle is contained within the satellite 
as of the publication date of this document. 
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Figure 4-2; Data Processing Tasks for Particle Navigation 
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Average intermixed solutions to obtain 
final satellite orbit 


Figure 4-3; Data Processing Tasks for Moon Navigation 
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1. Obtain, satellite attitude for each data set. 

2. Plot a histogram of all pulses resolved to the first can. 

3. Eliminate all identifiable star pulses 

4. Find three repeatable, properly spaced transits from the 
particle . 

5. Search all scans in each set for transits close to the 
manually determined particle transits. 

6. Filter all candidate transits on the basis of both pulse 
width and position within the scan. 

The steps outlined above were repeated for each set of data from short- 
ly after particle ejection to the time where the particle could no longer be 
observed. Identification of the particle transits became more difficult as 
it moved farther from the satellite, not only because of the lower detect- 
ability with distance as one would expect, but also because the particle was 
very nearly coincident with a transit from a faint star. 

Preliminary attitudes were obtained for many of the data sets. Based 
on these attitudes, a selection was. made manually of the bright stars thought 
to be most reliably detectable. Then, attitudes were recalculated using only 
these stars which reduced random errors introduced by misidentif ications , 

Slit parameters determined with data obtained from earlier experiments were 
used. The attitudes so obtained are given in Table 4-1. 

2. Moon Navigation Data 

Collection of experimental data for the moon navigation phase of this 
project was completed during September and early October. Two groups of days 
were favorable for observing the moon while the ATS-III spacecraft passed 
through the earth's shadow each day during the fall eclipse. Data was ob- 
tained on September 5 and 6 and from September 25 through October 4, The 
time base (29 days) over which observations were made was an important con- 
sideration. During any given night, the spacecraft moved through as much as 
1/30 of "an orbit during the observations and from September 25 through 
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TABLE 4-1 

ATTITUDES FOR SECOND PARTICLE NAVIGATION (DAY 202) 


Time 

?2 

? 2 

So 

s 

0001:06 

179.98510 

-0.94688 

103.92234 

598.92117 

0005:07 

180.03736 

-0.90411 

104.01863 

598.92306 

0008:06 

180.01195 

-0.92162 

103.98668 

598.92328 

0013:36 

180.06158 

-0.89467 

103.95315 

598.92403 

0014:37 

179.98780 

-0.91140 

104.00222 

598.92292 

0147:00 

180.03745 

-0.92439 

103.77701 

598.92436 

0148:20 

180.07191 

-0.88667 

103.80030 

598.92288 

0148:52 

180.05138 

-0.90771 

103.89069 

598.92580 

0217:01 

180.08580 

-0.83968 

103.72557 

598.91330 

0217:46 

180.07667 

-0.86864 

103.76770 

598.91382 

0218:01 

180.05829 

-0.90165 

103.66706 

598.91282 

0230:00 

180.08717 

-0.84285 

103.74656 

598.90948 

0230:15 

180.06163 

-0.91177 

103.62944 . 

598.90946 

0230:30 

180.07431 

-0.84206 

103.57260 

598.90930 

0231:00 

180.06994 

-0.89183 

103.74505 

598.90954 

0231:15 

180.09868 

-.0.80564 

103.65933 

598.90911 

0231:44 

180.07402 

-0.87156 

103.74375 



598.91048 
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October 4 the moon travelled through about 9/28 of its orbit. The relative 
motions of these two objects determine the final accuracy of the navigation 
results . 

Transit data was collected successfuly on ten days during the eclipse 
with the moon in the field of view. The amount and distribution of moon 
transits is tabulated in Table 4-II. 

Preliminary attitudes for most of the data sets on each of the 10 days 
indicated in Table 4-II were computed. 

The attitude solutions obtained each day are tabulated as Table 4-III. 
Notice the slow drift in right ascension with time. Clearly the spin rate 
does not represent the sunlit rate, nor can the rates be compared directly 
from one day to the next since the eclipse does not begin at precisely the 
same time each day and the spin rate is directly related to the cool-down 
experienced while the spacecraft was in the earth's shadow. The column 
headed "Residual" gives the root mean square of the functional residuals in 
microseconds for the data set in question (approximately 15 seconds of data). 

Figure 4-4 shows the approximate position of the moon with respect to 
the star field on the observation days. It is important to note the following 
features : 


1. Although both Mercury and Venus agpear in Q the same part of 
field of view (FOV from about -10 to -40 declination), no 
evidence of any Mercury detections can be found in the 
digital data. Since the visual magnitudes of these two 
planets are +0.70 and -3.40, respectively, and because of 
the dramatic cutoff with distance from the optical axis, 
this could have been anticipated. 

2. Because of the relative geometry of the Sun-Earth, Moon, 
and spacecraft, and since the sensor scanned in a backward 
direction in right ascension, all three slits leave the 
moon from a point on its limb (in contrast to the day-night 
line) on each observation day. Since the trailing edge of 
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TABLE 4-II 

MOON TRANSIT TIME DATA 



Time 


Numbei 

• of Transits 


Day 

From 

To 

Data Sets 

Slit 1 

Slit 2 

Slit 3 

Total 

249 

0301:53 

0313:41 

1-47 

562 

718 

698 

1978 

250 

0256:52 

0301:53 

1-21 

108 

288 

297 

693 


0316:16 

0317:51* 

22-67 

227 

569 

588 

1384 

269 

0255:27 

0258:06 

1-11 

81 

141 

141 

363 


0300:00 

0303:58 

12-27 

123 

201 

192 

516 


0313:58 

0326:15 

28-76 

413 

678 

661 

1752 

270 

0246:18 

0248:08 

1- 8 

54 

107 

103 

264 


0256:01 

0258:56 

9-20 

78 

153 

156 

387 


0306:01 

0308:56 

21-32 

71 

169 

160 

400 


0316:01 

0318:56 

33-44 

74 

167 

178 

419 


0321:01 

0325:50 

45-63 

97 

217 

226 

540 

271 

0239:02 

0241:47 

1-12 

105 

160 

167 

432 


0253:11 

0255:56 

13-24 

102 

163 

180 

445 


0307:00 

0309:45 

25-36 

98 

164 

176 

438 


0318:00 

0324:47 

37-64 

229 

388 

389 

1006 

272 

0240:03 

0242:49 

1-12 

176 

182 

175 

533 


0259:01 

0301:50 

13-24 

0 

0 

0 

0 


0308:31 

0313:38 

25-45 

163 

174 

174 

• 511 


0318:03 

0324:48 

46-73 

373 

395 

379 

1147 

273 

0231:02 

0233:47 

1-12 

74 

119 

127 

320 


0246:00 

0249:00 

13-25 

64 

141 

144 

349 


0306:00 

0308:45 

26-36 

84 

130 

128 

342 


0316:01 

0323:32 

37-64 

173 

369 

380 

922 

274 

0314:47 

0321:10 

1-15 

0 

148 

147 

295 

276 

0237:48 

0241:35 

1-16 

60 

172 

164 

396 


0245:00 

0248:45 

17-32 

155 

236 

259 

650 


0256:00 

0259:55 

33-48 

113 

185 

183 

481 


0307:01 

0310:45 

49-63 

35 

203 

195 

433 


0314:45 

0317:46 

64-76 

154 

186 

194 

• 534 

277 

0255:02 

0304:58 

1-39 

235 

524 

525 

1284 


0310:01 

0320:48 

40-82 

249 

605 

588 

1442 


Totals 4303 7483 7486 19272 


* Operational filter in place, transits not included in totals. 
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TABLE 4-III 


INITIAL ATTITUDES FOR FALL ECLIPSE MOON NAVIGATION 


Day 

Time 

Rt. Ascension 

Declination 

Spin Rate 

Period 

Residual 

Stars Used 

249 

0301:53 

207.664 

-89.162 

636.7697 

0.565354 

43 

1,16,17,73,81 


0313:26 

207.427 

-89.171 

636.8937 

0.565243 

45 

1,16,17,38,81 

250 

0256:52 

211.483 

-89.079 

636.7308 

0.565388 

30 

1,16,23,38,39,48,73,81 


0301:53 

211.359 

-89.084 

636.7917 

0.565334 

29 

1,16,23,38,48,73,81 

269 

0255:27 

215.599 

-88.959 

636.9445 

0.565198 

55 

1,16,17,23,38,39,48,73,81 


0325:59 

217.056 

-88.970 

637.1557 

0.565011 

50 

1,16,17,23,38,39,48,73,81 

270 

0246:18 

219.769 

-88.960 

636.8642 

0.565270 

27 

1,16,17,23,38,39,48,73,81 


0326:44 

219.140 

-88.975 

637.1554 

0.565011 

26 

1,16,17,23,38,39,48,81 

271 

0239:02 

219.361 

-88.950 

636.7721 

0.565351 

30 

1,16,17,23,38,39,48,73 


0324:47 

220.146 

-88.952 

637.1375 

0.565027 

28 

1,16,17,23,38,48,73 

27.2 

0240:02 

220.146 

-88.944 

636.7996 

0.565327 

26 

1,17,23,38,48,73 


0324:48 

220.277 

-88.948 

637.1499 

0.565016 

28 

1,17,23,38,48,73 

273 

0231:02 

220.674 

-88.946 

636.6723 

0.565440 

39 

1,17, 23, 48, Venus 


0323:32 

217.450 

-88.933 

637.1291 

0.565035 

46 

1,17, 23, 38, 48, 73, Venus 

274 

0314:47 

220.989 

-88.933 

637.0716 

0.565086 

45 

1,16, 17, 23, 38, Venus 


0321:10 

221.408 

-88.944 

637.1112 

0.565050 

34 

1,16,17 ,23 , Venus 

276 

0237:48 

222.504 

-88.924 

636.7313 

0.565388 

27 

1,16, 17, 23, 38, 73, 81, Venus 


0317:46 

223 .403 

-89.027 

637.0720 

0.565085 

47 

1,16, 17, 81, Venus 

277 

0255:03 

222.619 

-88.925 

636.8874 

0.565249 

30 

1,16, 17, 23, 38, 48, 81, Venus 


0320:48 

222.817 

-88.908 

637.0764 

0.565081 

35 

1,16, 17, 23, 38, 39, 48, 73, 81, Venus 
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each pulse determines the transit time, determination of 
the moon's center was simplified. 

3. Because the moon is a bright object, it completely oblit- 
erates that part of the scan approximately 25 degrees on 
either side of it. This means that on certain days one or 
more key stars are lost. Those stars which are regularly 
detected (when not in the vicinity of the moon) included 
Stars 1, 16, 17, 23, 38, 39, 48, and 81. From Day 273 on, 
the updated position of Venus was added to the star catalog 
and was treated as an additional star. 

A great deal of care was used in the process of selecting all possible 
moon transit times and in deleting those which were felt to be mis identifi- 
cations or inaccurate detections. The first step in obtaining these times 
and isolating them from the star transits was to generate an attitude solu- 
tion for as many data sets as possible on each day through the use of the 
program SCAT1. Among other information, this program lists all the repeat- 
able pulses not associated with any known star resolved by a multiple of the 
scan period to the first scan. Comparing these averages with the known posi- 
tion of the moon (from an ephemeris), the relative position within a scan 
could be estimated and the three transits punched on a tabulating card along 
with the data set number and spin period. These cards were used as input 
data along with the original digital tape to a second computer program to 
search for transits falling within a relatively large tolerance of these av- 
eraged times. Transits adjusted to a. trailing edge for each scan and their 
associated pulse widths are tabulated as Figure 4-5 for a typical data set 
from one particular day“(Day 270). The row labeled "Difference" is the 
slit-to-slit difference in transit times for A^ = slit^ - slit^ and A 2 = 
slit^ - slit^. The last entry in this row is the difference of the differ- 
ences, A = A 2 ”'A^- Notice that A is quite different from zero as would be 
expected since all three slits leave the moon from slightly different points 
on the limb. It is clear that slit 1 is not as regularly detectable as either 
slits 2 or 3 and its pulse width (shown in microseconds) varies considerably 
more than. the other two. 
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Data Set 270. 35 Period .565060 


Scan 

Slit 

One 

Slit 

Two 

Slit 

Three 


Time 

Width 

Time 

Width 

Time 

Width 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

.025378 

2673 

.030376 

817 

.035416 

1609 

2 

.025378 

2869 

.030374 

818 

.035413 

1610 

3 

.025377 

2616 

.030374 

820 

.035416 

1613 
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Figure 4-5: Computer Output of Moduloed Moon Transits for Day 270 
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Each day was scanned visually for data sets which have problems associ- 
ated with them. These were either corrected if possible by obtaining a better 
estimate of the location of the moon transits or the set was discarded. 

Those sets that remain were subjected to a third computer program for further 
processing. This program used an updated position (as obtained by -the average 
transit from the previous program) about which to search for- these same -moon 
transits except for the use of a smaller tolerance. The mean, and standard 
deviation of these moduloed transits were found and only those transits 
(trailing edges) within plus or' minus one standard deviation of the mean 
were retained which were then written on a separate moon transit tape along 
with the associated slit number. The result of the processing was a magnetic 
tape with moon trailing edge transit times and identifying slit number for 
each of the ten days of observation. The amount of available moon date from 
each day by slit number is shown in the last four columns of Table 4-II. 

3. Computer Program 

A FORTRAN computer program has been written and checked out for opera- 
tion on a Contral Data 6000 series machine to solve the navigation equations 
as described in Section IIIA. Actually, two separate programs exist for 
navigation; one for use with particle transits, and one for moon transits. 
Although the equations to be solved are very similar, the special require- 
ments of the moon navigation problem (such as the availability of an ephemeris 
for the sun and moon) made it necessary to create two different programs. 

For the purpose of this report, however, only the moon navigation program 
will be described here. The particle navigation program can be thought of 
as a subset of instructions from the moon program. 

(a) Input Parameters 

Some of the more important parameters to be preset or read in prior to 
operation of the navigation program are listed below in the order in which 
they are encountered: 

The earth gravitational field constants including fourth order 
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( 2 ) 

harmonics as given by Kozai^ . 

Gravitational constant in kilometers cubed per second squared. 

Earth and moon radii. 

Slit corrections for each observation day. 

Day and time zero (time of initial conditions) . 

Hour angle of 'the first point of Aries at time zero. 

Radiation pressure constants. 

Convergence criteria for six unknowns. 

i 

Moon ephemeris for the first six hours of each observation day 

and for 0000 hours only on all other days of interest. 

Moon horizontal parallax every 12 hours on the observation days 

and for 0000 hours only on all other days. 

Sun ephemeris at 0000 hours of each day. 

Sets of attitude parameters (eight unknowns). 

Initial estimate of the satellite position and velocity at time 

zero . 

Moon transit times and identifying slit numbers. 

(b) Program Operation 

The operation of the moon navigation program follows the lines of the 
flow chart in Figure 4-6. In order to follow the flow of the program, it is 
suggested that references be made to the equations developed in Section IA. 
Since the program is essentially one solving a complicated system of non- 
linear algebraic equations (as contrasted with a program such as SCAT1, the 
star identification program, which involves a great deal of data manipulation 
and logical decision making), it is rather straight forward with very little 
looping and decision making. 


2. Earth gravity constants were taken from NASA TN D-2103 "The drift of a 
24-hour equatorial satellite due to an earth gravity field through 4th 
order" by C. A, Wagner who cites Kozai's work as a reference. 
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Enter iteration loop on initial conditions 



Figure 4-6: Moon Navigation Computer Program Flow Chart 
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Cycle through 
all transit times 



Figure 4-6: (Continued) 
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A = B(AX) 


Figure 4-6; (Continued) 
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Figure 4-6: (Continued) 
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The program begins by setting up the array of gravitational constants 
and the slit corrections for the adjustment to transits from an extended 
source (for moon data only). Moon ephemeris data is read in from cards for 
Day 249 through Day 277. Right ascension and declination are read for the 
first six hours of each observation day and for just 0000 hours on all other 
days. Similarly, the horizontal parallax of the moon is read in for 0000 
hours and 1200 hours on observation days and 0000 hours on all others. Next 
the right ascension and declination of the sun at the beginning of all the 
•above days is read. in. 


A set of satellite attitudes (eight unknowns for the moon, four unknowns 
for the particle) is read in for each of the observation days. Finally, an 
estimate of the initial conditions of the satellite (position and velocity in 
kilometers and kilometers per second) at an arbitrary time zero (ejection 
time for the particle, a time preferably midway .between Days 249 and 277 for 
the moon) are read in. 

Next, moon transits are read in from magnetic tape and a subset of 
these is selected for use in the program. Transits are selected on the fol- 
lowing basis: 


Choose the first triplet (three consecutive transits for any 
one scan) encountered, provided the slit-to-slit transit 
separation is nearly the same for slits 1-2 as it is for slits 2-3. 

This test is required because slit one transits were questionable on two 
of the observation days and cannot always be trusted. 

Choose triplets from every third data set unless no such triplet 
exists, in which case choose the first pair of transits from 
slits 2 and 3. After finding these transits, the appropriate slit 
correction is applied. 

Next, the initial Kepler arc is set up and the appropriate arrays 
are initialized. Indexing parameters are set so that we may cycle 
from time zero to the earliest time, then from time zero to the 
latest transit time. -Indexing through the transits then begins. 

As each transit is considered, a check is made to see which side 
of time zero it is on. When the first transit forward of time zero is 
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encountered, some of the basic arrays are reinitialized so that 

the required integrations are performed in a positive sense. 

The first step after any given transit is chosen is to find the appro- 
priate satellite attitude solution (i.e,, the attitude prior to or equal to 

/ 

the current time). Then the H matrix and the product A^ C are defined. 

A test is made to see if this transit falls on the current Kepler arc. 
If not, a rectification takes place which involves forming the initial con- 
ditions for the next arc by integrating over the rectification interval 
(taken to be one hour) along the current Kepler arc. This integration is 
accomplished by dividing the time interval into nine steps and applying 
Simpson's rule to each subinterval. This process requires the use of the 
forces on the satellite which are obtained through the use of the subroutine 
FORCES . 

If the transit under consideration is on the current Kepler arc, the 
program will solve for a computed transit time within one-half a satellite 
rotational period of the measured time and will use the difference between 
these two times (computed - measured) to form the equations which adjust the 
present initial conditions to the next best set. A vector D is found by 

integrating from the time at the beginning of the Kepler arc to the computed 

» 

time. The vectors N and N are defined and the scalar U is computed. Next 
the transformation matrix f is defined and appropriately stored in the array 
V. The current row of the array B and vector A is then set up where A and B 
are defined as A = B AX with AX the correction vector to the current set of 
satellite initial position and velocity. This row in both A and B is then 
normalized by dividing by .the scalar U. This completes the processing for 
one transit time. 

When all n times have been considered, the generalized inverse of B 
(with dimensions' n rows and 8 columns) is found and the correction vector AX 
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is formed as 


T -IT 

AX = (B B) BA . 


These corrections are then added to the present values of the initial condi- 
tions to form the, next estimates 


X 


n+1 



+ AX . 


A test of the magnitude of the corrections 'is made and when all are smaller 
than some predetermined values, the iterative process is terminated and the 
final result printed. ■ ' 

In addition to the main program, several subroutines are required. In 
general, each subroutine serves a single, straight forward purpose suggested 
by its name . They are : 


GENINV calculates the generalized inverse of an N by M matrix 

ORBIT determines the orbital elements given a position and 

velocity vector 

PHIMAT sets up the 6x6 matrix § 

FORCES forms the forces for the satellite at a given time due 
to the earth's gravitational field, radiation pressure, 
and perturbations due to the sun and moon 

CVEC forms the integral I § ^ (s) F (s) ds 

XBAR sets up the position and velocity vector to the satellite 
given the' orbital elements 

ANOM computes the true anomaly at any given time from the 
orbital elements. 

BMAT sets up the transformation matrix (3 x 3) B 

EMAT sets up the Euler matrix (3 x 3) and its partial derivative 

with respect to the azimuth angle s. 

HMAT sets up the H matrix (3 x 3) 
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PMOON 


calculates the position and velocity vector to the moon 
at any given time from the tabulated values of the 
, ephemeris 
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II. NAVIGATION GEOMETRIES AND PERFORMANCE 
A. Ejected Particles 

Neither of the two ejected particles (probes) furnished useful naviga- 
tion information because of the short time interval over which each was 
observed; however, consideration will first be given to the geometry of the 
position of the probe as viewed from the spacecraft. 

The salient characteristics of the probe motion are obtained by assum- 
ing that the only force acting on the satellite and probe is an inverse- 
square central force field. 

In- order to most conveniently describe the vector from the spacecraft 
to the probe, R^ (t) - R^ (t) , a coordinate system is introduced which moves 
with the spacecraft (Figure 4-7). This system will be denoted as R, with unit 
vectors r, s, z, such that r = R^ (t)/jR^ (t)j , s is in the plane of the 
spacecraft so that R^ (t) * s > 0, and z = r x s. 


A 

r 


PROBE 



Figure 4-7: The Coordinate System R 
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Let the direction of ejection be determined' with respect to R at the 
time of ejection, t Q , by the angles A Q and E 0 shown in Figure 4-8 so that the 
velocity of ejection is 

Av = Av (cos E 0 cos A 0 r (t Q ) + cos E 0 sin A o s (t Q ) + sin E 0 z) . 



Figure 4-8: The Ejection Direction as Defined by A 0 

and E 0 

Also, let the generic position of the probe with respect to the satel- 
lite be given 

R^ (t) - R^ (t) = r r (t) + s s (t) + z z. 

Then, if the orbit of the spacecraft is circular and Av and t are small, it 
may be shown that 
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r = cos E 0 (cos A 0 sin v + 2 (1- cos v) sin A 0 ) 

AV / \ 

s = COS E ° \“ 2 cos A> C 1 " c °s v ) + (4 sin v - 3 v) sin A 0 J (18) 

AV 

z = — ; sin E. sin v 
v 

where v = angle from r (0) to r (t) (figure 4-8). 

Figure 4-9 is a plot of the probe's path as viewed from the moving 
system, R, for ejection such that &v = 1 mph, E 0 = 0, A 0 =0 and .90° (two 
cases). It is assumed that the orbital period of the satellite is 24 hours. 

If Aq =0, the ejection is in the radial direction. However, as the 
probe gains distance from the earth's center, its tangential velocity is 
insufficient to yield an angular rate equal to that of the satellite and it 
falls behind. After 6 hours, it achieves- its maximum radial distance from 
the satellite, approximately 3 miles. At 12 hours, its radial distance is 
zero, but its tangential distance is -15 miles. At 24 hours it returns to 
its origin position. The motion of the particle in the moving system is an 
ellipse . 

If A 0 = 90°, i.e., the ejection is in the tangential direction, then 
the ejection speed simply is added to the satellite speed to yield the 
particle -speed. This results in an orbit of the probe whose major axis 
length and period are significantly greater than that of the satellite. 

Hence, the probe falls behind the satellite a distance which increases approx- 
imately linearly with time. After 24 hours its tangential distance is -72 
miles. Since a small component of Av ■ s (0) results in a large separation 
between the spacecraft and probe for t > 12, ejection must take place so that 
Av • s (0) = 0. 


From (18) the direction to the -probe may be written 
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MILES 
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tan A = 


- 2 cos A p (1 - cos v) -f ( 4 sin v - 3 v) sin A, 
cos A 0 sin V + 2 (1 - cos v) sin A q 


TJ 

V 


sin E 


W 


y 


9 9 9 r 

(U + V ) cos E + W* 


W = sin E 0 sin v 


• (19) 


from small y, 

• A = A 0 - v + 0 (v 3 ) 

E = E 0 - ~ v 2 tan e 0 cos 2 % + 0 (v 3 ) 


•where 


A = azimuth of R^ (t) - R^ (t) with respect to R (t), 

E = elevation of R^ (t) - R^ (t) with respect to R (t). 


From (18) and (19) the following general conclusions may be drawn: 


1. The distance between the probe and spacecraft is a linear 
function of the ejection speed &v . 

2. z (t) is a function of 4v • z but is independent of A 0 . 
Moreover, z (t) is a periodic function of t. 

3. r (t) and s (t) are functions of 4v * r (t n ) and Av • s (t ) 
but are independent of "Sv ■ z. 

4. r (t) is a period function of t, but s (t) is a period 
function of t only if A 0 =0 of 180°, i.e., ejection is 
in the plane defined by the earth's center and z. 

5. The direction to the probe is independent of A v - 

6. The azimuth at any time. A, is independent of the initial 
elevation, E . 

J O 


Condition (5) is of particular interest. Recall that the proposed 
instrument measures transit times of targets, but these measurements are 
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quite similar to directional measurements. Hence, it must be inferred that 
the measurements are not effective in determining the ejection speed of the 
probe. Conversely, an error in assumed ejection speed has little effect upon 
the accuracy of the total system if the problem is formulated in terms of the 
eight unknowns (0) , e 0 , a 0 . 

Figure 4-10 gives the path of the probe as a function of time in terms 
of its right ascension, o’ , and declination, 8^. Here it is assumed that the 
satellite's orbit is circular, equatorial, and has a 24 hour period. The ele- 
vation of ejection E 0 = -45°, various ejection azimuths ar.e shown. Note that 
the path for 0 < t < 12 hours is a weak function of A 0 but a strong function 
of A q for 24 < t ^ 36 hours. 


In order to detect the probe, ejection must take place when the satel- 
lite is sunlit. However," the sun then limits the portion of the celestial 
sphere oyer which the probe and stars can be observed. In fact, an azimuthal 
regional of approximately + 100° centered at the sun's right . ascension is 
lost. This lost region is indicated in Figure 4-10. This figure shows that 
ejection must take place so that -8° < A 0 <8° or -172° < A o < 188° if the 
probe is to be in the unblocked field of view for 24 > t > 36 hours . 

Figure 4-11 and 4-12 plots the instrument magnitude as a function of time for 
Av = 3 ft/sec, -8° < A 0 < 2°, E 0 = -45°, -40°. In obtaining these figures, 
the attenuation of the optical system as a function of angle between the 
direction of the observed target and optic axis is taken into account. Also 
included is a correction for the a-priori estimate of instrument- magnitude as 
a function of probe distance. It turned out that the a-priori estimate of 
instrument magnitude is approximately 0.5 magnitude too small. 

The ejection mechanism was designed so that Av = 3 ft/sec., and E o =-45°. 
The limiting instrument is approximately 3.0; hence, if the probe is to be 
observed for t > 24 hours , then Figure 4-12 implies that - 8 ° < A o < 2 ° : 
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Figure 4-11: Magnitude of Probe as a Function of Time for Various Initial 

Azimuths, A 0 . Full Phase. E q = -40°, Av = 3 ft/sec 
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Figure 4-12: Magnitude of Probe as a Function of Time for Various Initial 

Azimuths, A 1 Full Phase. E = -45°, Av = 3 ft/sec 

O 0 
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If the probe could be ejected so that Av < 3 ft/sec., then the allowable 
interval for A<, would be increased. However, if Av is decreased, the azimuthal 
direction % of ejection is more difficult to control. Conversely, if Av is in- 
creased , the allowable interval for A Q decreases. Av = 3 ft/sec. is optimum. 
However, comparison of Figures 4-11 and 4-12 shows that E 0 = -40° generally 
yields greater probabilities of detection than does E 0 = -45°. 


In summary, the ejection mechanism was designed so that Av - 3 ft/sec., 
E o = -45°, and ejection direction was radial with respect to the satellite's 
center as shown in Figure 4-13. However, E 0 = -40° is to be preferred. In any 
event, the allowable range of azimuthal directions imposed by the present 
optical system is narrow (10°-). Neither probe was ejected within this range 
and, hence, the observation time of each was small. 

2. Navigation Performance 

The first probe was ejected with A o = 50° which is well outside the 
permissible azimuthal range for the experiment. No useful probe sightings 
were obtained. 

The second probe was observed until 2.5 hours after ejection. At this 
time it was calculated to be about 8.6 miles from the spacecraft. All attempts 
to obtain navigation information from these sighting? proved fruitless. 

Using the position of the velocity of the satellite at the time of 

ejection as supplied by Goddard Space Flight Center, along with 149 probe 
li m 

transits 1 47 (Day 202), the direction of the ejection was obtained as 

a 0 = 261.274° (right ascension) 
e o = 22.672° (declination). 

The resulting directions in terms of the azimuth and elevation with 
respect to R (0) was 
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I 


earth's center 


I 


SATELLITE ’’ 



Figure 4-13: Probe Velocity Due to Satellite Rotation and Desired- 

Ejection Component Perpendicular to Spin Axis 
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A 0 = 9.08° 

E 0 = -22.76°* 

A standard deviation of 34.4 x 10 seconds was obtained between the computed 
and measured times . 

3. Simulated Navigation Results 

Since no results were obtained using real data because each probe was 
ejected in an unfavorable direction, one is forced to rely upon simulated 
data to give some indication of the system's potential. 

Since the satellite is presently spinning faster than the nominal 100 
RPM (about 108 RPM) , it is unlikely that the third probe could be ejected in 
such a manner that it could be seen on the spacecraft's second orbit. Hope- 
fully, however, it will be ejected nearer to the desired direction than the 
first two probes. Assuming this to be the case, an error analysis for short 
arcs will be presented. The following is assumed; 

1. Equatorial, 24 hour, circular spacecraft orbit (Kepler 
force field) 

2. Ejection velocity given by 



Av = 3 ft/sec. 

With this information, a set of probe transits of three slits is con- 
structed for 0 < t < 6 hours. To each member of this set of correct times 

an error is added. Each error is chosen from a distribution which is normal 

— 6 

with a mean of zero and a standard deviation of 40 x 10 seconds. These 
erroneous times are then used to compute the initial position and velocity 
of the spacecraft along with the direction of ejection. Errors in the 
initial position and velocity were then computed as true quantity minus erron- 
eous quantity. Results are shown in Table 4-lV. 
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In Table 4-IV, E is the standard deviation of the difference between the 
input transit times and the computed transit times (computed with the converg- 
ed solution). Note that E is -less than 40 x 10 ^ seconds, as it should be. 
n is the number of transits used, and the last column is the distance between 
the probe and spacecraft at the end of the sighting duration, T, given in the 
first column. 

The greatest distance at which the second probe was sighted was 13.8 
km. Hence, it is optimistic to assume probe sighting after 4 hours during 
the first orbit. Hence, if the probe is sighted only on the first orbit, and 
no additional information is given, then 800 km is a lower bound of the 
initial position error. 

* 

Table 4-V was constructed using similar assumptions to those used for 
Table 4-IV except that a near optimum ejection direction was chosen. Sightings 
were taken for ,5 < t < 6 hours and for 24.4 < t < 26 hours. Column 1 in 
Table 4-V is the time of the last sightings used, and the distance given in 
the last column is the maximum distance between the probe and spacecraft in 
the sighting interval. Note that considerable improvement in- the position 
determination is afforded if sightings over the second orbit are obtained. 

In both Tables 4-IV and 4-V, the errors in any component of position or 
velocity appear to be correlated when looked upon as a function of the sight- 
ing duration. This seeming correlation occurs because only one list of 
erroneous times was generated. This list was then truncated as the sighting 
duration was shortened. 

It is also noted in Tables 4-IV and 4-V that the major error is in the 
direction perpendicular to the orbital plane. This conclusion is contra- 
dictory to an observer's intuition. A possible explanation may be related 
to the fact that the instrument is less efficient in determining elevation 
than azimuth. Tables 4-VI and 4-VII indicate the effect of errors in various 
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TABLE 4-IV 

OUTPUT ERROR AS A FUNCTION OF PROBE SIGHTING DURATION. SIMULATED DATA 
Aq = 9°, E 0 = -40° , Av = 3 FT/ SEC 


T(hrs) 

Position (km) 

Velocity (m/sec) 

Ejection Dir. 

£ 

(secx 10"6) 

n 

Dist. 

(km) 

6x 

6y 

6z 

6x 

6y 

6z 

6a 

6e 

6 

- 9.0 

-27.4 

- 71.6 

3.0 


3.6 

.020 

.009 

24.4 

639 

21.6 

5 

-15.3 

-16.9 

-131.0 

2.7 


8.0 

.020 

.009 

24.6 

534 

18.0 

4 

65.9 

169.3 

674.4 

-35.2 


21.5 

.021 

.004 

23.7 ■ 

429 

14.3 

3 








324 

10.7 


t = 0 is time of probe ejection 

6x = true x (0) — erroneous x (0) = error in x (0) 

6a = error in right ascension of ejection direction 
6e = error in declination of ejection direction 

£ = standard deviation of difference between input (true) times and computed times resulting 
from erroneous orbit 

n = number of transits used 

Dist, = maximum distance between probe and satellite in sighting duration 
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TABLE 4-V 

OUTPUT ERROR AS A FUNCTION OF PROBE SIGHTING DURATION. SIMULATED DATA 
A,, = -2° 19 ' , E 0 = -40°, Av = 3 FT/SEC. 


T(hrs) 

Position (km) 

Velocity (m/sec) 

Ejection Dir. 

(<*eg) 

E 

$ecxlO 

n 

Dist. 

(km) 

fix 

Sy 

fiz 

fix 

6y 

6z 

fia 


26 

- 5.6 

9.7 

46.3 

.398 

- .957 

- .560 

-.020 

9 

24.2 

810 

22.2 

25 

- 3.8 

11.4 

57.8 

.561 

- .970 

- .652 

-.020 

1 

24.4 

693 

22.2 

24.5 

.35 

- .42 

54.7 

- .049 

.003 

.060 

-.020 

-.006 

24.4 

639 

22.2 

6 

34.0 

54.1 

104.2 

2.41 

9.25 

2.75 

-.021 

-.002 

24.9 

594 

22.2 

5 

- 59.7 

66.8 

- 151.6 

- 6.75 

-12.6 

10.8 

-.020 

-.003 

24.1 

486 

18.3 

4 

-224.0 

524.1 

-1660 

*-73.5 

-94.4 

-6.10 

-.017 

-.008. 

23.3 

381 

14.4 

3 



No convergence 
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TABLE 4 -VI 


EFFECT OF VARIOUS ERRORS IN ELEVATION FOR CROSSING 
OF SLITS 1, 2 AND SLITS 1, 3* 6E (DEG) 


Elevation ' 
of Target , 
E 


sr « .1° 



6e 

- .1° 

6t^ = 40 p,sec 

6t 2 = 
40(j,sec 


5° 

0 

-.100 

' -.010 

0 

0 

-.121 




8° 

.028 


-.010 

.001 

0 

-.119 




11° 

.054 

1 

-.010 

.002 

0 

-.117 




14° 

.079 

-.096 

-.010 

.003 

0 

-.115 

l-t 

o 

o 



17° 

.103 

-.094 

-.010 

.004 

0 

-.112 

co 

3 

g 


19° 

.125 

-.092 

-.010 

.005 

0 

-.108 

o 

•H 

> 

H 

O 

o 

o 

22° 

.145 

-.089 

-.010 

.006 

0 

-.104 

0) 

ft 

<J 

rC 

Si 

4J 

c 

24° 

.162 

-.086 

-.010 

.007 

0 

-.100 

o 

4J 

a) 

5 

27° 

.178 

-.083 

-.010 

.008 

0 

-.096 

QJ 

CO 

CO 

m 

CO 

29° 

.191 

-.079 

-.010 

.008 

0 

-.092 

o 

Si 

4J 

o 

a) 

o 

<u 

32° 

.202 

-.076 

-.010 

.009 

0 

-.088 


w 

o 

Si 

m 

o 

rQ 

34° 

.212 

-.073 

-.010 

.010 

0 

-.083 

<D 

A-i 

<y 

4J 

D 

36° 

.221 

-.069 

-.010 

.011 

0 

.-079 

CO 

u 

.5 


38° 

.226 

-.066 

-.010 

.011 

0 

-.075 

o 

U 

u 

60 

<u 

ca 

60 

<0 

4>- 

O 

o 

.232 

-.063 

-.010 

.012 

0 

-.071 


13 

13 

Slits Used 

6E same 

for 1,2 or 1,3 

1,2 

1,3 

1,2 

1,3 

1,2 

1,3 


* See Volume 3, Figure 2 s for a definition of and § ; Figure 3 for r, cr, e; Figure 5 
slit numbering and orientation. 
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TABLE 4 -VI I 

EFFECT OF VARIOUS ERRORS ON AZIMUTH FOR CROSSING OF SLITS 1, 2 

AND SLITS 1, 3 6A (DEG) 






Input Errors 



Elevation 






6t_ = 

6t 0 = 


of Target, 






2 

3 


E 

6^ = .1° 

Se = .1° 


6t i = 

40 p,sec 

40 |j,sec 

40 |j,sec 

V r 

5° 

-.001 

0 

0 

0 

.012 

B| 

.012 

Errors of .1° in 

8° 

-.001 

-.005 

-.005 

0 

.012 

Wm 

.012 

these parameters 

11° 

1 

-.011 

nsi 

0 

.'012 

mm . 1 

.012 


o 

i— i 


-.016 

BE 

0 

.012 

.024 

.012 

yield negligible 

17° 

-.003 

-.021 

mm 

0 

.012 

.024 

.012 

azimuth errors 

19° 

-.004 

-.026 

-.028 

0 

.012 

.024 

.012 


22° 

-.004 

-.032 

-.033 

0 

.012 

.024 

.012 


24° 

-.005 

-.037 

-.037 

0 

.012 

.024 

.012 


27° 

-.005 

-.042 

-.044 

O’ 

.012 

.024 

.012 


29° 

-.006 

1 

-.050 

0 

.012 

.024 

.012 


32° 

-.006 


-.055 

0 

.012 

.024 

.012 


34° 

-.007 

-.058 

-.061 

0 

.012 

.024 

.012 


36° 

-.007 

-.063 

- .066 

0 

.012 

.024 

.012 


38° . 

-.008 

-.068 

-.072 

0 

.012 

.024 

.012 


o 

o 

r . 008 

-.074 

-.078 

0 

.012 

.024 

.012 


Slits Used 

1 ,2 or 1, 3 

H 

1,3 

H 


HM 


1,2 or 1,3 
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inputs upon elevation and azimuth error respectively. In constructing these 
tables it is assumed that transits occur in Slits #1 and #2, also Slits #1 
and #3 (#2 is center slit), but not in all three slits. The geometry of the 
slits and definition of symbols were discussed in Chapter 3. 

The first three columns depicting the output error, SE, show the effects 

of errors of 0.1° in r, cr, and g upon the elevation. For these parameters 

the same errors are obtained if Slits 1, 2, or 1, 3 are used. The next two 

columns show the effect of a 0.1° error in s if Slits 1, 2 or 1, 3 are used. 

The fifth column is the resulting error if an error in t 1 , the transit time 

6 ^ 

pertinent to the first slit, of 40 x 10 seconds is made and Slits 1, 2 are 
used. If Slits 1, 3 are used, the errors are multiplied by one-half as indi- 
cated in the sixth column. The last two columns give the effect of 40 x 10 
second errors in t^ and t^ respectively. 

The columns in Table 4-VImay be interpreted in a manner similar to those 
in Table 4-VII. Note that time errors produce roughly five times greater error 
in elevations than in azimuth. Moreover, the elevation is more sensitive to 
slit parameter errors than is the azimuth. 

4. Conclusions 

I 

Two probes were ejected, neither of which furnished information useful 
for navigation because of unsuitable ejection directions. Moreover, it is 
doubtful that the third probe, which is to be ejected in mid-summer of 1969, 
will yield good navigation results. This is because the spacecraft will be 
spinning faster than the nominal 100 rpm. Moreover, it is difficult to 
envision the ejection mechanism being improved with age. 

Even though the navigation system was unable to obtain its goals, the 
system is potentially useful if the ejection direction can be controlled 
within + 5°. This control would be more easily achieved if the spacecraft 
were spinning more slowly (say 50 rpm). Also more attention should be paid 
to obtaining an optical system whose response is more nearly uniform 
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throughout the field of view. The lens system used attenuated greatly the 
response from a target as this target was moved off the optical axis. The 
effect of this attenuation can be seen by comparing Figure 4-12 and Figure 4-14. 
Figure 4-12 was obtained by using the actual off-axis attenuation of the lens 
system. Figure 4-14 was obtained by assuming no off-axis attenuation. Note 
the improvement in probe detectability. 

B. Moon 

The moon was used as a secondary in the sightings listed in Table 4-VIII. 
These sightings yielded moderately accurate results. However, the accuracy 
of these results is about two orders of magnitude less than that which, can 
be obtained by earth-bound radar. The situation would probably be reversed 
for orbits about Mars if the moons of Mars, Deimos and Phobos were used as 
secondaries . 

* i 

No measurements were made 1 October 1968 because of a failure in the 
satellite-to-ground communications equipment. For the same reason, the 
measurement duration of 30 September 1968 is extremely short. 

1. Geometry of Sightings 

When the satellite is in the sunlight, approximately 190° of the azimuthal 
field of view is lost for all targets, including the moon. Also, the moon is 
sufficiently bright so that 50° of the azimuthal field is lost for stellar 
targets if the moon is in the field of view. Hence, if the satellite is sunlit, 
and the moon is in the field of view, then only 120° of the azimuthal field 
is available for stellar targets . A field this small is generally insufficient 
for an accurate attitude determination since the spacecraft attitude must be 
found before the orbit determination can be accomplished. The moon can be 
used as a secondary only if the satellite is in the earth's shadow. This 
was the case for all sighting times shown in Table 4-VIII. 

Figure 4-15 gives the declination and right ascension of the moon for 
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TABLE 4 -VIII 

TIMES AT WHICH SIGHTINGS OF THE MOON WERE MADE 


Day 

Initial Time 

Final Time 

Duration 

5 Sept. 68 

3.03 hr. 

3.22 hr. 

0.19 hr. 

6 Sept. 68 

2.95 hr . 

3.03 hr . 

.08 hr . 

25 Sept. 68 

2.93 hr. 

3.43 hr. 

- .50 hr. 

26 Sept. 68 

2.78 hr. 

3.43 hr. 

.55 hr. 

27 Sept. 68 

2.66 hr. 

3.41 hr. 

.74 hr. 

28 Sept . 68 

2.68 hr . 

3 .41 hr . 

.72 hr. 

29 Sept. 68 

2.54 hr. 

3.39 hr. 

- . 85 hr . 

30 Sept. 68 

3.26 hr. 

3.28 hr. 

.02 hr . 

2 Oct. 68 

2.63 hr. 

3.29 hr. 

.66 hr. 

3 Oct. 68 

2.95 hr. 

3.35 hr. 

.40 hr . 
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Figure 4-15: Path of Moon with Respect to Instrument Field of View 
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each of the 10 days on which sightings were made. No sightings were made 
from 7 September to 24 September because the moon was not within the instru- 
ment's field of view. 

Figure 4-16 is a plot of the path of the moon and the orbit of the 
satellite as projected onto the plane of the moon's orbit. Also shown in 
this figure are the in-plane positions of the moon and satellite at the 
times of the sightings. Note that all sightings are made when the satellite 
is at approximately the same position with respect to inertial space. This 
occurs because all sightings were taken when the satellite was in the earth.'s 
shadow . 

2. Correction to Transit Times 

The photomultiplier output electronics are designed for optimum opera- 
tion for stellar targets. Hence, these electronics are driven beyond the 
linear range of operation by a crossing of the moon. This results in- a meas- 
ured transit time which is dependent on the moon's intensity as viewed from 
the spacecraft. Therefore, the measured transit times must be corrected as a 
function of the moon's phase. In an actual application of the technique, this 
correction function must be obtained a-priori. This was not done. Instead, 
the positions of the moon and spacecraft were obtained from Goddard. From 
these positions a set of transit times was computed. This set was compared 
with the measured times to obtain the required corrections. A rationale was 
then given for these corrections and a correction was generated for each slit 
and each day. A more complete discussion of this point is contained in 
Chapter 2. It is mentioned here so that it is clear that a measure of hedg- 
ing is present when these corrections are used. Results will also be shown 
for cases ‘in which no corrections were used. In Table 4-IX these corrections 
are listed. 

On September 5 and 6 the transits obtained for Slit 1 could not be 
used because of an improper trigger level setting. 

>i 
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TABLE 4-IX 

TIME CORRECTIONS TO MOON TRANSIT TIMES 
AS A FUNCTION OF DAY AND SLIT NUMBER 




Time Correction 

— 6 

(Secs . x 10 ) 

Day 

Slit 1 

Slit 2 

■ Slit 3 

5 Sept . 68 

Could not be used 

101 

75 

6 Sept. 68 

these days 

150 

150 

25 Sept. 68 

-90 

-72 

-43 

26 Sept. 68 

- 2 

-30 

-69 

27 Sept. 68 

15 

-56 

-60 

28 Sept. 68 

11 

-62 

-71 

29 Sept. 68 

- 4 

-62 

-73 

30 Sept. 68 

0 

-20 

-68 

2 Oct. 68 

75 

50 

17 

3 Oct. 68 

55 

87 

98 
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Figure 4-16: Moon's Position and Component of Satellite's Position in Plane of Moon's Orbit 

at the Times of Sightings 
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3. Position of Spacecraft as Obtained by Moon Transits 

In this discussion, the moon's position is considered known (obtained 
from Ephemeris) and the spacecraft's position is to be determined. The moon 
was found to be a more useful target than the probe during the experiment, in 
that it was observable over a sufficiently long time interval to yield useful 
navigation data. 

One point concerning the moon's use was disappointing. Data obtained 
on 5 and 6 September could not be used. First, attempts were made to use data 
from all 10 days and thus obtain a solution valid on the 28 day data gathering 
interval. Convergence to a solution could not be obtained. The difficulty 
occurs because of the 19 day interval from 6 September to the next day of 
data gathering, 25 September (refer to Table 4-IX) . 

The linearization of the differential equations of motion over this 
long interval introduces errors of sufficient magnitude to cause nonconverg- 
ence. Attempts were also made to use only the 5 and 6 September data to 
obtain a solution and again convergence was not obtained.' This latter diffi- 
culty occurs because of the poor geometry defined by only these two days of 
sighting (Figure 4-16). 

Results. were obtained for data' gathered from 25 September to 3 October. 
This data was used to obtain solutions valid over the 8 day interval . 

Figures 4-17 and 4-18 give the components of spacecraft position error 
obtained from one set of 722 corrected moon transits. Define 

5x = x component of spacecraft position obtained from Goddard 

Spacecraft Ephemeris minus x component computed from moon 
transits' 

with similar definitions for 6y and 6z. Also, in these figures are the 
periods over which the' observations were made. Note that they comprise an 


164 


DEVIATIONS (km) 



TIME * 


Figure 4-17: Error in x and y Components of Spacecraft's Position 

as a Function of Time. Moon Used as Target 
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extremely small portion of the 8 day interval. 

Figure 4-16 shows that -6 <; 6x, §y < 58 km in the 8 day interval. How- 
ever, these errors are significantly smaller during the times observations 
were made, in that over these times 9 < fix, 6y < 27 km. 

/ 2 2 2 

Figure 4-17 gives 6z and V (§x) + (5y) + (fiz) . Note that the error 
in 5z dominates the other components of position error to the point that 

1 6z | ^VCSx) 2 + (6y) 2 + (6z) 2 . 

The fact that 6z is the largest of the component errors is not too surprising 
in view of the geometry shown in Figure 4-16. Note that all the sightings were 
taken at approximately the same point along the orbital path. Hence, the 
orbital path may be rotated about a line joining this point and the earth’s 
center without significantly affecting the transit times. Also note from 
Figure 4-18 that 6z is small over periods at which observations were made and 
also mid-way between these periods. 

Since transits were obtained on every revolution of the spacecraft 
during the observation periods, a sufficient number of transits were obtained 
so that three independent sets of transit times could be constructed. Figures 
4-17 and 4-18 are based on Set 1. Figures 4-19, 4-20, and 4-21 show the 
results obtained from these three independent sets. Only the results for 
Day 269 (25 September) are shown. 

Finally, Figure 4-22 is a plot of the position errors for Day 269, 
using Set 1 transit times, but without using the corrections shown in Table 
4-VI. Note that when the correction for the moon's intensity is not applied, 
the errors increase by a factor of more than 10." 
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TIME (HOURS) 

Error in x Component of Spacecraft Position for Day 269. 
Three Distinct Sets of Transit Times 
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Figure 4-19: 
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Figure 4-20: Error in y Component of Spacecraft Position for 

Day 269. Three Distinct Sets of Transit Times 
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Figure 4-21 t Error in z Component of Spacecraft Position for Day 
269. Three Distinct Sets of Transit Times 
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Figure 4-22: 



Error in Position Components of Spacecraft for Day 269. 
Data Set. No Correction of Transit Times 


First 
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4. Position of Moon as Obtained by Moon Transits 

Because of the nature of the problem, it .is easy to interchange the 
role of the primary and secondary satellites. Results will .now be shown for 
such an interchange; i.e., the program will be formulated so that the position 
of the spacecraft is given, and the position of the moon is to be determined. 

The geometry of Figure 4-15 indicates that the transits should "tie 
down" the moon's orbit better than the spacecraft's orbit. This is indeed 
the case. The convergence is excellent, in that initial guesses, several 
1000 km from the solution led to the right solution. 

A single solution was obtained from the data gathered over the 28 day 
data gathering interval. A comparison between the solution so obtained and 
the moon ephemeris supplied by Goddard could not be made over the complete 
interval because the Goddard ephemeris was incomplete. 

Figure 4-23 is a plot of the difference between the position computed 
from the transits, and that obtained from the Goddard Moon Ephemeris. Three 
different days were chosen at which these position differences were computed. 
These days were at the beginning, middle, and end of the 28 day interval. 

Note that the difference lies between -50 and 40 kms . This difference is 
less than the difference between the Goddard Moon Ephemeris and the American 
Nautical Almanac Moon Ephemeris. 

Finally, Figure 4-24 is a plot of the difference shown in the previous 
figure, except that no correction of transit times, due to the phase of the 
moon, were made'. The position errors are increased for this case in that 
they lie between -50 and 480 km; however, they are still relatively small. 
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Figure 4-23: Position Errors of Moon for Three Different Days. Transit Time Correction Used 
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C. Ejection Mechanism Performance 

As the ejector torque accelerates, the ball has an equal and opposite 
. reaction force on the ejector, tending to keep the ball in the cup (shown as 
F bl in Fi § ure 1_7 )* When the ejector reaches its extended position at x = 4.5 
degrees, this force decreases and the centripel force of the ejector becomes 
significant. The centrifugal reaction force of the ball F b2 is directed away 
from the pivot P q . The resultant of these forces, F be makes an angle T) with 
the x axis such that 0 < 7] < 90°. As long as F^ is directed against some 
portion of the ejector the relative static equilibrium between the ball and 
cup is maintained. However, for 7] < 7] , the static equilibrium is not main- 
tained and the ball will pivot about the contact point P^. Thus, in addition 
to the satellite tangential velocity and the accelerator velocity, the ball 
possesses a rolloff velocity about the- retaining fingers. 

The amount of rolloff velocity is a function of T) , which in turn 

c 

depends on the physical condition of the fingers. Since 17 = £ + v" , 71 is 

c c 'c 

well defined by G Q . For example, if a nominal finger length is assumed and 

contact between finger and ball is established at the finger midpoint located 

then C ~ 15 degrees, was shown in Figure 1-7. If this contact is at P 
c i 

C ^ 18 degrees. Under these conditions the components of rolloff velocity 

would be 

G_ = 15° ; v .= 0.84 ft/ seel v = -1.67 ft/sec. 
c x y 

C c “ 18 °; v x = 1.46 ft/sec. v y = -2.16 ft/sec. 

Tests were conducted using the SONS spare flight model to evaluate the ability 
of the ejection mechanism to provide a desired ball velocity and launch angle 
During these test shots, made under laboratory conditions, the ejector stop 
angle and spring force were adjusted so that the resultant elevation and 
velocity approximated 45 degrees and 3 ft/sec. However, if the finger, shown 
in Figures 1-6 and 1-7, is bent out only slightly, then during ejector accel- 
eration the ball will roll out to make contact with the fingers prior to the 
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time normally required for £ = £ . For example, if the ejector is originally 

calibrated at £ and the fingers are then bent such- that £ has a larger 
^ cai c 

value, £ c , the resultant conditions are those listed in Table 4-X. 

TABLE 4-X 

BALL LAUNCH CONDITIONS' DUE TO EJECTOR FINGER DEFECTS 


^cal 

^c 

V x-roll 

V y-roll 

a 

e 

v 

resultant 

15 

20 

1.94 

2.43 

-15.0 

33.7 

4.1 ft/sec 

18 

25 

4.11 

3.55 

-16.3 * 

23.6 

5.7 ft/sec 

19 

24 

3.25 

3.24 

-16.5 

25.0 

4.7 ft/sec 


a - azimuth angle (degrees) 
e - elevation angle (degrees) 
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TABLE H 

SUMMARY OF OPERATIONAL CONDITIONS 


DAY 

DATE 

EXPERIMENT 

ACQUISITION 

INTERVAL 

CONDITIONS 

CONSTRAINT GEOMETRY 

REFERENCE 

73 

13 MAR 

PRELIMINARY AND DEBUS 

0800- 

0600 

GMT 

SPRING ECLIPSE 
FULL MOON S»I4.3« 
SATELLITE IN EARTH'S 
SHADOW 

Ox" illuminctedT^ 

EARTH APPEARS AS ATHIN CRESCENT 

FIGURE 2 

98 

96' 

9T 

C, 

4 APR 
8 APR 
6 APR 

ATTITUDE DETERMINATION 

0800- 

0600 

GMT 

SPRING ECLIPSE 
MOON AT FIRST QUARTER 
8 a 27,3 0 -28° 
SATELLITE IN EARTH’S 
SHADOW 

Moon 

3B*-39° 

. 

FIGURE 3 

131 

THRU 

138 

10 MAY 
THRU 
14 MAY 

( 1) EQUIPMENT CHECK 

(2) MOON OBSERVATIONS 

2300- 

0600 

OMT 

SATELLITE IN FULL SUNLIGHT 
MOON ENTERING THE 
FIELD OF VIEW AND FULL 
8=3° (130 
S C 26.7*(I35) 

Moon ^ 

6° f7 « ( ,3 '> # 

''-'Moon (138) 

orbital path -- * 

— Sun 

DV^^\ EARTH SEEN AT 0400 
GMT FROM SATELLITE 

FIGURE 4 

ISI 

30 MAY 

EQUIPMENT MODIFICATION 
CHECK 

(BASE FOLLOWING 
THRESHOLD DETECTION) 

0000- 

OBOO 

GMT 

SATELLITE IN FULL SUNLIGHT 
MOON NOT IN FIELD 
OF VIEW 
8*26° 

C_ 

1 ' / 
i FO.V 

ORBITAL Moon "0> 

PATH ^ 7 

O' JO ' Sm 

EoHh 

L- — J EARTH SEEN AT 

\ . -_JSi My 0400 GMT FROM 

SATELLITE 

FIGURE 8 

170 

in 

18 JUN 

19 JUN 

FIRST PROBE EJECTION 

0000- 

0300 

OMT 

SATELLITE IN FULL SUNLIGHT 
MOON NOT IN FIELO 
OF VIEW 
8 = 0* 

c. 

} 

F.OA 

^ESgi ' - y - wm*m 

jpMT.LLi ^ jail 

FIGURE 8 
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DAY 

K53 

EXPERIMENT 

mi 

CONDITIONS 

193 

11 JUL 

MOON OBSERATIONS 


SATELLITE IN FULL SUNLIGHT 
MOON IN CENTER OF 
FIELD OF VIEW 
8 s -84.33* 

801 

THRU 

804 

B 

SECOND PROBE EJECTION 
(ACTUAL FIRING ON 
808)0000 GMT} 
ATTITUDE DATA 
ACQUISITION ON 
801,804 

8300- 

OBOO 

GMT 

SATELLITE IN FULL SUNUGHT 
MOON NOT IN FIELD 
OF VIEW 
S = + I3.6° (801) 

S = +8S.66«(804) 

S5 

o • 

5 SEP 
8 SEP 

MOON 08SERATI0N 

0300- 

0400 

GMT 

AUTUMN ECLIPSE 
SATELLITE IN EARTH'S 
SHADOW 70 MINUTES 
EVERY 84 HRS 
FULL MOON IN FIELD 
OF VIEW 

883 

l» SEP 

DETECTION COMPARISONS 
WITH SPRING ECLIPSE 
DATA 

PMT AGING 

0300- 

GMT 

AUTUMN ECLIPSE 
MOON NOT IN FIELD 
OF VIEW 







888 

THRU 

878 


84 SEP 
THRU 
4 OCT 


MOON NAVIGATION EXP. 


OOOO- 

0400 

OMT 


AUTUMN ECLIPSE 
MOON IN FIELD OF VIEW 
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50% 

1, Venus 

16, 17 

50% 

1, 16, 17 

46 (5) 

50% 

1, 16, 17 

. 48 (5) 

50% 

1, 16, 17 

130 (5) 

50% 

1, 16, 17 

214 (5) 

50% 

1, 16, 17 

81 (5) 

50% 

1, 16, 17 

- 

50% 

1, 16, 17 

48 (5) 

50% 

1, 16, 17 

130 (5) 

12.5% 

1, 16, 17 

23 (5) 

12.5% 

1, 16, 17 

130 (5) 

12.5% 

1, 16, 17 

214 (5) 

12.5% 

1, 16, 17 

81 (5) 

50% 

16, 48 

81, US 

50% 

None (7) 

Non. (7) 

50% 

16, 46, 61 

208, 37 


VISUAL 
f MAGNITUDE 
FAINTEST 
STAR 


-1.58 


2.54 


1.63 


1.29 


2.14 


2.14 


2.91 


3.33 


2.54 


1.29 


2.14 


2.91 


1.63 


2.91 


*3.33 


2.54 


2.84 


3.30 


2 H 6 50% 16, 48 81 2.54 
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OPERATING STATE 

STARS DETECTED 

VISUAL 

MAGNITUDE 

FAINTEST 

STAR 

THEORETICAL 

MAGNITUDE 

NOISE 

FREQUENCY 

(PPS) 

STAR FIELD 
GEOMETRY 

SENSOR 

OUTPUT 

NOTES 


DATE 

SYSTEM 

VOLTAGE 


THRES 




7-11 

2 

H 

6 

50% 

16, 17, 81 
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43.0 

Fig. 7 
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16, 17, 48 

37, 130, 95(9) 

2.91 

4.45 

32.2 
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H 

6 
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1, 16, 17, 23, 
* 39, 39 

73 

2.43 

4.45 

- 
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n 

H 

4 

50% 


81, 190 

3.21 

2.95 

12.1 

Flg.10 

Fig. 15 
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5 

50% 
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115(5). 73(9) 
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09, 125, 214 

3,33 

3.7 

34.5 

Fig. 10 

Fig, 16 
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1,16, 17, 23, 
38, 39, 48, 81, 
73(91. 102(9), 
115(9), 12S(S), 
173(9), 190(9), 
214(9), 17 

73, 247, 
319, 373 

3,82 

4.45 

20.5 

Ftg.10 

Fig. 17 

(9) 

260- 

278 

9- 24- 
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50% 

Venus 1, 16, 
17, 23, 30, 
39, 48, 81 
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- 

Fig. 11 

Fig. 18 
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NOTES 

(1) Regularly detected stars wero the only targets detected. 

(2) Detected more than 60% of time In one or more of slits. 

(3) Gain flxod at 50%. 

(4) No useful output Information. 

(5) Faintest star detected and Identified. Magnitude of 
faintest star Is based upon faintest star detected In 
all three illts at least 25% of time. 

(6) Star Identifications made through manual reduction of 
digital data. 

(7) Due to Interference frorr. sun and moon, detection of 
sufficient bright stars to perform target Identification 
was not possible. (Refer to Figure 4) 

(8) Automatic determ Inst Ion of star identification and attitude was utilized. 

(9) Detectod more than 10% of the time in at least one slit. 
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APPENDIX C 
THE DETERMINATION OF 

SPECTRAL ENERGT DISTRIBUTIONS AND INSTRUMENT MAGNITUDES 
FOR SOME SELECTED STARS 
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I . INTRODUCTION 

Spectral energy distributions were determined using a black-body 
approximation and color index information for 96 bright stars in an annular 
field of view between declinations of -10 and -40 degrees. Instrument 
magnitudes for these stars are obtained using the black-body curves and 
the known spectral response of the SCNS-ATS III system. 


Using a black-body distribution. 


I(X, T) 


2nh c ^ 

^(e hC/ “-U 


k = 1.38 x 10 ^ erg/°K 
-27 

h = 6.62 x 10 erg sec. 
c = 3 x 10^ cm/sec. 

T = temperature (°K) 

X = wavelength (cm). 


and blue (B) and visual (V) response curves the smoothed fluxes of radiation, 
F^ v and F-^, and given by 

p xv ' / v x • T > 

and 

F^ b = J B-^ . I(A, T) dA. 

respectively, with the integrations over the effective bandwidths of the B 

and V curves. Generation of F^ and F values for various temperatures (T) 

and use of the relationship between color index (Cl) and flux, Allen [l] 

?-■ • 1 "if , ' 

Cl — -2.5 log <E\ b /F X v ) + .71, 
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enables the determination of a temperature, T qJ for which the color index is 
zero. If a linear dependence of color index upon 1/T is assumed then the 
slope of the curve may be obtained by plotting the quantity 

-2.5 log <F XE /P W ) 

as a function of 1/T. Using the slope and the value of T q yields a relation 
between temperature and color index of the form 

T = 6533/CI + .51). 

Black-body temperatures may now be determined for the stars solely 
from color index information. A characteristic constant, C, can be obtained 
for each star by solving for G in the equation 

M v = -2.5 log (C J * I (A, T) d\), 

where is the visual magnitude of the star. Instrument magnitudes, M^, 
are given by 

M I = -2.5 log (C J *I(A , T) d\) 
where R. is the system spectral response. 
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II. CALCULATION OF BLACK-BODY TEMPERATURES , ENERGY DISTRIBUTION AND 
■ INSTRUMENT MAGNITUDES 

The spectral energy distribution of a star may be approximated by a 

black-body distribution (Planck function) in the region from 3800 R to 
o 

9000 A. The Planck function 


!(*, T) 


2rrh c^ 

, 5 , hc/kA.T 
X (e 


( 1 ) 


where • 

k = 1.38 x 10 -16 erg/°k 
-27 

h = 6,62 x 10 erg sec. 
c = 3 x 10^ cm/sec. 

T = temperature (°K) 

X = wavelength (cm) , 


requires an estimate of the black-body temperature T. Since only color 
index and visual magnitude information are available a relationship between 
color index and temperature is necessary. 


The color index of a star is dependent upon the fluxes of radiation 
in the blue and visual regions of the spectrum. If a black-body energy 
distribution is assumed and the blue and visual fluxes can be estimated, 
then the correct black-body temperature will make the measured and 
estimated color index agree. 

m 

The Arizona -Tonantzint la Catalogue , from which the color index and 

visual magnitude information for the stars was obtained, measures color 

r 3l 

index and visual magnitude with the Johnson and Morgan J UBV photometric 
system. Blue (B) and visual (V) response curves for this system are obtained 
by plotting the values as given by Allen with those from another source 
(unknown) . All values are normalized to one at the peak. A smooth curve 
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o 

is drawn through these points and response values are noted at 50 A 
intervals over the range of response (Figure C-l). 


The smoothed blue and visual fluxes for a star characterized by a 
black-body energy distribution can now be obtained by numerically integrating 
the following equations if the temperature, T, is known. 




F xv ' J v x • 


( 2 ) 


( 3 ) 


An approximate relation between color index and temperature is (Allen) 
B - V ~ 7300°/T - .52. (4) 


If this equation is solved for T the resulting function is as shown in curve 

[41 

A in Figure C-2. Black-body temperatures as given by Arp J are also shown 
(curve B) . The values agree well for small color indices but the disparity 
becomes greater as the color index increases. 


Using the color index values obtained from the Arizona-Tonantzintla 
Catalogue and Equation (4), black-body temperatures were estimated for 
several stars. Blue and visual fluxes were then calculated using Equations 
(2) and (3) and a two step Simpson's method of integration. 

A constant, C, was then determined for each star such that the equation, 

H v =-2.5 log (C J V x * I(X, T) dX), (5) 


191 




ARBITRARY UNITS 



3500 4000 4500 5000 5500 6000 6500 7000 
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Figure C-l: Johnson and Morgan Blue and Visual Response Curves 
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Figure C-2: Color Index-Temperature Relationships 
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was satisfied,' where is the measured visual magnitude of the star. A 
second constant, C 1 , may also be computed for each star satisfying the 
equation. 


M 


B 


-2.5 log 



I(X, T) d\). 


( 6 ) 


For each star C 1 differs from G by a factor proportional to the peak value 

of the unnormalized blue response curve. The ratio, C'/C, is thus the same 

for all stars and needs only to be determined once. Using the star Sirius 

(B - V = .01, M = -1.45) the ratio, C'/C, was computed as C . For all 
v s 

stars Equation (6) now can, be written as 

M b = -2.5 log (C s • C J B x • I(\, T) dX) (7) 

where the value of C is obtained from Equation (5) for each star. 


Blue magnitudes generated in the above manner are given in Table G-I 
for 32 stars in the field of view. Also given are measured blue magnitudes 
obtained from the Arizona-Tonantzintla Catalogue and the differences between 
the computed and measured values. The differences were found to vary approx- 
imately linearly with color index as shown in Figure C-3. To correct for the 
difference dependence upon color index a new relation between color index and 
temperature was determined. 

A relation between color index and blue and visual flux (Allen) 

B - V = 2.5 log (F XV /F XB ) + .71. 

By computing the quantity G, 

G = 2.5 log ( F ^ v / F j lB )s 
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TABLE C-X 

COMPARISON OF MEASURED AND COMPUTED BLUE MAGNITUDE ^ 



(1) Star numbers are obtained by consecutively numbering all stars of visual magnitude 
less than 4.0 ordered by visual (SAO) magnitude. 
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as a function of black-body temperature a value of T, T , for which the color 
index is zero may be determined by plotting T versus G as in Figure C-5. If 
the value of T for which G = -.71 is T q then from Figure C-4, T q is deter- 
mined to be about 12750°K.. If a linear dependence of color index upon tem- 
perature of the form, 

B - V = | + C, 

is assumedj the slope. A, can be obtained from a plot of G versus 1/T (Figure 

C-5. The constant A is evaluated to be 6533. C is thus equal to-A/T or 

o 

-.51. The new relation is therefore 


B - V 


6533 

T 


-.51. 


( 8 ) 


Solving Equation (8) for T yields 
T = 6533/ (Cl + .51). 


( 9 ) 


A plot of Equation (9) is shown as curve C in Figure C-2. These values 
agree more closely with Arp's values than do those of Allen. 

Blue magnitudes were again computed using temperatures obtained from 
Equation (9), the technique employed being the same as before. The differ- 
ences between the computed blue magnitudes and the measured blue magnitudes 
are shown in Figure C-6. The blue magnitudes and differences are tabulated 
in Table c-I. The values have a constant added to them so that' the mean 
difference is approximately zero. The linear trend has been removed and 
the magnitudes of the differences are considerably less than those previously 
obtained. The discrepancy might be further improved by including higher 
order terms in the approximation of Equation (8). 

For the present, the temperature-color index relation is significantly 
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Figure C-5: Determination, of Slope of G as a Function of Reciprocal Temperature 
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Figure C-6: Measured Minus Computed Blue Magnitude as a Function of Color Index (Equation 9) 
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accurate since' the differences vary by less than + .02 units of magnitude 
from the measured values. The black-body temperatures and constants, C, 
can now be computed for all stars of interest in the field of view. 
Instrument magnitudes are then calculated as follows: 


The spectral response of the SCNS/ATS-III instrument photo- 
multiplier is shown in Figure C-7. Instrument spectral 

o 

response values, R^, are determined at 50 A intervals. The 
response of the optical glass is approximated by a step 

„ o o 

function, rising at 3800 A and constant to beyond 9000 A. 
The instrument magnitudes are then obtained by numerically 
integrating the equation. 


= -2.5 log (C J • I(X, T) d\) - Z 


( 10 ) 


\ 1 = 3800 
\ 2 = 8700 


where Z is the instrument magnitude of a zero color index, 
zero visual magnitude star'. The instrument magnitudes are 
therefore relative. Instrument magnitudes and visual mag- 
nitudes taken from the Smithsonian Astrophysical Observatory 
(SA0) star catalogue [5] for the 96 stars in the field of 
view are shown in Table C-II. 
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Figure C-7: Spectral Response of the SONS Instrument (Photomultiplier and Optical System) 
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TABLE .C-II 

INSTRUMENT MAGNITUDES FOR STARS IN THE FIELD OF VIEW 


No. 

Mag (SAO) 

Ins. Mag. 

No. 

Mag (SAO) 

Ins. Mag. 

1 

-1.58 

-1.45 

203 

3.29 

3.83 

15 

1.21 

.80 

207 

3.30 

3.22 

16 

1.22 

1.54 

208 

3.30 

3.10 

17 

1.29 

1.21 

211 

3.32 

3.71 

23 

1.63 

1.36 

214 

3.33 

3.49 

26 

1.71 

1.50 

221 

3.37 

3.13 

37 

1.95 

1.83 

227 

3.41 

3.91 

38 

1.98 

2.21 

229 

3.42 

3.91 

39 

1.99 

1.82 

243 

3.47 

3.96 

48 

2.14 

1.93 

247 

3.48 

3.86 

57 

2.24 

2.54 

254 

3.51 

3.30 

61 

2.26 

2.58 

269 

3.59 

4.20 

66 

2.36 

2.87 

271 

3.59 

3.48 

73 

2.43 

2.36 

274 

3.60 

4.02 

78 

2.51 

2.28 

275 

3.61 

4.08 

81 

2.54 

2.26 

285 

3.64 

, 3.42 

89 

2.63 

2.45 

287 

3.64 

3.68 

92 

2.69 

2.71 

289 

3.64 

3.24 

94 

2.70 

2.58 

295 

3.65 

3.93 

95 

2.71 

2.65 

300 

3.67 

3.61 

98 

2.74 

3.35 

307 

3.68 

4.08 

100 

2.75 

2.53 

313 

3.70 

3.56 

102 

2.78 

2.50 

319 

3.72 

4.06 

107 

2.80 

2.56 

322 

3.72 

4.25 

114 

2.84 

3.31 

343 

3.77 

4.11 

115 

2.84 

3.14 

344 

3.77 

3.90 

119 

2.88 

3.10 

348 

3.78 

3.80 

124 

2.90- 

2.84 

352 

3.78 

4.19 

125 

2.90 

2.50 

356 

3.80 

3.54 

130 

2.91 

2.66 

359 

3.80 

3.88 

' 132 

2.91 

2.74 

360 

3.80 

3.89 

134 

2.94 

3.38 

362 

3.80 • 

4.25 

141 

2.96 

3.28 

373 

3.82 

4.14 

144 

2.98 

3.00 

375 

3.83 

4.15 

147 

3.00 

2.80 

380 

3.83 

3.73 ■ 

152 

3.02 

3.09 

393 

3.86 

4.28 

161 

3.07 

3.53 

398 

3.88 

4.34 

166 

3.09 

2.87 

404 

3.90 

4.33 

169 

3.10 

2.93 

409 

3.90 

4.32 

171 

3.10 

2.98 

415 

3.92 

4.30 

173 

3.11 

2.91 

416 

3.92 

4.44 

175 

3.12 

2.94 

429 

3.95 

4.35 

182 

3.16 

2.93 

430 

3.95 

4.18 

183 

3.16 

3.75 

434 

3.95 

4.08 

187 

3.19 

3.59 

441 

3.96 

4.47 

190 

3.21 

3.61 

444 

3.98 

4.34 

193 

3.22 

3.69 

446 

3.98 

4.42 

197 

3.25 

3.54 
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3.25 

3.80 
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' ATS -III/SCNS SLIT CHARACTERIZATION RESULTS 
I . INTRODUCTION 

In order to optimize the performance of the operational computer 
program for attitude determination, the characterization of the SCNS instru- 
ment and its relation to the satellite angular momentum vector must be as 
precise as possible. 

From the data reduction results covering the various periods of SCNS 
operation, it is clear that a single time invariant set of slit parameters is 
not sufficient to characterize this relation. If a fixed set of slit para- 
meters is employed, it is found that any set which yields essentially zero 
average residuals over all transits from each individual star/slit combina- 
tion used in the attitude solution for one day will not give similar zero 
averages with data from another day. What this means, in effect, is that 
the known pattern- of stars and the assumed relation between the spin axis 
(total angular momentum vector in this application) and the projection of 
the slits on the celestial sphere cannot be made to "fit" — to be internally 
consistent. The most plausible "explanation" for this apparent discrepancy 
is that the spin axis of the vehicle has undergone slight shifts with respect 
to the gross geometrical or figure axis of the spacecraft. A small asym- 
metric mass displacement could account for such a change in the principal 
axis relative to the apparent geometric symmetry of the satellite. It is, 
of course, also possible that the geometrical figure of the spacecraft itself 
could warp; this, however, seems much less likely. In any event, it is 
possible for the parameters which characterize the projection of the slits 
on the celestial sphere to vary. This appendix will review the slit projection 
and document the results obtained from operational slit parameter investiga- 
tions. 
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II. GEDMETKT'OF THE SLIT PROJECTION 

The basic geometry of the slit projection on the celestial sphere and 
parameters which characterize this projection are illustrated in Figure D-l. 
The three slit parameters are; 


r = the slit-to-slit angle as seen from the projected slit vertex 
(measured on the celestial sphere) , 

cr = the angle from the spin axis to the projected slit vertex as 
seen from the center of the celestial sphere, 

e = the angle between the center slit (2) and the great circle 
arc from the spin axis to the projected slit vertex as seen 
from that vertex (measured on the celestial sphere) . 



Figure D-l; Slit Configuration as Projected on the Celestial Sphere 
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It is important to note here that a shift in the spin axis relative to 
the figure axis can change e, cr, or both. The slit separation F, on the 
other hand, is independent of the spin axis. Indeed, r depends only on the 
manner in which the slits are projected from the image or reticle plane. onto 
the segment of the celestial sphere covered by the optical field. Unless 
the SONS instrument undergoes significant internal displacements, r should 
be essentially time invariant. 

Another observation may be made from reference to Figure A-l and to 
Figure D-2 --a SCNS star field plot on which the eight most solidly detect- 
able stars are emphasized. These eight stars have significantly higher 
detection frequencies than the remainder of the stellar targets in the field. 
Thus the eight appear to represent the best set to use for a dynamic determ- 
ination of the slit parameters directly from the operational data. 

It turns out that there are reasons for discarding Sirius and these 
will be discussed later. Nevertheless, with or without Sirius, the azimuth 
distribution of the basic star set is excellent. However, it is seen that the 
stars lie in a very compact range in declination — essentially a narrow band 
in zenith angle from the ATS-III spin axis. Thus the slit-to-slit separa- 
tions will be about the same for each of the stars in the set. Under this 
condition a and F appear to be highly interdependent; a change in cr coupled 
with a correlated change in r will leave the slit-to-slit separation in the 
center of the field almost unchanged. Turned around, this implies that a 
combination of locus in the (r, a) plane can be accurately determined but 
that the specific point along that locus is only weakly determined in a 
simultaneous solution for all three slit parameters. 

The problem of resolving F and a is not particularly critical in the 
attitude determination task since large shifts along the F, cr locus will 
propagate to very much smaller shifts in the attitude solution -- at least 
in cases with many stars and good azimuth geometry. However, when the 
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position of a single navigation target (such as the moon) must he inferred 
and when this target may cross the slits well away from the center of the 
field, the individual values of r and a become important. For this reason 
a single reliable value of T is established and used as a fixed parameter. 
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III. INDEPENDENT DETERMINATION OF T 

One approach to this task is to project the known slit separation angle 
T* in the reticle (focal) plane through the optical system to infer the angle 
T on the celestial sphere. The geometry of this projection is shown in 
Figure D-3 below; 


OPTICAL AXIS 



Figure D-3* Projection of the Slit Configuration from the 
Reticle Plane to the Celestial Sphere 
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The projected angle r depends on r', x, and f; it is determined by the 
relation 


tan r* 


sin F = 


2 f 

tan F* + 


2 ' 


,2 2 

f + x 



( 1 ) 


The design values for the reticle and optical variables are r* = 10°30', 
x = 1.0675 inches, and f = 76 mm. These yield |r| = 11?13250 as the pro- 
jected slit-to-slit separation. The nominal values of the reticle dimensions 
are all subject to fabrication tolerances and measurement uncertainties. 

Since a spare reticle was available a careful validation of the measurements 
was performed. The reticles were all formed as photocopies from the same 
master so the results obtained should apply equally well to the flight unit. 

The angle r’ was measured with the aid of a sine table. The reticle 
was fastened to the table and the tilt angle adjusted with a planar gauge 
until the first slit was parallel to the motion of a 100 power travelling 
microscope on precision ways. The height of the planar gauge was then 
measured and reduced to an angle of the sine table. This procedure was 
repeated for each of the other slits with no change in the relative alignment 
of the table and the microscope motion. Note that this method has the ad- 
vantage that motion of the microscope is required to be parallel to the slit 
over the entire slit length; this permits averaging over small fluctuations 
due to any slight random irregularities in the slit edges. The slit-to-slit 
angles from slit 1 to slit 2 and slit 2 to slit 3 were each determined and 
were found to agree within the uncertainty of the basic measurements. Thus, 
it appears that the assumption of slit symmetry is justified. The value of 
T* obtained was 10?454 with an uncertainty of less than one minute of arc. 

From this value of p' and a series of point settings with the measuring 
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microscope, four independent determinations of x were also made. These 
averaged to x = 1.075 inches with an uncertainty of approximately 0.003 inch. 

A direct nodal slide measurement of the effective focal length f of the 
Super Farron lens in combination with the actual reticle glass path was also 
performed. The average value obtained from a series of measurements was 
f = 75.42 mm. With the extremes discarded, the measurement spread 'was 0.25 
mm. 


If the set of measured values {r 1 - 10?454, x = 1.075 inches, and 
f = 75.42 mm} is inserted in Equation (1), the resulting value of T, correct- 
ed for sign is 


r = -11?1020 . 


for the projected slit-to-slit separation. The minus sign is simply a con- 
sequence of the way F is defined. Evaluation of possible errors in the F 
determination was obtained by varying each of the independent variables in 
Equation (1) one at a time. The results are as follows; 


AF‘ 

= ± 1 min (+0.0167) 

Ar = + 

Ax 

= + 0 o 01 inches 

-> Ar = + 

4f 

= + 1 mm 

-» AF = + 


0°0176, 

0?0115, 

0?016. 
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IV. SIMULTANEOUS SOLUTIONS FOR T, o, AND e FROM OPERATIONAL DATA 

For two dates on which SONS operations were conducted, the data 
acquisition conditions were ideal. On Day 96 (5 April -- spring eclipse) 
and again on Day 263 (19 September — fall eclipse) the annular region 
scanned by the SONS instrument is free of any major interference from the 
bright objects — sun, earth, or moon. On these days, the base line of the 
SONS video is unperturbed and all eight of the prime stellar targets are 
solidly detected. This data, then, represents the best available basis for 
a simultaneous dynamic solution for all three slit parameters -- r, c, and e. 

Complete nine unknown attitude/slit solutions have been obtained for 
each day and the results are listed in Table D-I, On Day 96, a relatively 
small quantity of data (five 15 second data sets) is available and this was 
all processed as a single composite set containing approximately 1,000 
identified transits. For Day 263, data was acquired over a much longer time 
span and it was partitioned and processed as five independent but inter- 
leaved composite data sets, each containing 1500 to 1800 identified transits. 
The spread in the results among these five solutions is a measure of the 
random component of the uncertainty in each unknown; these spreads are listed 
in parentheses under each corresponding average value. 

The average value of F obtained from these solutions of the nine unknown 
problem (six attitude parameters plus the three slit angles) is f = 11?0967 
where the results from Day 263 have been given a weighing factor of five 
relative to those from Day 96. This is in remarkably good agreement with 
the r obtained by projecting the reticle measurements through the optical 
system. These two independent F determinations differ by much less than the 
uncertainty inherent in either one of them. Since the reticle projection is 
the more certain of these two values, r = -11.1020 was adopted as the standard 
value. The dynamic solution is treated simply as strong corroboration. 
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TABLE D-I 

ATTITUDE /SLIT SOLUTIONS FOR SCNS/ATS-III 


Solution Type 

Day 

Number of 
Independent 
Solutions 

Pointing Direction 
«1 

Slit Parameters 
T O' e 

9 unknown with Sirius 

96 

1 • 

180.44242 0.09093 

-11.07725 

84.44870 

0.23424 

w/o Sirius 

96 

1 

180.44545 0.08541 

-11.10606 

84.38752 

0.15668 

with Sirius 

263 

5 ■ 

179.40277 -0.80276 

(0.00321) (0.00302) 

-11.09409 

(0.00811) 

84.21181 

(0.01469) 

0.45998 

(0.00576) 

w/o Sirius 

263 

5 

179.40753 -0.80552 

(0.00385) (0.00448) 

-11.10133 

(0.05106) 

84.19855 

(0.09971) 

0.38148 

(0.05106) 

8 unknown with Sirius 
(ct, e at 
fixed f) 

263 

5 

179.40288 -0.80276 

(0.00309) (0.00297) 

-11.0925 

fixed 

84.21488 

(0.01055) 

0.45991 

(0.00329) 

7 unknown with Sirius 
(Center slit, 

263 

5 

179.39640 -0.80419 

(0.00627) (0.00227) 

-- 

84.21181 

fixed 

0.47875 

(0.00962) 

e only) , c . . 

J w/o Sirius 

263 

5 

179.40287 -0.80675 

(0.00835) (0.00349) 

-- 

84.21181 

fixed 

0.38975 

(0.05508) 
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The detailed solutions of the nine unknown problem reveal the correla- 
tion between r and CT quite clearly. In Figure B-4 (r, CT ) solution points 
are plotted for each of the twelve independent reductions in the nine unknown 
problem. Also plotted are the values of cr obtained at various fixed p values 
from eight unknown reductions of the data on Day 193 (11 July) . Note that 
two independent interleaved data sets were reduced and exhibit excellent 
agreement with one another. Observe also that both lunar and solar inter- 
ference were encountered on Day 193; the number and distribution of available 
solid stellar targets is insufficient to solve for all three slit parameters 
simultaneously but the solution for (cr, e) for a fixed p goes very well. 
Finally, it is seen that the (r, cr) loci in Figure D-4 are all essentially 
straight lines on any given day and that the traces for the various days are 
parallel to one another. However, at fixed p, the zenith angle cr varies with 
time (Day) as expected. 

An additional observation of some importance may be made with reference 
to Table D-l, The angles and £ describing the pointing direction are 
found to have much smaller uncertainty than do the slit parameters. In effect, 
the attitude solutions are somewhat better than the slit solutions provided 
that the slit solutions all fall on the correct (r, cr) locus. As a point of 
interest, the internal consistency of the five independent point direction 
determinations on Day 263 is better than 15 arc seconds. 
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( degrees) 

Figure D-4 (f, a) Loci 
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